The James Webb Space Telescope (JWST) is expected to observe galaxies at $z>10$ that are presently inaccessible. Here, we use a self-consistent empirical model, the UniverseMachine, to generate mock galaxy catalogues and lightcones over the redshift range $z=0-15$. These data include realistic galaxy properties (stellar masses, star formation rates, and UV luminosities), galaxy-halo relationships, and galaxy-galaxy clustering. Mock observables are also provided for different model parameters spanning observational uncertainties at $z<10$. We predict that Cycle 1 JWST surveys will very likely detect galaxies with $M_*>10^7 M_\odot$ and/or $M_{1500}<-17$ out to at least $z\sim 13.5$. Number density uncertainties at $z>12$ expand dramatically, so efforts to detect $z>12$ galaxies will provide the most valuable constraints on galaxy formation models. The faint-end slopes of the stellar mass/luminosity functions at a given mass/luminosity threshold steepen as redshift increases. This is because observable galaxies are hosted by haloes in the exponentially falling regime of the halo mass function at high redshifts. Hence, these faint-end slopes are robustly predicted to become shallower below current observable limits ($M_\ast < 10^7M_\odot$ or $M_\mathrm{1500}>-17$). For reionization models, extrapolating luminosity functions with a constant faint-end slope from $M_{1500}=-17$ down to $M_{1500}=-12$ gives the most reasonable upper limit for the total UV luminosity and cosmic star formation rate up to $z\sim 12$. We compare to three other empirical models and one semi-analytic model, showing that the range of predicted observables from our approach encompasses predictions from other techniques. Public catalogues and lightcones for common fields are available online.

We present MINESweeper, a tool to measure stellar parameters by jointly fitting observed spectra and broadband photometry to model isochrones and spectral libraries. This approach enables the measurement of spectrophotometric distances, in addition to stellar parameters such as Teff, log(g), [Fe/H], [a/Fe], and radial velocity. MINESweeper employs a Bayesian framework and can easily incorporate a variety of priors, including Gaia parallaxes. Mock data are fit in order to demonstrate how the precision of derived parameters depends on evolutionary phase and SNR. We then fit a selection of data in order to validate the model outputs. Fits to a variety of benchmark stars including Procyon, Arcturus, and the Sun result in derived stellar parameters that are in good agreement with the literature. We then fit combined spectra and photometry of stars in the open and globular clusters M92, M13, M3, M107, M71, and M67. Derived distances, [Fe/H], [a/Fe], and log(g)-Teff, relations are in overall good agreement with literature values, although there are trends between metallicity and log(g), within clusters that point to systematic uncertainties at the ~0.1 dex level. Finally, we fit a large sample of stars from the H3 Spectroscopic Survey in which high quality Gaia parallaxes are also available. These stars are fit without the Gaia parallaxes so that the geometric parallaxes can serve as an independent test of the spectrophotometric distances. Comparison between the two reveals good agreement within their formal uncertainties after accounting for the Gaia zero point uncertainties.

We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time flexibility in the assumed stellar physics. We include a model for nebular line emission that can employ either precomputed Cloudy lookup tables (for efficiency), or direct photoionization calculations for all young stars (for flexibility). The dust content follows either observationally-motivated prescriptions, direct modeling from galaxy formation simulations, or a novel approach that includes the dust content via learning-based algorithms from the SIMBA cosmological galaxy formation simulation. AGN can additionally be included via a range of prescriptions. The output of these models are broadband SEDs, as well as filter-convolved images. Powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models, and seamlessly interfaces with GIZMO, AREPO, GASOLINE, CHANGA, and ENZO. We demonstrate the capabilities of the code via three applications: a model for the star formation rate (SFR) - infrared luminosity relation in galaxies (including the impact of AGN); the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs; and the impact of galaxy inclination angle on dust attenuation laws.

In the $\Lambda$CDM paradigm the Galactic stellar halo is predicted to harbor the accreted debris of smaller systems. To identify these systems, the H3 Spectroscopic Survey, combined with $Gaia$, is gathering 6D phase-space and chemical information in the distant Galaxy. Here we present a comprehensive inventory of structure within 50 kpc from the Galactic center using a sample of 5684 giants at $|b|>40^{\circ}$ and $|Z|>2$ kpc. We identify known structures including the high-$\alpha$ disk, the in-situ halo (disk stars heated to eccentric orbits), Sagittarius (Sgr), $Gaia$-Sausage-Enceladus (GSE), the Helmi Streams, Sequoia, and Thamnos. Additionally, we identify the following new structures: (i) Aleph ([Fe/H]$=-0.5$), a low eccentricity structure that rises a surprising 10 kpc off the plane, (ii, iii) Arjuna ([Fe/H]$=-1.2$) and I'itoi ([Fe/H]$<-2$), which comprise the high-energy retrograde halo along with Sequoia, and (iv) Wukong ([Fe/H]$=-1.6$), a prograde phase-space overdensity chemically distinct from GSE. For each structure we provide [Fe/H], [$\alpha$/Fe], and orbital parameters. Stars born within the Galaxy are a major component at $|Z|\sim$2 kpc ($\approx$60$\%$), but their relative fraction declines sharply to $\lesssim$5$\%$ past 15 kpc. Beyond 15 kpc, $>$80$\%$ of the halo is built by two massive ($M_{\star}\sim10^{8}-10^{9}M_{\odot}$) accreted dwarfs: GSE ([Fe/H]$=-1.2$) within 25 kpc, and Sgr ([Fe/H]$=-1.0$) beyond 25 kpc. This explains the relatively high overall metallicity of the halo ([Fe/H]$\approx-1.2$). We attribute $\gtrsim$95$\%$ of the sample to one of the listed structures, pointing to a halo built entirely from accreted dwarfs and heating of the disk.

The primary method for inferring the stellar mass ($M_*$) of a galaxy is through spectral energy distribution (SED) modeling. However, the technique rests on assumptions such as the galaxy star formation history and dust attenuation law that can severely impact the accuracy of derived physical properties from SED modeling. Here, we examine the effect that the assumed star formation history (SFH) has on the stellar properties inferred from SED fitting by ground truthing them against mock observations of high-resolution cosmological hydrodynamic galaxy formation simulations. Classically, SFHs are modeled with simplified parameterized functional forms, but these forms are unlikely to capture the true diversity of galaxy SFHs and may impose systematic biases with under-reported uncertainties on results. We demonstrate that flexible nonparametric star formation histories outperform traditional parametric forms in capturing variations in galaxy star formation histories, and as a result, lead to significantly improved stellar masses in SED fitting. We find a decrease in the average bias of 0.4 dex with a delayed-$\tau$ model to a bias of just under 0.05 dex for the nonparametric model. Similarly, using nonparametric star formation histories in SED fitting result in increased accuracy in recovered galaxy star formation rates (SFRs) and stellar ages.

Previous studies have shown that the large, diffuse galaxies NGC1052-DF2 and NGC1052-DF4 both have populations of unusually luminous globular clusters as well as a very low dark matter content. Here we present newly-obtained deep Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) imaging of one of these galaxies, NGC1052-DF4. We use these data to measure the distance of the galaxy from the location of the tip of the red giant branch (TRGB). We find a rapid increase in the number of detected stars fainter than $m_{F814W} \sim 27.3$, which we identify as the onset of the red giant branch. Using a forward modeling approach that takes the photometric uncertainties into account, we find a TRGB magnitude of $m_{F814W,\rm TRGB}=27.47 \pm 0.16$. The inferred distance, including the uncertainty in the absolute calibration, is $D_{\rm TRGB}=20.0 \pm 1.6$ Mpc. The TRGB distance of NGC1052-DF4 is consistent with the previously-determined surface brightness fluctuation distance of $D_{\rm SBF}=18.7\pm 1.7$ Mpc to NGC1052-DF2 and is consistent with the distance of the bright elliptical galaxy NGC1052. We conclude that the unusual properties of these galaxies cannot be explained by distance errors.

The archaeological record of stars in the Milky Way opens a uniquely detailed window into the early formation and assembly of galaxies. Here we use 11,000 main-sequence turn-off stars with well-measured ages, [Fe/H], [$\alpha$/Fe], and orbits from the H3 Survey and Gaia to time the major events in the early Galaxy. Located beyond the Galactic plane, $1\lesssim |Z|/\rm kpc \lesssim4$, this sample contains three chemically distinct groups: a low metallicity population, and low-$\alpha$ and high-$\alpha$ groups at higher metallicity. The age and orbit distributions of these populations show that: 1) the high-$\alpha$ group, which includes both disk stars and the in-situ halo, has a star-formation history independent of eccentricity that abruptly truncated $8.3\pm0.1$ Gyr ago ($z\simeq1$); 2) the low metallicity population, which we identify as the accreted stellar halo, is on eccentric orbits and its star formation truncated $10.2.^{+0.2}_{-0.1}$ Gyr ago ($z\simeq2$); 3) the low-$\alpha$ population is primarily on low eccentricity orbits and the bulk of its stars formed less than 8 Gyr ago. These results suggest a scenario in which the Milky Way accreted a satellite galaxy at $z\approx2$ that merged with the early disk by $z\approx1$. This merger truncated star formation in the early high-$\alpha$ disk and perturbed a fraction of that disk onto halo-like orbits. The merger enabled the formation of a chemically distinct, low-$\alpha$ disk at $z\lesssim1$. The lack of any stars on halo-like orbits at younger ages indicates that this event was the last significant disturbance to the Milky Way disk.

We present spatially-resolved measurements of star formation histories (SFHs), metallicities, and distances in three nearby elliptical galaxies and the bulge of M31 derived using the pixel color-magnitude diagram (pCMD) technique. We compute pCMDs from archival $\textit{HST}$ photometry of M87, M49, NGC 3377 and M31, and fit the data using the new code $\texttt{PCMDPy}$. We measure distances to each system that are accurate to $\sim 10\%$. The recovered non-parametric SFHs place reasonable ($\pm 1$ dex) constraints on the recent (< 2 Gyr) star formation in M31 and NGC 3377, both of which show evidence of inside-out growth. The SFHs in M87 and M49 are constrained only at the oldest ages. The pCMD technique is a promising new avenue for studying the evolutionary history of the nearby universe, and is highly complementary to existing stellar population modeling techniques.

We present a description of the Dragonfly Wide Field Survey (DWFS), a deep photometric survey of a wide area of sky. The DWFS covers 330 $\mathrm{deg}^2$ in the equatorial GAMA fields and the Stripe 82 fields in the SDSS $g$ and $r$ bands. It is carried out with the 48-lens Dragonfly Telephoto Array, a telescope that is optimized for the detection of low surface brightness emission. The main goal of the survey is to study the dwarf galaxy population beyond the Local Group. In this paper, we describe the survey design and show early results. We reach $1\sigma$ depths of $\mu_g\approx 31$ mag arcsec$^{-2}$ on arcminute scales and show that Milky Way satellites such as Sextans, Bootes, and Ursa Major should be detectable out to $D\gtrsim 10$ Mpc. We also provide an overview of the elements and operation of the 48-lens Dragonfly telescope and a detailed description of its data reduction pipeline. The pipeline is fully automated, with individual frames subjected to a rigorous series of quality tests. The sky subtraction is performed in two stages, ensuring that emission features with spatial scales up to $\sim 0.^{\circ}9 \times 0.^{\circ}6$ are preserved. The DWFS provides unparalleled sensitivity to low surface brightness features on arcminute scales.

We present SPECULATOR - a fast, accurate, and flexible framework for emulating stellar population synthesis (SPS) models for predicting galaxy spectra and photometry. For emulating spectra, we use principal component analysis to construct a set of basis functions, and neural networks to learn the basis coefficients as a function of the SPS model parameters. For photometry, we parameterize the magnitudes (for the filters of interest) as a function of SPS parameters by a neural network. The resulting emulators are able to predict spectra and photometry under both simple and complicated SPS model parameterizations to percent-level accuracy, giving a factor of $10^3$-$10^4$ speed up over direct SPS computation. They have readily-computable derivatives, making them amenable to gradient-based inference and optimization methods. The emulators are also straightforward to call from a GPU, giving an additional order-of-magnitude speed-up. Rapid SPS computations delivered by emulation offers a massive reduction in the computational resources required to infer the physical properties of galaxies from observed spectra or photometry and simulate galaxy populations under SPS models, whilst maintaining the accuracy required for a range of applications.

We present an in-depth study of surface brightness fluctuations (SBFs) in low-luminosity stellar systems. Using the MIST models, we compute theoretical predictions for absolute SBF magnitudes in the LSST, HST, and proposed WFIRST filter systems. We compare our calculations to the observed SBF magnitudes of dwarf galaxies that have independent distance measurements from the tip of the red giant branch method. Consistent with previous studies, we find that single-age population models show excellent agreement with the observed SBF-color relation of low-mass galaxies with $0.5 \lesssim g - i \lesssim 0.9$. For bluer galaxies, the observed relation is better fit by models with composite stellar populations. To study SBF recovery from low-luminosity galaxies, we perform detailed image simulations in which we inject fully populated model galaxies into deep ground-based images from real observations. We demonstrate that measurements of SBF magnitudes from these simulated data correspond to the theoretical values with negligible bias ($\lesssim0.01$ mag). We then use the simulations to show that LSST will provide data of sufficient quality and depth to measure SBF distances with precisions of ${\sim}10$-20% to ultra-faint $\left(\mathrm{10^4 \leq M_\star/M_\odot \leq 10^5}\right)$ and low-mass classical ($\leq10^7$ M$_\odot$) dwarf galaxies out to ${\sim}4$ Mpc and ${\sim}25$ Mpc, respectively, within the first few years of its deep-wide-fast survey. Many systematic uncertainties remain, including an irreducible "sampling scatter" in the SBFs of ultra-faint dwarfs due to their undersampled stellar mass functions. We nonetheless conclude that SBFs in the new generation of wide-field imaging surveys have the potential to play a critical role in the efficient confirmation and characterization of dwarf galaxies in the nearby universe.

The protagonists of cosmic reionization remain elusive. Faint star-forming galaxies are leading candidates because they are numerous and may have significant ionizing photon escape fractions ($f_{esc}$). Here we update this picture via an empirical model that successfully predicts latest observations (e.g., the drop in star-formation density at z>8). We generate an ionizing spectrum for each galaxy in our model and constrain $f_{esc}$ using latest measurements of the reionization timeline (e.g., Ly$\alpha$ damping of quasars and galaxies at z>7). Assuming a constant $f_{esc}$, we find $M_{UV}$<-13.5 galaxies need $f_{esc}=0.21^{+0.06}_{-0.04}$ to complete reionization. The inferred IGM neutral fraction is [0.9, 0.5, 0.1] at z=[8.2, 6.8, 6.2]$\pm$0.2, i.e., the bulk of reionization transpires in 300 Myrs. Inspired by the emergent sample of Lyman Continuum (LyC) leakers that overwhelmingly displays higher-than-average star-formation surface density ($\Sigma$), we propose a model relating $f_{esc}$ to $\Sigma$ and find $f_{esc}\propto\Sigma^{0.4\pm0.1}$. Since $\Sigma$ falls by ~2.5 dex between z=8 and z=0, our model explains the humble upper limits on $f_{esc}$ at lower redshifts and its required evolution to ~0.2 at z>6. Within this model, strikingly, <5% of galaxies with $M_{UV}$<-18 (the `oligarchs') account for >80% of the reionization budget. In fact, faint sources ($M_{UV}$>-16) must be relegated to a limited role to ensure high neutral fractions at z=7-8. Shallow faint-end slopes of the UV luminosity function ($\alpha$>-2) and/or $f_{esc}$ distributions skewed toward bright galaxies produce the required late and rapid reionization. We predict LyC leakers like COLA1 (z=6.6, $f_{esc}$~30%, $M_{UV}$=-21.5) become increasingly common towards z~6 and that the drivers of reionization do not lie hidden across the faint-end of the luminosity function, but are already known to us. (abridged)

There has been a long-standing factor-of-two tension between the observed star formation rate density and the observed stellar mass buildup after $z\sim2$. Recently we have proposed that sophisticated panchromatic SED models can resolve this tension, as these methods infer systematically higher masses and lower star formation rates than standard approaches. In a series of papers we now extend this analysis and present a complete, self-consistent census of galaxy formation over $0.2 < z < 3$ inferred with the \texttt{Prospector} galaxy SED-fitting code. In this work, Paper I, we present the evolution of the galaxy stellar mass function using new mass measurements of $\sim$10$^5$ galaxies in the 3D-HST and COSMOS-2015 surveys. We employ a new methodology to infer the mass function from the observed stellar masses: instead of fitting independent mass functions in a series of fixed redshift intervals, we construct a continuity model that directly fits for the redshift evolution of the mass function. This approach ensures a smooth picture of galaxy assembly and makes use of the full, non-Gaussian uncertainty contours in our stellar mass inferences. The resulting mass function has higher number densities at a fixed stellar mass than almost any other measurement in the literature, largely owing to the older stellar ages inferred by \texttt{Prospector}. The stellar mass density is $\sim$50% higher than previous measurements, with the offset peaking at $z\sim1$. The next two papers in this series will present the new measurements of star-forming main sequence and the cosmic star formation rate density, respectively.

The $100^\circ$-long thin stellar stream in the Milky Way halo, GD-1, has an ensemble of features that may be due to dynamical interactions. Using high-resolution MMT/Hectochelle spectroscopy we show that a spur of GD-1-like stars outside of the main stream are kinematically and chemically consistent with the main stream. In the spur, as in the main stream, GD-1 has a low intrinsic radial velocity dispersion, $\sigma_{V_r}\lesssim1\,\rm km\,s^{-1}$, is metal-poor, $\rm [Fe/H]\approx-2.3$, with little $\rm [Fe/H]$ spread and some variation in $\rm [\alpha/Fe]$ abundances, which point to a common globular cluster progenitor. At a fixed location along the stream, the median radial velocity offset between the spur and the main stream is smaller than $0.5\,\rm km\,s^{-1}$, comparable to the measurement uncertainty. A flyby of a massive, compact object can change orbits of stars in a stellar stream and produce features like the spur observed in GD-1. In this scenario, the radial velocity of the GD-1 spur relative to the stream constrains the orbit of the perturber and its current on-sky position to $\approx5,000\,\rm deg^2$. The family of acceptable perturber orbits overlaps the stellar and dark-matter debris of the Sagittarius dwarf galaxy in present-day position and velocity. This suggests that GD-1 may have been perturbed by a globular cluster or an extremely compact dark-matter subhalo formerly associated with Sagittarius.

The extended main sequence turn offs (eMSTOs) of several young to intermediate age clusters are examined in the Magellanic Clouds and the Milky Way. We explore the effects of extended star formation (eSF) and a range of stellar rotation rates on the behavior of the color-magnitude diagram (CMD), paying particular attention to the MSTO. We create synthetic stellar populations based on MESA stellar models to simulate observed Hubble Space Telescope and Gaia star cluster data. We model the effect of rotation as a non-parametric distribution, allowing for maximum flexibility. In our models the slow rotators comprise the blueward, and fast rotators the redward portion of the eMSTO. We simulate data under three scenarios: non-rotating eSF, a range of rotation rates with a single age, and a combination of age and rotation effects. We find that two of the five clusters (the youngest and oldest) favor an age spread, but these also achieve the overall worst fits of all clusters. The other three clusters show comparable statistical evidence between rotation and an age spread. In all five cases, a rotation rate distribution alone is capable of qualitatively matching the observed eMSTO structure. In future work, we aim to compare our predicted Vsin(i) with observations in order to better constrain the physics related to stellar rotation.

We investigate the morphology of the colour-magnitude diagram (CMD) of the open cluster NGC 2509 in comparison with other Galactic open clusters of similar age using Gaia photometry. At $\sim900$ Myr Galactic open clusters in our sample all show an extended main sequence turn off (eMSTO) with the exception of NGC 2509, which presents an exceptionally narrow CMD. Our analysis of the Gaia data rules out differential extinction, stellar density, and binaries as a cause for the singular MSTO morphology in this cluster. We interpret this feature as a consequence of the stellar rotation distribution within the cluster and present the analysis with MIST stellar evolution models that include the effect of stellar rotation on which we based our conclusion. In particular, these models point to an unusually narrow range of stellar rotation rates ($\Omega/\Omega_{\rm{crit,ZAMS}} = [0.4, 0.6]$) within the cluster as the cause of this singular feature in the CMD of NGC 2509. Interestingly, models that do not include rotation are not as good at reproducing the morphology of the observed CMD in this cluster.

We present extremely deep imaging of the well-known edge-on galaxy NGC 4565 (the "Needle Galaxy"), observed as part of the Dragonfly Edge-on Galaxies Survey. Our imaging traces starlight over scales of ~100 kpc to surface brightnesses of ~32 mag/arcsec^2 in azimuthally averaged radial profiles. In narrow slice profiles along the major axis (with bin sizes ranging from 1.7 x 0.5 kpc to 1.7 x 7.8 kpc), we achieve a depth of ~29-30 mag/arcsec^2. We confirm the previously observed asymmetric disc truncation. In addition, the sharp northwest truncation turns over to a shallower component that coincides with a fan-like feature seen to wrap around the northwest disc limb. We propose that the fan may be a tidal ribbon, and qualitatively replicate the fan with simple simulations. Alternative explanations of the fan and the disc's asymmetry are considered as well. We conclude that we are most likely seeing evidence for accretion-based outer disk growth in NGC 4565.

The Galactic stellar halo is predicted to have formed at least partially from the tidal disruption of accreted dwarf galaxies. This assembly history should be detectable in the orbital and chemical properties of stars. The H3 Survey is obtaining spectra for 200,000 stars, and, when combined with Gaia data, is providing detailed orbital and chemical properties of Galactic halo stars. Unlike previous surveys of the halo, the H3 target selection is based solely on magnitude and Gaia parallax; the survey therefore provides a nearly unbiased view of the entire stellar halo at high latitudes. In this paper we present the distribution of stellar metallicities as a function of Galactocentric distance and orbital properties for a sample of 4232 kinematically-selected halo giants to 100 kpc. The stellar halo is relatively metal-rich, [Fe/H]=-1.2, and there is no discernable metallicity gradient over the range $6<R_{\rm gal}<100$ kpc. However, the halo metallicity distribution is highly structured including distinct metal-rich and metal-poor components at $R_{\rm gal}<10$ kpc and $R_{\rm gal}>30$ kpc, respectively. Metal-poor stars with [Fe/H]$<-2$ are a small population of the halo at all distances and orbital categories. We associate the "in-situ" stellar halo with stars displaying thick-disk chemistry on halo-like orbits; such stars are confined to $|z|<10$ kpc. The majority of the stellar halo is resolved into discrete features in orbital-chemical space, suggesting that the bulk of the stellar halo formed from the accretion and tidal disruption of dwarf galaxies. (ABRIDGED)

The timing argument provides a lower limit on the mass of the Milky Way. We find, using a sample of 32 stars at $R > 60$ kpc drawn from the H3 Spectroscopic Survey and mock catalogs created from published numerical simulations, that M$_{200} > 0.91\times 10^{12}$ M$_\odot$ with 90% confidence. We recommend using this limit to refine the allowed prior mass range in more complex and sophisticated statistical treatments of Milky Way dynamics. The use of such a prior would have significantly reduced many previously published uncertainty ranges. Our analysis suggests that the most likely value of M$_{200}$ is $\sim 1.4 \times 10^{12}$ M$_\odot$, but establishing this as the Milky Way mass requires a larger sample of outer halo stars and a more complete analysis of the inner halo stars in H3. The imminent growth in the sample of outer halo stars due to ongoing and planned surveys will make this possible.

In 2008 it was reported that the stellar stream of the edge-on spiral NGC5907 loops twice around the galaxy, enveloping it in a giant corkscrew-like structure. Here we present imaging of this iconic object with the Dragonfly Telephoto Array, reaching a $1\sigma$ surface brightness level of $\mu_g\approx 30.5$ mag/arcsec$^2$ on spatial scales of 1' (the approximate width of the stream). We find a qualitatively different morphology from that reported in the 2008 study. The Dragonfly data do not show two loops but a single curved stream with a total length of 45' (220 kpc). The surface brightness of the stream ranges from $\mu_g \approx 27.6$ mag/arcsec$^2$ to $\mu_g\approx 28.8$ mag/arcsec$^2$, and it extends significantly beyond the region where tidal features had previously been detected. We find a density enhancement near the luminosity-weighted midpoint of the stream which we identify as the likely remnant of a nearly-disrupted progenitor galaxy. A restricted N-body simulation provides a qualitative match to the detected features. In terms of its spatial extent and stellar mass the stream is similar to Sagittarius, and our results demonstrate the efficacy of low surface brightness-optimized telescopes for obtaining maps of such large streams outside the Local Group. The census of these rare, relatively high mass events complements the census of common, low mass ones that is provided by studies of streams in the Milky Way halo.

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) is a concept for a NASA probe-class space mission. It is the spectroscopic follow-up mission to WFIRST, boosting its scientific return by obtaining deep NIR & MIR slit spectroscopy for most of the galaxies imaged by the WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for ~200M galaxies out to z=7 and beyond, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS and WFIRST together will produce a definitive 3D map of the Universe over 2000 sq deg. ATLAS Science Goals are: (1) Discover how galaxies have evolved in the cosmic web of dark matter from cosmic dawn through the peak era of galaxy assembly. (2) Discover the nature of cosmic acceleration. (3) Probe the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy. (4) Discover the bulk compositional building blocks of planetesimals formed in the outer Solar System. These flow down to the ATLAS Scientific Objectives: (1A) Trace the relation between galaxies and dark matter with less than 10% shot noise on relevant scales at 1<z<7. (1B) Probe the physics of galaxy evolution at 1<z<7. (2) Obtain definitive measurements of dark energy and tests of General Relativity. (3) Measure the 3D structure and stellar content of the inner Milky Way to a distance of 25 kpc. (4) Detect and quantify the composition of 3,000 planetesimals in the outer Solar System. ATLAS is a 1.5m telescope with a FoV of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000, and a wavelength range of 1-4 microns. ATLAS has an unprecedented spectroscopic capability based on DMDs, with a spectroscopic multiplex factor ~6,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope.

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.

We present a method to flexibly and self-consistently determine individual galaxies' star formation rates (SFRs) from their host haloes' potential well depths, assembly histories, and redshifts. The method is constrained by galaxies' observed stellar mass functions, SFRs (specific and cosmic), quenched fractions, UV luminosity functions, UV-SM relations, IRX-UV relations, auto- and cross-correlation functions (including quenched and star-forming subsamples), and quenching dependence on environment; each observable is reproduced over the full redshift range available, up to 0<z<10. Key findings include: galaxy assembly correlates strongly with halo assembly; quenching at z>1 correlates strongly with halo mass; quenched fractions at fixed halo mass decrease with increasing redshift; massive quenched galaxies reside in higher-mass haloes than star-forming galaxies at fixed galaxy mass; star-forming and quenched galaxies' star formation histories at fixed mass differ most at z<0.5; satellites have large scatter in quenching timescales after infall, and have modestly higher quenched fractions than central galaxies; Planck cosmologies result in up to 0.3 dex lower stellar mass-halo mass ratios at early times; and, nonetheless, stellar mass-halo mass ratios rise at z>5. Also presented are revised stellar mass-halo mass relations for all, quenched, star-forming, central, and satellite galaxies; the dependence of star formation histories on halo mass, stellar mass, and galaxy SSFR; quenched fractions and quenching timescale distributions for satellites; and predictions for higher-redshift galaxy correlation functions and weak lensing surface densities. The public data release (DR1) includes the massively parallel (>10^5 cores) implementation (the UniverseMachine), the newly compiled and remeasured observational data, derived galaxy formation constraints, and mock catalogues including lightcones.

Modern theories of galaxy formation predict that the Galactic stellar halo was hierarchically assembled from the accretion and disruption of smaller systems. This hierarchical assembly is expected to produce a high degree of structure in the combined phase and chemistry space; this structure should provide a relatively direct probe of the accretion history of our Galaxy. Revealing this structure requires precise 3D positions (including distances), 3D velocities, and chemistry for large samples of stars. The Gaia satellite is delivering proper motions and parallaxes for >1 billion stars to G~20. However, radial velocities and metallicities will only be available to G~15, which is insufficient to probe the outer stellar halo (>10 kpc). Moreover, parallaxes will not be precise enough to deliver high-quality distances for stars beyond ~10 kpc. Identifying accreted systems throughout the stellar halo therefore requires a large ground-based spectroscopic survey to complement Gaia. Here we provide an overview of the H3 Stellar Spectroscopic Survey, which will deliver precise stellar parameters and spectrophotometric distances for 200,000 stars to r=18. Spectra are obtained with the Hectochelle instrument at the MMT, which is configured for the H3 Survey to deliver resolution R~23,000 spectra covering the wavelength range 5150A-5300A. The survey is optimized for stellar halo science and therefore focuses on high Galactic latitude fields (|b|>30 deg.), sparsely sampling 15,000 sq. degrees. Targets are selected on the basis of Gaia parallaxes, enabling very efficient selection of bone fide halo stars. The survey began in the Fall of 2017 and has collected 88,000 spectra to-date. All of the data, including the derived stellar parameters, will eventually be made publicly available via the survey website: h3survey.rc.fas.harvard.edu.

The UVJ color-color diagram is a popular and efficient method to distinguish between quiescent and star-forming galaxies through their rest-frame $U-V$ vs. $V-J$ colors. Here we explore the information content of this color-color space using the Bayesian inference machine Prospector. We fit the same physical model to two datasets: (i) UVJ fluxes alone, and (ii) full UV-mid IR (MIR) broadband SEDs from the 3D-HST survey. Notably this model uses both nonparametric SFHs and a flexible dust attenuation curve, both of which have the potential to `break' the typical correlations observed in UVJ color-color space. Instead, these fits confirm observed trends between UVJ colors and observed galaxy properties, including specific star formation rate (sSFR), dust attenuation, stellar age, and stellar metallicity. They also demonstrate that UVJ colors do not, on their own, constrain stellar age or metallicity; the observed trends in the UVJ diagram are instead driven by galaxy scaling relationships and thus will evolve with cosmological time. We also show that UVJ colors 'saturate' below $\log(\mathrm{sSFR/yr}^{-1})\lesssim -10.5$, i.e. changing sSFR no longer produces substantial changes in UVJ colors. We show that far-UV and/or MIR fluxes continue to correlate with sSFR down to low sSFRs and can be used in color-color diagrams to efficiently target galaxies with much lower levels of ongoing star formation. We provide selection criteria in these new color-color spaces as a function of desired sample sSFR.

We recommend a conceptual design study for a spectroscopic facility in the southern hemisphere comprising a large diameter telescope, fiber system, and spectrographs collectively optimized for massively-multiplexed spectroscopy. As a baseline, we propose an 11.4-meter aperture, optical spectroscopic survey telescope with a five square degree field of view. Using current technologies, the facility could be equipped with 15,000 robotically-controlled fibers feeding spectrographs over 360<lambda<1330 nm with options for fiber-fed spectrographs at high resolution and a panoramic IFU at a separate focus. This would enable transformational progress via its ability to access a larger fraction of objects from Gaia, LSST, Euclid, and WFIRST than any currently funded or planned spectroscopic facility. An ESO-sponsored study (arXiv:1701.01976) discussed the scientific potential in ambitious new spectroscopic surveys in Galactic astronomy, extragalactic astronomy, and cosmology. The US community should establish links with European and other international communities to plan for such a powerful facility and maximize the potential of large aperture multi-object spectroscopy given the considerable investment in deep imaging surveys.

The chemical composition of galaxies has been measured out to z~4. However, nearly all studies beyond z~0.7 are based on strong-line emission from HII regions within star-forming galaxies. Measuring the chemical composition of distant quiescent galaxies is extremely challenging, as the required stellar absorption features are faint and shifted to near-infrared wavelengths. Here, we present ultra-deep rest-frame optical spectra of five massive quiescent galaxies at z~1.4, all of which show numerous stellar absorption lines. We derive the abundance ratios [Mg/Fe] and [Fe/H] for three out of five galaxies; the remaining two galaxies have too young luminosity-weighted ages to yield robust measurements. Similar to lower-redshift findings, [Mg/Fe] appears positively correlated with stellar mass, while [Fe/H] is approximately constant with mass. These results may imply that the stellar mass-metallicity relation was already in place at z~1.4. While the [Mg/Fe]-mass relation at z~1.4 is consistent with the z<0.7 relation, [Fe/H] at z~1.4 is ~0.2 dex lower than at z<0.7. With a [Mg/Fe] of 0.44(+0.08,-0.07) the most massive galaxy may be more alpha-enhanced than similar-mass galaxies at lower redshift, but the offset is less significant than the [Mg/Fe] of 0.6 previously found for a massive galaxy at z=2.1. Nonetheless, these results combined may suggest that [Mg/Fe] in the most massive galaxies decreases over time, possibly by accreting low-mass, less alpha-enhanced galaxies. A larger galaxy sample is needed to confirm this scenario. Finally, the abundance ratios indicate short star-formation timescales of 0.2-1.0 Gyr.

We present a model for the interaction of the GD-1 stellar stream with a massive perturber that naturally explains many of the observed stream features, including a gap and an off-stream spur of stars. The model involves an impulse by a fast encounter, after which the stream grows a loop of stars at different orbital energies. At specific viewing angles, this loop appears offset from the stream track. A quantitative comparison of the spur and gap features prefers models where the perturber is in the mass range of $10^6\,\rm M_\odot$ to $10^8\,\rm M_\odot$. Orbit integrations back in time show that the stream encounter could not have been caused by any known globular cluster or dwarf galaxy with a determined orbit, and mass, size and impact-parameter arguments show that it could not have been caused by a molecular cloud in the Milky Way disk. The most plausible explanation for the gap-and-spur structure is an encounter with a dark-matter substructure, like those predicted to populate galactic halos in LCDM cosmology. However, the expected densities of LCDM subhalos in this mass range and in this part of the Milky Way are $2-3\,\sigma$ lower than the inferred density of the GD-1 perturber. This observation opens up the possibility that detailed observations of streams could measure the mass spectrum of dark-matter substructures and even identify individual substructures and their orbits in the Galactic halo.

Galaxy observations are influenced by many physical parameters: stellar masses, star formation rates (SFRs), star formation histories (SFHs), metallicities, dust, black hole activity, and more. As a result, inferring accurate physical parameters requires high-dimensional models which capture or marginalize over this complexity. Here we re-assess inferences of galaxy stellar masses and SFRs using the 14-parameter physical model Prospector-$\alpha$ built in the Prospector Bayesian inference framework. We fit the photometry of 58,461 galaxies from the 3D-HST catalogs at $0.5 < z < 2.5$. The resulting stellar masses are $\sim0.1-0.3$ dex larger than the fiducial masses while remaining consistent with dynamical constraints. This change is primarily due to the systematically older SFHs inferred with Prospector. The SFRs are $\sim0.1-1+$ dex lower than UV+IR SFRs, with the largest offsets caused by emission from "old" ($t>100$ Myr) stars. These new inferences lower the observed cosmic star formation rate density by $\sim0.2$ dex and increase the observed stellar mass growth by $\sim 0.1$ dex, finally bringing these two quantities into agreement and implying an older, more quiescent Universe than found by previous studies at these redshifts. We corroborate these results by showing that the Prospector-$\alpha$ SFHs are both more physically realistic and are much better predictors of the evolution of the stellar mass function. Finally, we highlight examples of observational data which can break degeneracies in the current model; these observations can be incorporated into priors in future models to produce new & more accurate physical parameters.

In simple models of the Milky Way, tidally disrupting satellites produce long and thin---nearly one-dimensional---stellar streams. Using astrometric data from the Gaia second data release and photometry from the Dark Energy Survey, we demonstrate that the Jhelum stream, a stellar stream in the inner halo, is a two-dimensional structure. The spatial distribution of highly probable Jhelum members reveals a dense thin component and an associated diffuse, spatially offset component. These two spatial components have indistinguishable proper motions (at $\sigma\sim1\,\rm mas\,yr^{-1}$ level) and a similar ratio of blue straggler to blue horizontal branch stars, which indicates a common origin for the two components. The best-fit orbit to the narrow component (pericenter $8\,\rm kpc$, apocenter $24\,\rm kpc$), however, does not explain the wide component of the Jhelum stream. On the other hand, an older orbital wrap of Jhelum's orbit traces the Indus stream, indicating a possible connection between these two structures and additional complexity in Jhelum's formation. Substructure in the Jhelum progenitor or precession of its tidal debris in the Milky Way potential may explain the observed structure of Jhelum. Future spectroscopic data will enable discrimination between these "nature" and "nurture" formation scenarios. Jhelum adds to the growing list of cold stellar streams that display complex morphologies and promise to reveal the dynamical history of the Milky Way.

We present The Payne, a general method for the precise and simultaneous determination of numerous stellar labels from observed spectra, based on fitting physical spectral models. The Payne combines a number of important methodological aspects: it exploits the information from much of the available spectral range; it fits all labels (stellar parameters and element abundances) simultaneously; it uses spectral models, where the atmosphere structure and the radiative transport are consistently calculated to reflect the stellar labels. At its core The Payne has an approach to accurate and precise interpolation and prediction of the spectrum in high-dimensional label-space, which is flexible and robust, yet based on only a moderate number of ab initio models (O(1000) for 25 labels). With a simple neural-net-like functional form and a suitable choice of training labels, this interpolation yields a spectral flux prediction good to $10^{-3}$ rms across a wide range of $T_{\rm eff}$ and log g (including dwarfs and giants). We illustrate the power of this approach by applying it to the APOGEE DR14 data set, drawing on Kurucz models with recently improved line lists: without recalibration, we obtain physically sensible stellar parameters as well as 15 element abundances that appear to be more precise than the published APOGEE DR14 values. In short, The Payne is an approach that for the first time combines all these key ingredients, necessary for progress towards optimal modelling of survey spectra; and it leads to both precise and accurate estimates of stellar labels, based on physical models and without re-calibration. Both the codes and catalog are made publicly available online.

Early type galaxies (ETGs) frequently show emission from warm ionized gas. These Low Ionization Emission Regions (LIERs) were originally attributed to a central, low-luminosity active galactic nuclei. However, the recent discovery of spatially-extended LIER emission suggests ionization by both a central source and an extended component that follows a stellar-like radial distribution. For passively-evolving galaxies with old stellar populations, hot post-Asymptotic Giant Branch (AGB) stars are the only viable extended source of ionizing photons. In this work, we present the first prediction of LIER-like emission from post-AGB stars that is based on fully self-consistent stellar evolution and photoionization models. We show that models where post-AGB stars are the dominant source of ionizing photons reproduce the nebular emission signatures observed in ETGs, including LIER-like emission line ratios in standard optical diagnostic diagrams and H$\alpha$ equivalent widths of order 0.1-3 angstroms. We test the sensitivity of LIER-like emission to the details of post-AGB models, including the mass loss efficiency and convective mixing efficiency, and show that line strengths are relatively insensitive to post-AGB timescale variations. Finally, we examine the UV-optical colors of the models and the stellar populations responsible for the UV-excess observed in some ETGs. We find that allowing as little as 3% of the HB population to be uniformly distributed to very hot temperatures (30,000 K) produces realistic UV colors for old, quiescent ETGs.

Gaia is currently revolutionizing modern astronomy. However, much of the Galactic plane, center and the spiral arm regions are obscured by interstellar extinction, rendering them inaccessible because Gaia is an optical instrument. An all-sky near infrared (NIR) space observatory operating in the optical NIR, separated in time from the original Gaia would provide microarcsecond NIR astrometry and millimag photometry to penetrate obscured regions unraveling the internal dynamics of the Galaxy.

Nonparametric star formation histories (SFHs) have long promised to be the `gold standard' for galaxy spectral energy distribution (SED) modeling as they are flexible enough to describe the full diversity of SFH shapes, whereas parametric models rule out a significant fraction of these shapes {\it a priori}. However, this flexibility is not fully constrained even with high-quality observations, making it critical to choose a well-motivated prior. Here, we use the SED-fitting code \texttt{Prospector} to explore the effect of different nonparametric priors by fitting SFHs to mock UV-IR photometry generated from a diverse set of input SFHs. First, we confirm that nonparametric SFHs recover input SFHs with less bias and return more accurate errors than do parametric SFHs. We further find that, while nonparametric SFHs robustly recover the overall shape of the input SFH, the primary determinant of the size and shape of the posterior star formation rate (SFR) as a function of time is the choice of prior, rather than the photometric noise. As a practical demonstration, we fit the UV-IR photometry of $\sim$6000 galaxies from the GAMA survey and measure inter-prior scatters in mass (0.1 dex), SFR$_{100\; \mathrm{Myr}}$ (0.8 dex), and mass-weighted ages (0.2 dex), with the bluest star-forming galaxies showing the most sensitivity. An important distinguishing characteristic for nonparametric models is the characteristic timescale for changes in SFR(t). This difference controls whether galaxies are assembled in bursts or in steady-state star formation, corresponding respectively to (feedback-dominated/accretion-dominated) models of galaxy formation and to (larger/smaller) confidence intervals derived from SED-fitting. High-quality spectroscopy has the potential to further distinguish between these proposed models of SFR(t).

It is challenging to reliably identify stars that were born together outside of actively star-forming regions and bound stellar systems. However, co-natal stars should be present throughout the Galaxy, and their demographics can shed light on the clustered nature of star formation and the dynamical state of the disk. In previous work we presented a set of simulations of the Galactic disk that followed the clustered formation and dynamical evolution of 4 billion individual stars over the last 5 Gyr. The simulations predict that a high fraction of co-moving stars with physical and 3D velocity separation of $\Delta r < 20$ pc and $\Delta v < 1.5$ km s$^{-1}$ are co-natal. In this \textit{Letter}, we use \textit{Gaia} DR2 and LAMOST DR4 data to identify and study co-moving pairs. We find that the distribution of relative velocities and separations of pairs in the data is in good agreement with the predictions from the simulation. We identify 111 co-moving pairs in the Solar neighborhood with reliable astrometric and spectroscopic measurements. These pairs show a strong preference for having similar metallicities when compared to random field pairs. We therefore conclude that these pairs were very likely born together. The simulations predict that co-natal pairs originate preferentially from high-mass and relatively young ($< 1$ Gyr) star clusters. \textit{Gaia} will eventually deliver well-determined metallicities for the brightest stars, enabling the identification of thousands of co-natal pairs due to disrupting star clusters in the solar neighborhood.

We present a comprehensive study of the applications of the pixel color-magnitude diagram (pCMD) technique for measuring star formation histories (SFHs) and other stellar population parameters of galaxies, and demonstrate that the technique can also constrain distances. SFHs have previously been measured through either the modeling of resolved-star CMDs or of integrated-light SEDs, yet neither approach can easily be applied to galaxies in the "semi-resolved regime". The pCMD technique has previously been shown to have the potential to measure stellar populations and star formation histories in semi-resolved galaxies. Here we present Pixel Color-Magnitude Diagrams with Python (pcmdpy), a GPU-accelerated package that makes significant computational improvements to the original code and including more realistic physical models. These advances include the simultaneous fitting of distance, modeling a Gaussian metallicity-distribution function, and an observationally-motivated dust model. GPU-acceleration allows these more realistic models to be fit roughly 7x faster than the simpler models in the original code. We present results from a suite of mock tests, showing that with proper model assumptions, the code can simultaneously recover SFH, [Fe/H], distance, and dust extinction. Our results suggest the code, applied to observations with HST-like resolution, should constrain these properties with high precision within 10 Mpc and can be applied to systems out to as far as 100 Mpc. pCMDs open a new window to studying the stellar populations of many galaxies that cannot be readily studied through other means.

The ultra-diffuse galaxy NGC1052-DF2 has a very low velocity dispersion, indicating that it has little or no dark matter. Here we report the discovery of a second galaxy in this class, residing in the same group. NGC1052-DF4 closely resembles NGC1052-DF2 in terms of its size, surface brightness, and morphology; has a similar distance of $D=19.9\pm 2.8$ Mpc; and also has a population of luminous globular clusters extending out to 7 kpc from the center of the galaxy. Accurate radial velocities of the diffuse galaxy light and seven of the globular clusters were obtained with the Low Resolution Imaging Spectrograph on the Keck I telescope. The velocity of the diffuse light is identical to the median velocity of the clusters, $v_{\rm sys}=\langle v_{\rm gc} \rangle=1445$ km/s, and close to the central velocity of the NGC1052 group. The rms spread of the observed velocities is very small at $\sigma_{\rm obs}=5.8$ km/s. Taking observational uncertainties into account we determine an intrinsic velocity dispersion of $\sigma_{\rm intr}=4.2^{+4.4}_{-2.2}$ km/s, consistent with the expected value from the stars alone ($\sigma_{\rm stars}\approx 7$ km/s) and lower than expected from a standard NFW halo ($\sigma_{\rm halo}\sim 30$ km/s). We conclude that NGC1052-DF2 is not an isolated case but that a class of such objects exists. The origin of these large, faint galaxies with an excess of luminous globular clusters and an apparent lack of dark matter is, at present, not understood.

Galaxies form and evolve in the context of their local and large-scale environments. Their baryonic content that we observe with imaging and spectroscopy is intimately connected to the properties of their dark matter halos, and to their location in the "cosmic web" of large-scale structure. Very large spectroscopic surveys of the local universe (e.g., SDSS and GAMA) measure galaxy positions (location within large-scale structure), statistical clustering (a direct constraint on dark matter halo masses), and spectral features (measuring physical conditions of the gas and stars within the galaxies, as well as internal velocities). Deep surveys with the James Webb Space Telescope (JWST) will revolutionize spectroscopic measurements of redshifts and spectral properties for galaxies out to the epoch of reionization, but with numerical statistics and over cosmic volumes that are too small to map large-scale structure and to constrain halo properties via clustering. Here, we consider advances in understanding galaxy evolution that would be enabled by very large spectroscopic surveys at high redshifts: very large numbers of galaxies (outstanding statistics) over large co-moving volumes (large-scale structure on all scales) over broad redshift ranges (evolution over most of cosmic history). The required observational facility can be established as part of the probe portfolio by NASA within the next decade.

Over the past decade, empirical constraints on the galaxy-dark matter halo connection have significantly advanced our understanding of galaxy evolution. Past techniques have focused on connections between halo properties and galaxy stellar mass and/or star formation rates. Empirical techniques in the next decade will link halo assembly histories with galaxies' circumgalactic media, supermassive black holes, morphologies, kinematics, sizes, colors, metallicities, and transient rates. Uncovering these links will resolve many critical uncertainties in galaxy formation and will enable much higher-fidelity mock catalogs essential for interpreting observations. Achieving these results will require broader and deeper spectroscopic coverage of galaxies and their circumgalactic media; survey teams will also need to meet several criteria (cross-comparisons, public access, and covariance matrices) to facilitate combining data across different surveys. Acting on these recommendations will continue enabling dramatic progress in both empirical modeling and galaxy evolution for the next decade.

The velocity dispersion of the ultra diffuse galaxy NGC1052-DF2 was found to be $\sigma_{\rm gc}=7.8^{+5.2}_{-2.2} \ \mathrm{kms^{-1}}$, much lower than expected from the stellar mass -- halo mass relation and nearly identical to the expected value from the stellar mass alone. This result was based on the radial velocities of ten luminous globular clusters that were assumed to be associated with the galaxy. A more precise measurement is possible from high resolution spectroscopy of the diffuse stellar light. Here we present an integrated spectrum of the diffuse light of NGC1052-DF2 obtained with the Keck Cosmic Web Imager, with an instrumental resolution of $\sigma_{\rm instr}\approx 12 \ \mathrm{kms^{-1}}$. The systemic velocity of the galaxy is $v_{\rm sys}=1805\pm 1.1 \ \mathrm{kms^{-1}}$, in very good agreement with the average velocity of the globular clusters ($\langle v_{\rm gc}\rangle = 1803\pm 2 \ \mathrm{kms^{-1}}$). There is no evidence for rotation within the KCWI field of view. We find a stellar velocity dispersion of $\sigma_{\rm stars}=8.4 \pm 2.1 \ \mathrm{kms^{-1}}$, consistent with the dispersion that was derived from the globular clusters. The implied dynamical mass within the half-light radius $r_{1/2}=2.7 \ \mathrm{kpc}$ is $M_{\rm dyn}= (1.3 \pm 0.8) \times 10^8 $ M$_{\odot}$, similar to the stellar mass within that radius ($M_{\rm stars}=(1.0 \pm 0.2) \times 10^8 \ \mathrm{M}_{\odot}$). With this confirmation of the low velocity dispersion of NGC1052-DF2, the most urgent question is whether this "missing dark matter problem" is unique to this galaxy or applies more widely.

Parametric models for galaxy star-formation histories (SFHs) are widely used, though they are known to impose strong priors on physical parameters. This has consequences for measurements of the galaxy stellar-mass function (GSMF), star-formation-rate density (SFRD) and star-forming main sequence (SFMS). We investigate the effects of the exponentially declining, delayed exponentially declining, lognormal and double power law SFH models using BAGPIPES. We demonstrate that each of these models imposes strong priors on specific star-formation rates (sSFRs), potentially biasing the SFMS, and also imposes a strong prior preference for young stellar populations. We show that stellar mass, SFR and mass-weighted age inferences from high-quality mock photometry vary with the choice of SFH model by at least 0.1, 0.3 and 0.2 dex respectively. However the biases with respect to the true values depend more on the true SFH shape than the choice of model. We also demonstrate that photometric data cannot discriminate between SFH models, meaning it is important to perform independent tests to find well-motivated priors. We finally fit a low-redshift, volume-complete sample of galaxies from the Galaxy and Mass Assembly (GAMA) Survey with each model. We demonstrate that our stellar masses and SFRs at redshift, $z\sim0.05$ are consistent with other analyses. However, our inferred cosmic SFRDs peak at $z\sim0.4$, approximately 6 Gyr later than direct observations suggest, meaning our mass-weighted ages are significantly underestimated. This makes the use of parametric SFH models for understanding mass assembly in galaxies challenging. In a companion paper we consider non-parametric SFH models.

The clustered nature of star formation should produce a high degree of structure in the combined phase and chemical space in the Galactic disk. To date, observed structure of this kind has been mostly limited to bound clusters and moving groups. In this paper we present a new dynamical model of the Galactic disk that takes into account the clustered nature of star formation. This model predicts that the combined phase and chemical space is rich in substructure, and that this structure is sensitive to both the precise nature of clustered star formation and the large-scale properties of the Galaxy. The model self-consistently evolves 4 billion stars over the last 5 Gyr in a realistic potential that includes an axisymmetric component, a bar, spiral arms, and giant molecular clouds (GMCs). All stars are born in clusters with an observationally-motivated range of initial conditions. As direct \textit{N}-body calculations for billions of stars is computationally infeasible, we have developed a method of initializing star cluster particles to mimic the effects of direct \textit{N}-body effects, while the actual orbit integrations are treated as test particles within the analytic potential. We demonstrate that the combination of chemical and phase space information is much more effective at identifying truly co-natal populations than either chemical or phase space alone. Furthermore, we show that co-moving pairs of stars are very likely to be co-natal if their velocity separation is $< 2$ km s$^{-1}$ and their metallicity separation is $< 0.05$ dex. The results presented here bode well for harnessing the synergies between \textit{Gaia} and spectroscopic surveys to reveal the assembly history of the Galactic disk.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

It has been suggested that the dark matter-deficient galaxies NGC1052-DF2 and NGC1052-DF4 might not be members of the NGC1052 group but in the foreground at $\sim 13$ Mpc, and satellites of the bright spiral galaxy NGC1042. We previously showed that the CMDs of the galaxies are inconsistent with this hypothesis, and derived distances of 19-20 Mpc from their surface brightness fluctuation signals. Here we note that NGC1042 is almost certainly a member of the NGC1052 group as well, based on its radial velocity, the HI distribution in the NGC1052/NGC1042 system, and the Tully-Fisher relation.

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) Probe is a concept for a NASA probe-class space mission. It is the follow-up space mission to WFIRST, boosting its scientific return by obtaining deep IR slit spectroscopy for 70% of all galaxies imaged by a 2000 sq deg WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for 200M galaxies out to z < 7, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS Probe science spans four broad categories: (1) Revolutionizing galaxy evolution studies by tracing the relation between galaxies and dark matter from galaxy groups to cosmic voids and filaments, from the epoch of reionization through the peak era of galaxy assembly; (2) Opening a new window into the dark Universe by weighing the dark matter filaments using 3D weak lensing with spectroscopic redshifts, and obtaining definitive measurements of dark energy and modification of General Relativity using galaxy clustering; (3) Probing the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy; and (4) Exploring the formation history of the outer Solar System by characterizing Kuiper Belt Objects. ATLAS Probe is a 1.5m telescope with a field of view of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000 over 1-4 microns, and a spectroscopic multiplex factor >5,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope; it has a single instrument, a telescope aperture that allows for a lighter launch vehicle, and mature technology. ATLAS Probe will lead to transformative science over the entire range of astrophysics: from galaxy evolution to the dark Universe, from Solar System objects to the dusty regions of the Milky Way.

Stellar population synthesis (SPS) models have long struggled to reproduce observed optical through near-IR (NIR) spectral energy distributions (SED) of massive quiescent galaxies. We revisit this issue using a novel approach that combines the diagnostic power of full-spectrum fitting with recently updated stellar spectral libraries. First, we perform full-spectrum fitting of continuum-normalized stacked SDSS spectra in bins of velocity dispersion to infer their stellar population properties, such as the elemental abundances and age. Next, we use the resulting best-fit parameters to compute $ugriz$ colors, which are then compared to observed colors of the same galaxies. With this approach we are able to predict the $ugriz$ SEDs of low and high mass galaxies at the $\lesssim 0.03$ mag level in nearly all cases. We find that the full optical through NIR SEDs of quiescent galaxies can be reproduced only when the spectrum is fit with a flexibility that is able to capture the behavior of the entire optical absorption line spectrum. The models include variations in individual elemental abundances, nebular emission lines, and the presence of young stellar components. The successful prediction of the SED shape from continuum-normalized spectra implies that the continuum information is largely contained in the narrow absorption features. These results also imply that attempts to model broadband photometry of quiescent systems will suffer from potentially significant biases if the detailed abundance patterns are not taken into account.

While we observe a large amount of cold interstellar gas and dust in a subset of the early-type galaxy (ETG) population, the source of this material remains unclear. The two main, competing scenarios are external accretion of lower mass, gas-rich dwarfs and internal production from stellar mass loss and/or cooling from the hot interstellar medium (ISM). We test these hypotheses with measurements of the stellar and nebular metallicities of three ETGs (NGC 2768, NGC 3245, and NGC 4694) from new long-slit, high signal-to-noise ratio spectroscopy from the Multi-Object Double Spectographs (MODs) on the Large Binocular Telescope (LBT). These ETGs have modest star formation rates and minimal evidence of nuclear activity. We model the stellar continuum to derive chemical abundances and measure gas-phase abundances with standard nebular diagnostics. We find that the stellar and gas-phase abundances are very similar, which supports internal production and is very inconsistent with the accretion of smaller, lower metallicity dwarfs. All three of these galaxies are also consistent with an extrapolation of the mass-metallicity relation to higher mass galaxies with lower specific star formation rates. The emission line flux ratios along the long-slit, as well as global line ratios clearly indicate that photoionization dominates and ionization by alternate sources including AGN activity, shocks, cosmic rays, dissipative magnetohydrodynamic waves, and single degenerate Type Ia supernovae progenitors do not significantly affect the line ratios.

We explore the connection between the UV luminosity functions (LFs) of high-$z$ galaxies and the distribution of stellar masses and star-formation histories (SFHs) in their host dark matter halos. We provide a baseline for a redshift-independent star-formation efficiency model to which observations and models can be compared. Our model assigns a star-formation rate (SFR) to each dark matter halo based on the growth rate of the halo and a redshift-independent star-formation efficiency. The dark matter halo accretion rate is obtained from a high-resolution $N$-body simulation in order to capture the stochasticity in accretion histories and to obtain spatial information for the distribution of galaxies. The halo mass dependence of the star-formation efficiency is calibrated at $z=4$ by requiring a match to the observed UV LF at this redshift. The model then correctly predicts the observed UV LF at $z=5-10$. We present predictions for the UV luminosity and stellar mass functions, \textit{JWST} number counts, and SFHs. In particular, we find a stellar-to-halo mass relation at $z=4-10$ that scales with halo mass at $M_{\rm h}<10^{11}~M_{\odot}$ as $M_{\star}\propto M_{\rm h}^2$, with a normalization that is higher than the relation inferred at $z=0$. The average SFRs increase as a function of time to $z=4$, although there is significant scatter around the average: about 6\% of the $z=4$ galaxies show no significant mass growth. Using these SFHs, we present redshift-dependent UV-to-SFR conversion factors, mass return fractions, and mass-to-light ratios for different intial mass functions and metallicities, finding that current estimates of the cosmic SFR density at $z\sim10$ may be overestimated by $\sim0.1-0.2~\mathrm{dex}$.

We present a comparison of the observed, spatially integrated stellar and ionized gas velocity dispersions of $\sim1000$ massive ($\log M_{\star}/M_{\odot}\gtrsim\,10.3$) galaxies in the Large Early Galaxy Astrophysics Census (LEGA-C) survey at $0.6\lesssim\,z\lesssim1.0$. The high $S/N\sim20{\rm\AA^{-1}}$ afforded by 20 hour VLT/VIMOS spectra allows for joint modeling of the stellar continuum and emission lines in all galaxies, spanning the full range of galaxy colors and morphologies. These observed integrated velocity dispersions (denoted as $\sigma'_{g, int}$ and $\sigma'_{\star, int}$) are related to the intrinsic velocity dispersions of ionized gas or stars, but also include rotational motions through beam smearing and spectral extraction. We find good average agreement between observed velocity dispersions, with $\langle\log(\sigma'_{g, int}/\sigma'_{\star, int})\rangle=-0.003$. This result does not depend strongly on stellar population, structural properties, or alignment with respect to the slit. However, in all regimes we find significant scatter between $\sigma'_{g, int}$ and $\sigma'_{\star, int}$, with an overall scatter of 0.13 dex of which 0.05 dex is due to observational uncertainties. For an individual galaxy, the scatter between $\sigma'_{g, int}$ and $\sigma'_{\star, int}$ translates to an additional uncertainty of $\sim0.24\rm{dex}$ on dynamical mass derived from $\sigma'_{g, int}$, on top of measurement errors and uncertainties from Virial constant or size estimates. We measure the $z\sim0.8$ stellar mass Faber-Jackson relation and demonstrate that emission line widths can be used to measure scaling relations. However, these relations will exhibit increased scatter and slopes that are artificially steepened by selecting on subsets of galaxies with progressively brighter emission lines.

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 the \prospector\ spectral energy distribution (SED) fitting code, we analyze the properties of 19 Extreme Emission Line Galaxies (EELGs) identified in the bluest composite SED in the \zfourge\ survey at $2.5 \leq z \leq 4$. \prospector\ includes a physical model for nebular emission and returns probability distributions for stellar mass, stellar metallicity, dust attenuation, and nonparametric star formation history (SFH). The EELGs show evidence for a starburst in the most recent 50 Myr, with the median EELG having a specific star formation rate (sSFR) of 4.6 Gyr$^{-1}$ and forming 15\% of its mass in this short time. For a sample of more typical star-forming galaxies (SFGs) at the same redshifts, the median SFG has a sSFR of 1.1 Gyr$^{-1}$ and forms only $4\%$ of its mass in the last 50 Myr. We find that virtually all of our EELGs have rising SFHs, while most of our SFGs do not. From our analysis, we hypothesize that many, if not most, star-forming galaxies at $z \geq 2.5$ undergo an extreme H$\beta$+$[\hbox{{\rm O}\kern 0.1em{\sc iii}}]$ emission line phase early in their lifetimes. In a companion paper, we obtain spectroscopic confirmation of the EELGs as part of our {\sc MOSEL} survey. In the future, explorations of uncertainties in modeling the UV slope for galaxies at $z>2$ are needed to better constrain their properties, e.g. stellar metallicities.

Galactic globular clusters (GCs) show overwhelming photometric and spectroscopic evidence for the existence of multiple stellar populations. The question of whether or not there exists a GC that represents a true 'simple stellar population' remains open. Here we focus on Ruprecht 106 (R106), a halo GC with [Fe/H]=-1.5 and [alpha/Fe]~0. A previous spectroscopic study found no sign of the Na-O anticorrelation among 9 of its brightest red giants, which led to the conclusion that R106 is a true simple stellar population GC. Here we present new Hubble Space Telescope (HST) Wide Field Camera 3 photometry of R106 that, when combined with archival HST images spanning a 6-year baseline, allows us to create proper motion cleaned color-magnitude diagrams spanning the ultraviolet (F336W) to the near-infrared (F814W). These data allow us to construct the pseudo-color C_{U,B,I} that is sensitive to the presence of light-element abundance spreads. We find no evidence of a split along the red giant branch (RGB) in the C_{U,B,I} diagram but the width of the RGB (sigma_CUBI = 0.015) is marginally broader than expected from artificial star tests (sigma_CUBI = 0.009). The observed spread in C_{U,B,I} is smaller than any other Galactic GC studied to date. Our results raise important questions about the role of formation environment and primordial chemical composition in the formation of multiple stellar populations in GCs.

We recently inferred that the galaxy NGC1052-DF2 has little or no dark matter and a rich system of unusual globular clusters. We assumed that the galaxy is a satellite of the luminous elliptical galaxy NGC1052 at ~20 Mpc, on the basis of its surface brightness fluctuations (SBF) distance of $19.0 \pm 1.7$ Mpc, its radial velocity of ~1800 km/s, and its projected position. Here we analyze the color-magnitude diagram (CMD) of NGC1052-DF2, following the suggestion by Trujillo et al. (2018) that the tip of the red giant branch (TRGB) can be detected in currently available HST data and the galaxy is at ~13 Mpc. Using fully populated galaxy models we show that the CMD is strongly influenced by blends. These blends produce a "phantom" TRGB ~2 times brighter than the true TRGB, which can lead to erroneous distance estimates ~1.4 times smaller than the actual distance. We compare NGC1052-DF2 to model images as well as other galaxies in our HST sample, and show that the large population of unblended RGB stars expected for distances of ~13 Mpc is not detected. We also provide a new distance measurement to NGC1052-DF2 that is free of calibration uncertainties, by anchoring it to a satellite of the megamaser host galaxy NGC4258. From a megamaser-TRGB-SBF distance ladder we obtain $D=18.7 \pm 1.7$ Mpc, consistent with our previous measurement and with the distance to the elliptical galaxy NGC1052.

There is a long history of using optical emission and absorption lines to constrain the metallicity and ionization parameters of gas in galaxies. However, comparable diagnostics are less well-developed for the UV. Here, we assess the diagnostic potential of both absorption and emission features in the UV and evaluate the diagnostics against observations of local and high redshift galaxies. We use the CloudyFSPS nebular emission model of Byler et al. 2017, extended to include emission predictions in the UV, to evaluate the metallicity sensitivity of established UV stellar absorption indices, and to identify those that include a significant contribution from nebular emission. We present model UV emission line fluxes as a function of metallicity and ionization parameter, assuming both instantaneous bursts and constant SFRs. We identify combinations of strong emission lines that constrain metallicity and ionization parameter, including [CIII] 1907, CIII] 1909, OIII] 1661,1666, SiIII]1883,1892, CIV 1548,1551, NII] 1750,1752, and MgII 2796, and develop UV versions of the canonical "BPT" diagram. We quantify the relative contribution from stellar wind emission and nebular line emission to diagnostic line ratios that include the CIV 1548,1551 lines, and also develop an observationally motivated relationship for N and C enrichment that improves the performance of photoionization models. We summarize the best diagnostic choices and the associated redshift range for low-, mid-, and high-resolution rest-UV spectroscopy in preparation for the launch of the James Webb Space Telescope.

The discovery of the ultra-diffuse galaxy NGC 1052-DF2 and its peculiar population of star clusters has raised new questions about the connections between galaxies and dark matter halos at the extremes of galaxy formation. In light of debates over the measured velocity dispersion of its star clusters and the associated mass estimate, we constrain mass models of DF2 using its observed kinematics with a range of priors on the halo mass. Models in which the galaxy obeys a standard stellar-halo mass relation are in tension with the data and also require a large central density core. Better fits are obtained when the halo mass is left free, even after accounting for increased model complexity. The dynamical mass-to-light ratio for our model with a weak prior on the halo mass is $1.7^{+0.7}_{-0.5} \ M_\odot / L_{\odot, V}$, consistent with the stellar population estimate for DF2. We use tidal analysis to find that the low-mass models are consistent with the undisturbed isophotes of DF2. Finally we compare with Local Group dwarf galaxies and demonstrate that DF2 is an outlier in both its spatial extent and its relative dark matter deficit.

A large population of ultra-diffuse galaxies (UDGs) was recently discovered in the Coma cluster. Here we present optical spectra of three such UDGs, DF7, DF44 and DF17, which have central surface brightnesses of $\mu_g \approx 24.4-25.1$ mag arcsec$^{-2}$. The spectra were acquired as part of an ancillary program within the SDSS-IV MaNGA Survey. We stacked 19 fibers in the central regions from larger integral field units (IFUs) per source. With over 13.5 hours of on-source integration we achieved a mean signal-to-noise ratio (S/N) in the optical of $9.5$\AA$^{-1}$, $7.9$\AA$^{-1}$ and $5.0$\AA$^{-1}$, respectively, for DF7, DF44 and DF17. Stellar population models applied to these spectra enable measurements of recession velocities, ages and metallicities. The recession velocities of DF7, DF44 and DF17 are $6599^{+40}_{-25}$km/s, $6402^{+41}_{-39}$km/s and $8315^{+43}_{-43}$km/s, spectroscopically confirming that all of them reside in the Coma cluster. The stellar populations of these three galaxies are old and metal-poor, with ages of $7.9^{+3.6}_{-2.5}$Gyr, $8.9^{+4.3}_{-3.3}$Gyr and $9.1^{+3.9}_{-5.5}$Gyr, and iron abundances of $\mathrm{[Fe/H]}$ $-1.0^{+0.3}_{-0.4}$, $-1.3^{+0.4}_{-0.4}$ and $-0.8^{+0.5}_{-0.5}$, respectively. Their stellar masses are $3$-$6\times10^8 M_\odot$. The UDGs in our sample are as old or older than galaxies at similar stellar mass or velocity dispersion (only DF44 has an independently measured dispersion). They all follow the well-established stellar mass$-$stellar metallicity relation, while DF44 lies below the velocity dispersion-metallicity relation. These results, combined with the fact that UDGs are unusually large for their stellar mass, suggest that stellar mass plays a more important role in setting stellar population properties for these galaxies than either size or surface brightness.

Brightest Cluster Galaxies (BCGs) in massive dark matter halos are shaped by complex merging processes. We present a detailed stellar population analysis in the central region of Abell 3827 at $z\sim0.1$, including five-nucleus galaxies involved in a BCG assembly. Based on deep spectroscopy from Multi Unit Spectroscopic Explorer (MUSE), we fit the optical spectra of 13 early-type galaxies (ETGs) in the central 70 kpc of the cluster. The stellar populations in the central $R=1$ kpc of these ETGs are old (6-10 Gyr). Their [Fe/H] increases with $\sigma_{\star}$ and stellar mass. More importantly, [$\alpha$/Fe] of galaxies close to the cluster center do not seem to depend on $\sigma_{\star}$ or stellar mass, indicating that the cluster center shapes the [$\alpha$/Fe]-$\sigma_{\star}$ and [$\alpha$/Fe]-$M_{\star}$ relations differently than other environments where [$\alpha$/Fe] is observed to increase with increasing $\sigma_{\star}$ or $M_{\star}$. Our results reveal the coordinated assembly of BCGs: their building blocks are different from the general low mass populations by their high [$\alpha$/Fe]. Massive galaxies thus grow by accreting preferentially high [$\alpha$/Fe] systems. The radial profiles also bear the imprint of the coordinated assembly. Their declining [Fe/H] and flat [$\alpha$/Fe] radial profiles confirm that the accreted systems have low metallicity and high [$\alpha$/Fe] stellar contents.

The relation between galaxies and dark matter halos reflects the combined effects of many distinct physical processes. Observations indicate that the $z=0$ stellar mass-halo mass (SMHM) relation has remarkably small scatter in stellar mass at fixed halo mass ($\lesssim$ 0.2 dex) with little dependence on halo mass. We investigate the origins of this scatter by combining N-body simulations with observational constraints on the SMHM relation. We find that at the group and cluster scale ($M_{\rm vir}>10^{14}{M_\odot}$) the scatter due purely to hierarchical assembly is $\approx0.16$ dex, which is comparable to recent direct observational estimates. At lower masses, mass buildup since $z\approx2$ is driven largely by in-situ growth. We include a model for the in-situ buildup of stellar mass and find that an intrinsic scatter in this growth channel of $0.2$ dex produces a relation between scatter and halo mass that is consistent with observations from $10^{12}{M_\odot}<M_{\rm vir}<10^{14.75}{M_\odot}$. The approximately constant scatter across a wide range of halo masses at $z=0$ thus appears to be a coincidence as it is determined largely by in-situ growth at low masses and by hierarchical assembly at high masses. These results indicate that the scatter in the SMHM relation can provide unique insight into the regularity of the galaxy formation process.

Stellar photometric variability offers a novel probe of the interior structure and evolutionary state of stars. Here we present a census of stellar variability on day to decade timescales across the color-magnitude diagram for 73,000 stars brighter than $M_I$=-5 in the Whirlpool Galaxy (M51). Our Cycle 24 HST program acquired V and I-band images over 34 epochs spanning one year with pseudo-random cadences enabling sensitivity to periods from days to months. We supplement these data with archival V and I-band HST data obtained in 1995 and 2005 providing sensitivity to variability on decade timescales. At least 50% of stars brighter than $M_I$=-7 show strong evidence for variability within our Cycle 24 data; among stars with V-I>2 the variability fraction rises to ~100%. Large amplitude variability (>0.3 mag) on decade timescales is restricted to red supergiants and very luminous blue stars. Both populations display fairly smooth variability on month-year timescales. The Cepheid instability strip is clearly visible in our data, although the variability fraction within this region never exceeds ~10%. The location of variable stars across the color magnitude diagram broadly agrees with theoretical sources of variability, including the instability strip, red supergiant pulsational instabilities, long-period fundamental mode pulsations, and radiation-dominated envelopes in massive stars. Our data can be used to place stringent constraints on the precise onset of these various instabilities and their lifetimes and growth rates.

We present HST/ACS imaging of twenty-three very low surface brightness ($\mu_{e,V}$=25-27.5) galaxies detected in the fields of four nearby galaxy groups. These objects were selected from deep optical imaging obtained with the Dragonfly Telephoto Array. Seven are new discoveries, while most of the others had been identified previously in visual surveys of deep photographic plates and more recent surveys. Few have previously been studied in detail. From the ACS images, we measure distances to the galaxies using both the tip of the red giant branch method and the surface brightness fluctuations method. We demonstrate that the two methods are consistent with each other in the regime where both can be applied. The distances to 15 out of 20 galaxies with stable measurements are consistent with that of the targeted group within errors. This suggests that assuming group membership based solely on projected proximity is ~75% successful in this regime. The galaxies are nearly round, with a median axis ratio of 0.85, and visually resemble dwarf spheroidal galaxies. The objects have a range of sizes, from $R_e$=0.4 kpc to $R_e$=1.8 kpc, with a median of 1.0 kpc. They range in luminosity from $M_V$=-11.4 to $M_V$=-15.6, with a median of -12.4. Galaxies with $R_e$~1 kpc and $M_V$~-12 are fairly rare in the Local Group but we find many of them in this relatively small sample. Four of the objects fall in the class of ultra diffuse galaxies (UDGs), with $R_e$>1.5 kpc and $\mu_{0,V}$>24 mag/arcsec^2, including the recently identified dark matter deficient galaxy NGC1052-DF2.

The Hyades, Praesepe, and Pleiades are well studied stellar clusters that anchor important secondary stellar age indicators. Recent studies have shown that main sequence turn off-based ages for these clusters may depend on the degree of rotation in the underlying stellar models. Rotation induces structural instabilities that can enhance the chemical mixing of a star, extending its fuel supply. In addition, rotation introduces a modulation of the star's observed magnitude and color due to the effects of gravity darkening. We aim to investigate the extent to which stellar rotation affects the age determination of star clusters. We utilize the MESA stellar evolution code to create models that cover a range of rotation rates corresponding to $\Omega/\Omega_c=0.0$ to $0.6$ in $0.1$ dex steps, allowing the assessment of variations in this dimension. The statistical analysis package, MATCH, is employed to derive ages and metallicities by fitting our MESA models to Tycho $B_T$, $V_T$ and 2MASS $J$, $K_s$ color-magnitude diagrams. We find that the derived ages are relatively insensitive to the effects of rotation. For the Hyades, Praesepe, and Pleiades, we derive ages based on synthetic populations that model a distribution of rotation rates or a fixed rate. Across each case, derived ages tend to agree roughly within errors, near $680$, $590$, and $110-160$ Myr for the Hyades, Praesepe, and Pleiades, respectively. These ages are in agreement with Li depletion boundary-based ages and previous analyses that used non-rotating isochrones. Our methods do not provide a strong constraint on the metallicities of these clusters.

We use the framework developed as part of the MESA Isochrones and Stellar Tracks (MIST) project to assess the utility of several types of observables in jointly measuring the age and 1D stellar model parameters in star clusters. We begin with a pedagogical overview summarizing the effects of stellar model parameters, such as the helium abundance, mass-loss efficiency, and the mixing length parameter, on observational diagnostics including the color-magnitude diagram, mass-radius relation, and surface abundances, amongst others. We find that these parameters and the stellar age influence observables in qualitatively distinctive, degeneracy-breaking ways. To assess the current state of affairs, we use the recent Gaia Data Release 2 (DR2) along with data from the literature to investigate three well-studied old open clusters---NGC6819, M67, NGC6791---as case studies. Although there is no obvious tension between the existing observations and the MIST models for NGC6819, there are interesting discrepancies in the cases of M67 and NGC6791. At this time, parallax zero point uncertainties in Gaia DR2 remain one of the limiting factors in the analysis of these clusters. With a combination of exquisite photometry, parallax distances, and cluster memberships from Gaia at the end of its mission, we anticipate precise and accurate ages for these and other star clusters in the Galaxy.

We present a color-based method for identifying red giants based on Pan-STARRS grz and WISE W1 and W2 photometry. We utilize a subsample of bright stars with precise parallaxes from Gaia DR2 to verify that the color-based selection reliably separates dwarfs from giants. The selection is conservative in the sense that contamination is small (~30%) but not all giants are included (the selection primarily identifies K giants). The color-based selection can be applied to stars brighter than $W1\approx16$, more than two magnitudes fainter than techniques relying on shallower 2MASS photometry. Many streams and clouds are visible in the resulting sky maps, especially when binned by Gaia DR2 proper motions, including the Sagittarius stream, the Hercules-Aquila Cloud, the Eastern Banded Structure, Monoceros, and the Virgo Overdensity. In addition to the characterization of new and known stellar streams, we expect that this method for selecting red giants will enable detailed analysis of the diffuse stellar halo to distances exceeding 100 kpc.

We use deep $HST$ WFC3/IR imaging to study the Initial Mass Function (IMF) of the ultra faint dwarf galaxy Coma Berenices (Com Ber). Our observations reach the lowest stellar mass ever probed in a resolved galaxy, with 50\% completeness at $\sim 0.17$ M$_{\odot}$. Unresolved background galaxies however limit our purity below $\sim 0.23$ M$_{\odot}$. If modeled with a single power law, we find that the IMF slope is $-1.45^{+0.29}_{-0.3}$ (68\% credible intervals), compared to a Milky Way value of $-2.3$. For a broken power law, we obtain a low-mass slope of $-1.18_{-0.33}^{+0.49}$, a high-mass slope of $-1.88_{-0.49}^{+0.43}$ and a break mass of $0.57_{-0.08}^{+0.12}$ M$_{\odot}$, compared to $-1.3$, $-2.3$ and 0.5 M$_{\odot}$ for a Kroupa IMF. For a log-normal IMF model we obtain values of $0.33_{-0.16}^{+0.15}$ M$_{\odot}$ for the location parameter and of $0.68_{-0.12}^{+0.17}$ for $\sigma$ (0.22 M$_{\odot}$ and 0.57 for the Chabrier system IMF). All three parametrizations produce similar agreement with the data. Our results agree with previous analysis of shallower optical HST data. However analysis of similar optical data of other dwarfs finds IMFs significantly more bottom-light than in the Milky Way. These results suggest two, non mutually exclusive, possibilities: that the discrepancy of the dwarf galaxies IMF with respect to the Milky Way is, at least partly, an artifact of using a single power law model, and that there is real variance in the IMF at low masses between the currently studied nearby dwarfs, with Com Ber being similar to the Milky Way, but other dwarfs differing significantly.

We recently published velocity measurements of luminous globular clusters in the galaxy NGC1052-DF2, concluding that it lies far off the canonical stellar mass - halo mass relation. Here we present a revised velocity for one of the globular clusters, GC-98, and a revised velocity dispersion measurement for the galaxy. We find that the intrinsic dispersion $\sigma=5.6^{+5.2}_{-3.8}$ km/s using Approximate Bayesian Computation, or $\sigma=7.8^{+5.2}_{-2.2}$ km/s using the likelihood. The expected dispersion from the stars alone is ~7 km/s. Responding to a request from the Editors of ApJ Letters and RNAAS, we also briefly comment on the recent analysis of our measurements by Martin et al. (2018).

In this paper, we provide a physical model for the origin of variations in the shapes and bump strengths of dust attenuation laws in galaxies by combining a large suite of cosmological "zoom-in" galaxy formation simulations with 3D Monte Carlo dust radiative transfer calculations. We model galaxies over 3 orders of magnitude in stellar mass, ranging from Milky Way like systems through massive galaxies at high-redshift. Critically, for these calculations we employ a constant underlying dust extinction law in all cases, and examine how the role of geometry and radiative transfer effects impact the resultant attenuation curves. Our main results follow. Despite our usage of a constant dust extinction curve, we find dramatic variations in the derived attenuation laws. The slopes of normalized attenuation laws depend primarily on the complexities of star-dust geometry. Increasing fractions of unobscured young stars flatten normalized curves, while increasing fractions of unobscured old stars steepen curves. Similar to the slopes of our model attenuation laws, we find dramatic variation in the 2175 Angstrom ultraviolet (UV) bump strength, including a subset of curves with little to no bump. These bump strengths are primarily influenced by the fraction of unobscured O and B stars in our model, with the impact of scattered light having only a secondary effect. Taken together, these results lead to a natural relationship between the attenuation curve slope and 2175 Angstrom bump strength. Finally, we apply these results to a 25 Mpc/h box cosmological hydrodynamic simulation in order to model the expected dispersion in attenuation laws at integer redshifts from z=0-6. A significant dispersion is expected at low redshifts, and decreases toward z=6. We provide tabulated results for the best fit median attenuation curve at all redshifts.

Oxygen is the most abundant "metal" element in stars and in the cosmos. But determining oxygen abundances in stars has proven challenging, because of the shortage of detectable atomic oxygen lines in their optical spectra as well as observational and theoretical complications with these lines (e.g., blends, 3D, non-LTE). Nonetheless, Ting et al. (2017) were recently able to demonstrate that oxygen abundances can be determined from low-resolution (R$\simeq$2000) optical spectra. Here we investigate the physical processes that enable such a measurement for cool stars, such as K-giants. We show that the strongest spectral diagnostics of oxygen come from the CNO atomic-molecular network, but are manifested in spectral features that do not involve oxygen. In the outer atmosphere layers most of the carbon is locked up in CO, and changes to the oxygen abundance directly affect the abundances of all other carbon-bearing molecules, thereby changing the strength of CH, CN, and C$_2$ features across the optical spectrum. In deeper atmosphere layers most of the carbon is in atomic form, and any change in the oxygen abundance has little effect on the other carbon-bearing molecules. The key physical effect enabling such oxygen abundance measurements is that spectral features in the optical arise from both the CO-dominant and the atomic carbon-dominant regions, providing non-degenerate constraints on both C and O. Beyond the case at hand, the results show that physically sound abundances measurements need not be limited to those elements that have observable lines themselves.

The effective temperature (Teff) distribution of stellar evolution models along the red giant branch (RGB) is sensitive to a number of parameters including the overall metallicity, elemental abundance patterns, the efficiency of convection, and the treatment of the surface boundary condition. Recently there has been interest in using observational estimates of the RGB Teff to place constraints on the mixing length parameter, a_MLT, and possible variation with metallicity. Here we use 1D MESA stellar evolution models to explore the sensitivity of the RGB Teff to the treatment of the surface boundary condition. We find that different surface boundary conditions can lead to +/- 100 K metallicity-dependent offsets on the RGB relative to one another in spite of the fact that all models can reproduce the properties of the Sun. Moreover, for a given atmosphere T-tau relation, we find that the RGB Teff is also sensitive to the optical depth at which the surface boundary condition is applied in the stellar model. Nearly all models adopt the photosphere as the location of the surface boundary condition but this choice is somewhat arbitrary. We compare our models to stellar parameters derived from the APOGEE-Kepler sample of first ascent red giants and find that systematic uncertainties in the models due to treatment of the surface boundary condition place a limit of ~100 K below which it is not possible to make firm conclusions regarding the fidelity of the current generation of stellar models.

We report tentative evidence for a cold stellar stream in the ultra-diffuse galaxy NGC1052-DF2. If confirmed, this stream (which we refer to as "The Maybe Stream") would be the first cold stellar stream detected outside of the Local Group. The candidate stream is very narrow and has an unusual and highly curved shape.

We recently found an ultra diffuse galaxy (UDG) with a half-light radius of R_e = 2.2 kpc and little or no dark matter. The total mass of NGC1052-DF2 was measured from the radial velocities of bright compact objects that are associated with the galaxy. Here we analyze these objects using a combination of HST imaging and Keck spectroscopy. Their average size is <r_h> = 6.2+-0.5 pc and their average ellipticity is <{\epsilon}> = 0.18+-0.02. From a stacked Keck spectrum we derive an age >9 Gyr and a metallicity of [Fe/H] = -1.35+-0.12. Their properties are similar to {\omega} Centauri, the brightest and largest globular cluster in the Milky Way, and our results demonstrate that the luminosity function of metal-poor globular clusters is not universal. The fraction of the total stellar mass that is in the globular cluster system is similar to that in other UDGs, and consistent with "failed galaxy" scenarios where star formation terminated shortly after the clusters were formed. However, the galaxy is a factor of ~1000 removed from the relation between globular cluster mass and total galaxy mass that has been found for other galaxies, including other UDGs. We infer that a dark matter halo is not a prerequisite for the formation of metal-poor globular cluster-like objects in high redshift galaxies.

Studies of galaxy surveys in the context of the cold dark matter paradigm have shown that the mass of the dark matter halo and the total stellar mass are coupled through a function that varies smoothly with mass. Their average ratio M_{halo}/M_{stars} has a minimum of about 30 for galaxies with stellar masses near that of the Milky Way (approximately 5x10^{10} solar masses) and increases both towards lower masses and towards higher masses. The scatter in this relation is not well known; it is generally thought to be less than a factor of two for massive galaxies but much larger for dwarf galaxies. Here we report the radial velocities of ten luminous globular-cluster-like objects in the ultra-diffuse galaxy NGC1052-DF2, which has a stellar mass of approximately 2x10^8 solar masses. We infer that its velocity dispersion is less than 10.5 kilometers per second with 90 per cent confidence, and we determine from this that its total mass within a radius of 7.6 kiloparsecs is less than 3.4x10^8 solar masses. This implies that the ratio M_{halo}/M_{stars} is of order unity (and consistent with zero), a factor of at least 400 lower than expected. NGC1052-DF2 demonstrates that dark matter is not always coupled with baryonic matter on galactic scales.

Forward modeling of the full galaxy SED is a powerful technique, providing self-consistent constraints on stellar ages, dust properties, and metallicities. However, the accuracy of these results is contingent on the accuracy of the model. One significant source of uncertainty is the contribution of obscured AGN, as they are relatively common and can produce substantial mid-IR (MIR) emission. Here we include emission from dusty AGN torii in the Prospector SED-fitting framework, and fit the UV-IR broadband photometry of 129 nearby galaxies. We find that 10% of the fitted galaxies host an AGN contributing $>$10% of the observed galaxy MIR luminosity. We demonstrate the necessity of this AGN component in the following ways. First, we compare observed spectral features to spectral features predicted from our model fit to the photometry. We find that the AGN component greatly improves predictions for observed H$\alpha$ and H$\beta$ luminosities, as well as mid-infrared Akari and Spitzer/IRS spectra. Second, we show that inclusion of the AGN component changes stellar ages and SFRs by up to a factor of 10, and dust attenuations by up to a factor of 2.5. Finally, we show that the strength of our model AGN component correlates with independent AGN indicators, suggesting that these galaxies truly host AGN. Notably, only 46% of the SED-detected AGN would be detected with a simple MIR color selection. Based on these results, we conclude that SED models which fit MIR data without AGN components are vulnerable to substantial bias in their derived parameters.

We discuss the approach of searching low mass dwarf galaxies, $\lesssim10^6\textrm{ M}_{\odot}$, in the general field, using integrated light surveys. By exploring the limiting surface brightness-spatial resolution ($\mu_{\textrm{eff,lim}}-\theta$) parameter space, we suggest that faint field dwarfs in the Local Volume, between $3$ and $10 \textrm{ Mpc}$, are expected to be detected effectively and in large numbers using integrated light photometric surveys, complementary to the classical star counts method. We use a sample of Local Group dwarf galaxies to construct relations between their photometric and structural parameters, $\textrm{M}_{*}$-$\mu_{\textrm{eff,V}}$ and $\textrm{M}_{*}$-$\textrm{R}_{\textrm{eff}}$. We use these relations, along with assumed functional forms for the halo mass function and the stellar mass-halo mass relation, to calculate the lowest detectable stellar masses in the Local Volume and the expected number of galaxies as a function of the limiting surface brightness and spatial resolution. The number of detected galaxies depends mostly on the limiting surface brightness for distances $>3\textrm{ Mpc}$ while spatial resolution starts to play a role at distances $>8\textrm{ Mpc}$. Surveys with $\mu_{\textrm{eff,lim}}\sim30\textrm{ mag arcsec}^{-2}$ should be able to detect galaxies with stellar masses down to $\sim10^4 \textrm{ M}_{\odot}$ in the Local Volume. Depending on the assumed stellar mass-halo mass relation, the expected number of galaxies between $3$ and $10\textrm{ Mpc}$ is $0.04-0.35\textrm{ deg}^{-2}$, assuming a limiting surface brightness of $\sim29-30\textrm{ mag arcsec}^{-2}$ and a spatial resolution $<4''$. We currently look for field dwarf galaxies by performing a blank wide-field survey with the Dragonfly Telephoto Array, optimized for the detection of ultra-low surface brightness structures.

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 develop a data-driven spectral model for identifying and characterizing spatially unresolved multiple-star systems and apply it to APOGEE DR13 spectra of main-sequence stars. Binaries and triples are identified as targets whose spectra can be significantly better fit by a superposition of two or three model spectra, drawn from the same isochrone, than any single-star model. From an initial sample of $\sim$20,000 main-sequence targets, we identify $\sim$2,500 binaries in which both the primary and secondary star contribute detectably to the spectrum, simultaneously fitting for the velocities and stellar parameters of both components. We additionally identify and fit $\sim$200 triple systems, as well as $\sim$700 velocity-variable systems in which the secondary does not contribute detectably to the spectrum. Our model simplifies the process of simultaneously fitting single- or multi-epoch spectra with composite models and does not depend on a velocity offset between the two components of a binary, making it sensitive to traditionally undetectable systems with periods of hundreds or thousands of years. In agreement with conventional expectations, almost all the spectrally-identified binaries with measured parallaxes fall above the main sequence in the color-magnitude diagram. We find excellent agreement between spectrally and dynamically inferred mass ratios for the $\sim$600 binaries in which a dynamical mass ratio can be measured from multi-epoch radial velocities. We obtain full orbital solutions for 64 systems, including 14 close binaries within hierarchical triples. We make available catalogs of stellar parameters, abundances, mass ratios, and orbital parameters.

We investigate the evolution of dust content in galaxies from redshifts z=0 to z=9.5. Using empirically motivated prescriptions, we model galactic-scale properties -- including halo mass, stellar mass, star formation rate, gas mass, and metallicity -- to make predictions for the galactic evolution of dust mass and dust temperature in main sequence galaxies. Our simple analytic model, which predicts that galaxies in the early Universe had greater quantities of dust than their low-redshift counterparts, does a good job at reproducing observed trends between galaxy dust and stellar mass out to z~6. We find that for fixed galaxy stellar mass, the dust temperature increases from z=0 to z=6. Our model forecasts a population of low-mass, high-redshift galaxies with interstellar dust as hot as, or hotter than, their more massive counterparts; but this prediction needs to be constrained by observations. Finally, we make predictions for observing 1.1-mm flux density arising from interstellar dust emission with the Atacama Large Millimeter Array.

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.

SDSS-V will be an all-sky, multi-epoch spectroscopic survey of over six million objects. It is designed to decode the history of the Milky Way, trace the emergence of the chemical elements, reveal the inner workings of stars, and investigate the origin of planets. It will also create an integral-field spectroscopic map of the gas in the Galaxy and the Local Group that is 1,000x larger than the current state of the art and at high enough spatial resolution to reveal the self-regulation mechanisms of galactic ecosystems. SDSS-V will pioneer systematic, spectroscopic monitoring across the whole sky, revealing changes on timescales from 20 minutes to 20 years. The survey will thus track the flickers, flares, and radical transformations of the most luminous persistent objects in the universe: massive black holes growing at the centers of galaxies. The scope and flexibility of SDSS-V will be unique among extant and future spectroscopic surveys: it is all-sky, with matched survey infrastructures in both hemispheres; it provides near-IR and optical multi-object fiber spectroscopy that is rapidly reconfigurable to serve high target densities, targets of opportunity, and time-domain monitoring; and it provides optical, ultra-wide-field integral field spectroscopy. SDSS-V, with its programs anticipated to start in 2020, will be well-timed to multiply the scientific output from major space missions (e.g., TESS, Gaia, eROSITA) and ground-based projects. SDSS-V builds on the 25-year heritage of SDSS's advances in data analysis, collaboration infrastructure, and product deliverables. The project is now refining its science scope, optimizing the survey strategies, and developing new hardware that builds on the SDSS-IV infrastructure. We present here an overview of the current state of these developments as we seek to build our worldwide consortium of institutional and individual members.

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.

The observable spectrum of an unresolved binary star system is a superposition of two single-star spectra. Even without a detectable velocity offset between the two stellar components, the combined spectrum of a binary system is in general different from that of either component, and fitting it with single-star models may yield inaccurate stellar parameters and abundances. We perform simple experiments with synthetic spectra to investigate the effect of unresolved main-sequence binaries on spectral fitting, modeling spectra similar to those collected by the APOGEE, GALAH, and LAMOST surveys. We find that fitting unresolved binaries with single-star models introduces systematic biases in the derived stellar parameters and abundances that are modest but certainly not negligible, with typical systematic errors of $300\,\rm K$ in $T_{\rm eff}$, 0.1 dex in $\log g$, and 0.1 dex in $[\rm Fe/H]$ for APOGEE-like spectra of solar-type stars. These biases are smaller for spectra at optical wavelengths than in the near-infrared. We show that biases can be corrected by fitting spectra with a binary model, which adds only two labels to the fit and includes single-star models as a special case. Our model provides a promising new method to constrain the Galactic binary population, including systems with single-epoch spectra and no detectable velocity offset between the two stars.

We present the properties of NGC 4993, the host galaxy of GW170817, the first gravitational wave (GW) event from the merger of a binary neutron star (BNS) system and the first with an electromagnetic (EM) counterpart. We use both archival photometry and new optical/near-IR imaging and spectroscopy, together with stellar population synthesis models to infer the global properties of the host galaxy. We infer a star formation history peaked at $\gtrsim 10$ Gyr ago, with subsequent exponential decline leading to a low current star formation rate of 0.01 M$_{\odot}$ yr$^{-1}$, which we convert into a binary merger timescale probability distribution. We find a median merger timescale of $11.2^{+0.7}_{-1.4}$ Gyr, with a 90% confidence range of $6.8-13.6$ Gyr. This in turn indicates an initial binary separation of $\approx 4.5$ R$_{\odot}$, comparable to the inferred values for Galactic BNS systems. We also use new and archival $Hubble$ $Space$ $Telescope$ images to measure a projected offset of the optical counterpart of $2.1$ kpc (0.64$r_{e}$) from the center of NGC 4993 and to place a limit of $M_{r} \gtrsim -7.2$ mag on any pre-existing emission, which rules out the brighter half of the globular cluster luminosity function. Finally, the age and offset of the system indicates it experienced a modest natal kick with an upper limit of $\sim 200$ km s$^{-1}$. Future GW$-$EM observations of BNS mergers will enable measurement of their population delay time distribution, which will directly inform their viability as the dominant source of $r$-process enrichment in the Universe.

The LAMOST survey has acquired low-resolution spectra (R=1,800) for 5 million stars across the Milky Way, far more than any current stellar survey at a corresponding or higher spectral resolution. It is often assumed that only very few elemental abundances can be measured from such low-resolution spectra, limiting their utility for Galactic archaeology studies. However, Ting et al. (2017) used ab initio models to argue that low-resolution spectra should enable precision measurements of many elemental abundances, at least in theory. Here we verify this claim in practice by measuring the relative abundances of 14 elements from LAMOST spectra with a precision of $\lesssim$ 0.1 dex for objects with ${\rm S/N}_{\rm LAMOST}$ > 30 (per pixel). We employ a spectral modeling method in which a data-driven model is combined with priors that the model gradient spectra should resemble ab initio spectral models. This approach assures that the data-driven abundance determinations draw on physically sensible features in the spectrum in their predictions and do not just exploit astrophysical correlations among abundances. Our analysis is constrained to the number of elemental abundances measured in the APOGEE survey, which is the source of the training labels. Obtaining high quality/resolution spectra for a subset of LAMOST stars to measure more elemental abundances as training labels and then applying this method to the full LAMOST catalog will provide a sample with more than 20 elemental abundances that is an order of magnitude larger than current high-resolution surveys, substantially increasing the sample size for Galactic archaeology.

We present Gemini/GNIRS cross-dispersed near-infrared spectra of 12 nearby early-type galaxies, with the aim of testing commonly used stellar population synthesis models. We select a subset of galaxies from the atlas/ sample which span a wide range of ages (SSP-equivalent ages of 1--15~Gyr) at approximately solar metallicity. We derive star formation histories using four different stellar population synthesis models, namely those of [bruzual_stellar_2003], Conroy, Gunn \& White (2009; 2010), [maraston_stellar_2011] and [vazdekis_uv-extended_2016]. We compare star formation histories derived from near-infrared spectra with those derived from optical spectra using the same models. We find that while all models agree in the optical, the derived star formation histories vary dramatically from model to model in the near-infrared. We find that this variation is largely driven by the choice of stellar spectral library, such that models including high quality spectral libraries provide the best fits to the data, and are the most self-consistent when comparing optically-derived properties with near-infrared ones. We also find the impact of age variation in the near-infrared to be subtle, and largely encoded in the shape of the continuum, meaning that the common approach of removing continuum information with a high-order polynomial greatly reduces our ability to constrain ages in the near-infrared.

We analyze the optical continuum of star-forming galaxies in SDSS by fitting stacked spectra with stellar population synthesis models to investigate the relation between stellar mass, stellar metallicity, dust attenuation and star formation rate. We fit models calculated with star formation and chemical evolution histories that are derived empirically from multi-epoch observations of the stellar mass---star formation rate and the stellar mass---gas-phase metallicity relations, respectively. We also fit linear combinations of single burst models with a range of metallicities and ages. Star formation and chemical evolution histories are unconstrained for these models. The stellar mass---stellar metallicity relations obtained from the two methods agree with the relation measured from individual supergiant stars in nearby galaxies. These relations are also consistent with the relation obtained from emission line analysis of gas-phase metallicity after accounting for systematic offsets in the gas-phase-metallicity. We measure dust attenuation of the stellar continuum and show that its dependence on stellar mass and star formation rate is consistent with previously reported results derived from nebular emission lines. However, stellar continuum attenuation is smaller than nebular emission line attenuation. The continuum-to-nebular attenuation ratio depends on stellar mass and is smaller in more massive galaxies. Our consistent analysis of stellar continuum and nebular emission lines paves the way for a comprehensive investigation of stellar metallicities of star-forming and quiescent galaxies.

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 Hubble Space Telescope imaging of two ultra diffuse galaxies (UDGs) with measured stellar velocity dispersions in the Coma cluster. The galaxies, Dragonfly 44 and DFX1, have effective radii of 4.7 kpc and 3.5 kpc and velocity dispersions of $47^{+8}_{-6}$ km/s and $30^{+7}_{-7}$ km/s, respectively. Both galaxies are associated with a striking number of compact objects, tentatively identified as globular clusters: $N_{\rm gc}=74\pm 18$ for Dragonfly 44 and $N_{\rm gc}=62\pm 17$ for DFX1. The number of globular clusters is far higher than expected from the luminosities of the galaxies but is consistent with expectations from the empirical relation between dynamical mass and globular cluster count defined by other galaxies. Combining our data for these two objects with previous HST observations of Coma UDGs we find that UDGs have a factor of $6.9^{+1.0}_{-2.4}$ more globular clusters than other galaxies of the same luminosity, in contrast to a recent study of a similar sample by Amorisco et al. (2017), but consistent with earlier results for individual galaxies. The Harris et al. (2017) relation between globular cluster count and dark matter halo mass implies a median halo mass of $M_{\rm halo}\sim 1.5\times 10^{11}\,{\rm M}_{\odot}$ for the sixteen Coma UDGs that have been observed with HST so far, with the largest and brightest having $M_{\rm halo}\sim 5\times 10^{11}\,{\rm M}_{\odot}$.

Understanding the evolution of the Milky Way calls for the precise abundance determination of many elements in many stars. A common perception is that deriving more than a few elemental abundances ([Fe/H], [$\alpha$/Fe], perhaps [C/H], [N/H]) requires medium-to-high spectral resolution, $R \gtrsim 10,000$, mostly to overcome the effects of line blending. In recent work (Rix et al. 2016; Ting et al. 2016) we presented an efficient and practical way to model the full stellar spectrum, even when fitting a large number of stellar labels simultaneously. In this paper we quantify to what precision the abundances of many different elements can be recovered, as a function of spectroscopic resolution and wavelength range. In the limit of perfect spectral models and spectral normalization, we show that the precision of elemental abundances is nearly independent of resolution, for a fixed exposure time and number of detector pixels; low-resolution spectra simply afford much higher S/N per pixel and generally larger wavelength range in a single setting. We also show that estimates of most stellar labels are not strongly correlated with one another once $R \gtrsim 1,000$. Modest errors in the line spread function, as well as small radial velocity errors, do not affect these conclusions, and data driven models indicate that spectral (continuum) normalization can be achieved well enough in practice. These results, to be confirmed with an analysis of observed low-resolution data, open up new possibilities for the design of large spectroscopic stellar surveys and for the re-analysis of archival low-resolution datasets.

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.

We utilise a series of high-resolution cosmological zoom simulations of galaxy formation to investigate the relationship between the ultraviolet (UV) slope, beta, and the ratio of the infrared luminosity to UV luminosity (IRX) in the spectral energy distributions (SEDs) of galaxies. We employ dust radiative transfer calculations in which the SEDs of the stars in galaxies propagate through the dusty interstellar medium. Our main goals are to understand the origin of, and scatter in the IRX-beta relation; to assess the efficacy of simplified stellar population synthesis screen models in capturing the essential physics in the IRX-beta relation; and to understand systematic deviations from the canonical local IRX-beta relations in particular populations of high-redshift galaxies. Our main results follow. Galaxies that have young stellar populations with relatively cospatial UV and IR emitting regions and a Milky Way-like extinction curve fall on or near the standard Meurer relation. This behaviour is well captured by simplified screen models. Scatter in the IRX-beta relation is dominated by three major effects: (i) older stellar populations drive galaxies below the relations defined for local starbursts due to a reddening of their intrinsic UV SEDs; (ii) complex geometries in high-z heavily star forming galaxies drive galaxies toward blue UV slopes owing to optically thin UV sightlines; (iii) shallow extinction curves drive galaxies downward in the IRX-beta plane due to lowered NUV/FUV extinction ratios. We use these features of the UV slopes of galaxies to derive a fitting relation that reasonably collapses the scatter back toward the canonical local relation. Finally, we use these results to develop an understanding for the location of two particularly enigmatic populations of galaxies in the IRX-beta plane: z~2-4 dusty star forming galaxies, and z>5 star forming galaxies.

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 use the first Gaia data release, combined with RAVE and APOGEE spectroscopic surveys, to investigate the origin of halo stars within <~3 kpc from the Sun. We identify halo stars kinematically, as moving with a relative speed of at least 220 km/s with respect to the local standard of rest. These stars are in general more metal-poor than the disk, but surprisingly, half of our halo sample is comprised of stars with [Fe/H]>-1. The orbital directions of these metal-rich halo stars are preferentially aligned with the disk rotation, in sharp contrast with the isotropic orbital distribution of the more metal-poor halo stars. We find similar properties in the Latte cosmological zoom-in simulation of a Milky Way-like galaxy from the FIRE project. In Latte, metal-rich halo stars formed primarily inside of the solar circle, while lower-metallicity halo stars preferentially formed at larger distances (extending beyond the virial radius). This suggests that metal-rich halo stars in the Solar neighborhood in fact formed in-situ within the Galactic disk rather than having been accreted from satellite systems. These stars, currently on halo-like orbits, therefore have likely undergone substantial radial migration/heating.

Broadband photometry of galaxies measures an unresolved mix of complex stellar populations, gas, and dust. Interpreting these data is a challenge for models: many studies have shown that properties derived from modeling galaxy photometry are uncertain by a factor of two or more, and yet answering key questions in the field now requires higher accuracy than this. Here, we present a new model framework specifically designed for these complexities. Our model, Prospector-$\alpha$, includes dust attenuation and re-radiation, a flexible attenuation curve, nebular emission, stellar metallicity, and a 6-component nonparametric star formation history. The flexibility and range of the parameter space, coupled with MCMC sampling within the Prospector inference framework, is designed to provide unbiased parameters and realistic error bars. We assess the accuracy of the model with aperture-matched optical spectroscopy, which was excluded from the fits. We compare spectral features predicted solely from fits to the broadband photometry to the observed spectral features. Our model predicts H$\alpha$ luminosities with a scatter of $\sim$0.18 dex and an offset of $\sim$0.1 dex across a wide range of morphological types and stellar masses. This agreement is remarkable, as the H$\alpha$ luminosity is dependent on accurate star formation rates, dust attenuation, and stellar metallicities. The model also accurately predicts dust-sensitive Balmer decrements, spectroscopic stellar metallicities, PAH mass fractions, and the age- and metallicity-sensitive features D$_{\mathrm{n}}$4000 and H$\delta$. Although the model passes all these tests, we caution that we have not yet assessed its performance at higher redshift or the accuracy of recovered stellar masses.

In the era of large stellar spectroscopic surveys, there is emphasis on deriving not only stellar abundances but also ages for millions of stars. In the context of Galactic archeology, stellar ages provide a direct probe of the formation history of the Galaxy. We use the stellar evolution code MESA to compute models with atomic diffusion--with and without radiative acceleration--and extra mixing in the surface layers. The extra mixing consists of both density-dependent turbulent mixing and envelope overshoot mixing. Based on these models we argue that it is important to distinguish between initial, bulk abundances (parameters) and current, surface abundances (variables) in the analysis of individual stellar ages. In stars that maintain radiative regions on evolutionary timescales, atomic diffusion modifies the surface abundances. We show that when initial, bulk metallicity is equated with current, surface metallicity in isochrone age analysis the resulting stellar ages can be systematically over-estimated by up to 20%. The change of surface abundances with evolutionary phase also complicates chemical tagging, the concept that dispersed star clusters can be identified through unique, high-dimensional chemical signatures. Stars from the same cluster, but in different evolutionary phases, will show different surface abundances. We speculate that calibration of stellar models may allow us to estimate not only stellar ages but also initial abundances for individual stars. In the meantime, analyzing the chemical properties of stars in similar evolutionary phases is essential to minimize the effects of atomic diffusion in the context of chemical tagging.

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 investigate the integrated properties of massive (>10 Msun), rotating, single-star stellar populations for a variety of initial rotation rates (v/vcrit=0.0, 0.2, 0.4, 0.5, and 0.6). We couple the new MESA Isochrone and Stellar Tracks (MIST) models to the Flexible Stellar Population Synthesis (FSPS) package, extending the stellar population synthesis models to include the contributions from very massive stars (>100 Msun), which can be significant in the first ~4 Myr after a starburst. These models predict ionizing luminosities that are consistent with recent observations of young nuclear star clusters. We also construct composite stellar populations assuming a distribution of initial rotation rates. Even in low-metallicity environments where rotation has a significant effect on the evolution of massive stars, we find that stellar population models require a significant contribution from fast-rotating (v/vcrit>0.4) stars in order to sustain the production of ionizing photons beyond a few Myr following a starburst. These results have potentially important implications for cosmic reionization by massive stars and the interpretation of nebular emission lines in high-redshift star-forming galaxies.

Using group catalogs from the SDSS DR7, we attempt to measure galactic conformity in the local universe. We measure the quenched fraction of neighbor galaxies around isolated primary galaxies, dividing the isolated sample into star-forming and quiescent objects. We restrict our measurements to scales $>1$ Mpc to probe the correlations between the formation histories of distinct halos. Over the stellar mass range $10^{9.7} \le M_\ast/M_\odot \le 10^{10.9}$, we find minimal statistical evidence for conformity. We further compare these data to predictions of the halo age-matching model, in which the oldest galaxies are associated with the oldest halos at fixed $M_\ast$. For models with strong correlations between halo and stellar age, the conformity signal is too large to be consistent with the data. For weaker implementations of age-matching, galactic conformity is not a sensitive diagnostic of halo assembly bias, and would not produce a detectable signal in SDSS data. We reproduce the results of Kauffmann et al 2013, in which the star formation rates of neighbor galaxies are significantly reduced around primary galaxies when the primaries are themselves low star formers. However, we find this result is mainly driven by contamination in the isolation criterion, when using our group catalog to remove the small fraction of satellite galaxies in the sample, the conformity signal largely goes away. Lastly, we show that small conformity signals, i.e., 2-5% differences in the quenched fractions of neighbor galaxies, can be produced by mechanisms other than halo assembly bias. For example, if passive galaxies occupy more massive halos than star forming galaxies of the same stellar mass, a conformity signal that is consistent with recent measurements from PRIMUS (Berti et al 2016) can be produced.

The Universe is almost totally unexplored at low surface brightness levels. In spite of great progress in the construction of large telescopes and improvements in the sensitivity of detectors, the limiting surface brightness of imaging observations has remained static for about forty years. Recent technical advances have at last begun to erode the barriers preventing progress. In this Chapter we describe the technical challenges to low surface brightness imaging, describe some solutions, and highlight some relevant observations that have been undertaken recently with both large and small telescopes. Our main focus will be on discoveries made with the Dragonfly Telephoto Array (Dragonfly), which is a new telescope concept designed to probe the Universe down to hitherto unprecedented low surface brightness levels. We conclude by arguing that these discoveries are probably only scratching the surface of interesting phenomena that are observable when the Universe is explored at low surface brightness levels.

Unlike spiral galaxies such as the Milky Way, the majority of the stars in massive elliptical galaxies were formed in a short period early in the history of the Universe. The duration of this formation period can be measured using the ratio of magnesium to iron abundance ([Mg/Fe]), which reflects the relative enrichment by core-collapse and type Ia supernovae. For local galaxies, [Mg/Fe] probes the combined formation history of all stars currently in the galaxy, including younger and metal-poor stars that were added during late-time mergers. Therefore, to directly constrain the initial star-formation period, we must study galaxies at earlier epochs. The most distant galaxy for which [Mg/Fe] had previously been measured is at z~1.4, with [Mg/Fe]=0.45(+0.05,-0.19). A slightly earlier epoch (z~1.6) was probed by stacking the spectra of 24 massive quiescent galaxies, yielding an average [Mg/Fe] of 0.31+/-0.12. However, the relatively low S/N of the data and the use of index analysis techniques for both studies resulted in measurement errors that are too large to allow us to form strong conclusions. Deeper spectra at even earlier epochs in combination with analysis techniques based on full spectral fitting are required to precisely measure the abundance pattern shortly after the major star-forming phase (z>2). Here we report a measurement of [Mg/Fe] for a massive quiescent galaxy at z=2.1. With [Mg/Fe]=0.59+/-0.11, this galaxy is the most Mg-enhanced massive galaxy found so far, having twice the Mg enhancement of similar-mass galaxies today. The abundance pattern of the galaxy is consistent with enrichment exclusively by core-collapse supernovae and with a star-formation timescale of 0.1-0.5 Gyr - characteristics that are similar to population II stars in the Milky Way. With an average past SFR of 600-3000 Msol/yr, this galaxy was among the most vigorous star-forming galaxies in the Universe.

Accounting for nebular emission when modeling galaxy spectral energy distributions (SEDs) is important, as both line and continuum emission can contribute significantly to the total observed flux. In this work, we present a new nebular emission model integrated within the Flexible Stellar Population Synthesis code that computes the total line and continuum emission for complex stellar populations using the photoionization code Cloudy. The self-consistent coupling of the nebular emission to the matched ionizing spectrum produces emission line intensities that correctly scale with the stellar population as a function of age and metallicity. This more complete model of galaxy SEDs will improve estimates of global gas properties derived with diagnostic diagrams, star formation rates based on H$\alpha$, and stellar masses derived from NIR broadband photometry. Our models agree well with results from other photoionization models and are able to reproduce observed emission from H II regions and star-forming galaxies. Our models show improved agreement with the observed H II regions in the Ne III/O II plane and show satisfactory agreement with He II emission from $z=2$ galaxies when including rotating stellar models. Models including post-asymptotic giant branch stars are able to reproduce line ratios consistent with low-ionization emission regions (LIERs).

Large dust grains can fluctuate dramatically in their local density, relative to gas, in neutral, turbulent disks. Small, high-redshift galaxies (before reionization) represent ideal environments for this process. We show via simple arguments and simulations that order-of-magnitude fluctuations are expected in local abundances of large grains under these conditions. This can have important consequences for star formation and stellar abundances in extremely metal-poor stars. Low-mass stars could form in dust-enhanced regions almost immediately after some dust forms, even if the galaxy-average metallicity is too low for fragmentation to occur. The abundances of these 'promoted' stars may contain interesting signatures, as the CNO abundances (concentrated in large carbonaceous grains and ices) and Mg and Si (in large silicate grains) can be enhanced or fluctuate independently. Remarkably, otherwise puzzling abundance patterns of some metal-poor stars can be well-fit by standard core-collapse SNe yields, if we allow for fluctuating dust-to-gas ratios. We also show that the observed log-normal-like distribution of enhancements in these species agrees with our simulations. Moreover, we confirm Mg and Si are correlated in these stars, with abundance ratios similar to those in local silicate grains. Meanwhile [Mg/Ca], predicted to be nearly invariant from pure SNe yields, shows large enhancements as expected in the dust-promoted model, preferentially in the [C/Fe]-enhanced metal-poor stars. This suggests that (1) dust exists in second-generation star formation, (2) dust-to-gas ratio fluctuations occur and can be important for star formation, and (3) light element abundances of these stars may be affected by the chemistry of dust where they formed, rather than directly tracing nucleosynthesis.

There is evidence that the well-established mass-metallicity relation in galaxies is correlated with a third parameter: star formation rate (SFR). The strength of this correlation may be used to disentangle the relative importance of different physical processes (e.g., infall of pristine gas, metal-enriched outflows) in governing chemical evolution. However, all three parameters are susceptible to biases that might affect the observed strength of the relation between them. We analyze possible sources of systematic error, including sample bias, application of signal-to-noise ratio cuts on emission lines, choice of metallicity calibration, uncertainty in stellar mass determination, aperture effects, and dust. We present the first analysis of the relation between stellar mass, gas phase metallicity, and SFR using strong line abundance diagnostics from Dopita et al. (2013) for ~130,000 star-forming galaxies in the Sloan Digital Sky Survey and provide a detailed comparison of these diagnostics in an appendix. Using these abundance diagnostics yields a 30-55% weaker anti-correlation between metallicity and SFR at fixed stellar mass than that reported by Mannucci et al. (2010). We find that, for all abundance diagnostics, the anti-correlation with SFR is stronger for the relatively few galaxies whose current SFRs are elevated above their past average SFRs. This is also true for the new abundance diagnostic of Dopita et al. (2016), which gives anti-correlation between metallicity and SFR only in the high specific star formation rate (sSFR) regime, in contrast to the recent results of Kashino et al. (2016). The poorly constrained strength of the relation between stellar mass, metallicity, and SFR must be carefully accounted for in theoretical studies of chemical evolution.

We describe the redmonster automated redshift measurement and spectral classification software designed for the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digital Sky Survey IV (SDSS-IV). We describe the algorithms, the template standard and requirements, and the newly developed galaxy templates to be used on eBOSS spectra. We present results from testing on early data from eBOSS, where we have found a 90.5% automated redshift and spectral classification success rate for the luminous red galaxy sample (redshifts 0.6$\lesssim z \lesssim$1.0). The redmonster performance meets the eBOSS cosmology requirements for redshift classification and catastrophic failures, and represents a significant improvement over the previous pipeline. We describe the empirical processes used to determine the optimum number of additive polynomial terms in our models and an acceptable $\Delta\chi_r^2$ threshold for declaring statistical confidence. Statistical errors on redshift measurement due to photon shot noise are assessed, and we find typical values of a few tens of km s$^{-1}$. An investigation of redshift differences in repeat observations scaled by error estimates yields a distribution with a Gaussian mean and standard deviation of $\mu\sim$0.01 and $\sigma\sim$0.65, respectively, suggesting the reported statistical redshift uncertainties are over-estimated by $\sim$54%. We assess the effects of object magnitude, signal-to-noise ratio, fiber number, and fiber head location on the pipeline's redshift success rate. Finally, we describe directions of ongoing development.

Long dynamical timescales in the outskirts of galaxies preserve the information content of their accretion histories, for example in the form of stellar population gradients. We present a detailed analysis of the stellar halo properties of a statistically representative sample of early-type galaxies from the Illustris simulation and show that stellar population gradients at large radii can indeed be used to infer basic properties of galactic accretion histories. We measure metallicity, age, and surface-brightness profiles in quiescent Illustris galaxies ranging from $\mathrm{M}_\star = 10^{10} - 2\times 10^{12}\;\mathrm{M}_\odot$ and show that they are in reasonable agreement with observations. At fixed mass, galaxies that accreted little of their stellar halo material tend to have steeper metallicity and surface-brightness profiles between 2 - 4 effective radii (R$_e$) than those with larger accreted fractions. Profiles of metallicity and surface-brightness in the stellar halo typically flatten from z = 1 to the present. This suggests that the accretion of stars into the stellar halo tends to flatten metallicity and surface-brightness profiles, a picture which is supported by the tight correlation between the two gradients in the stellar halo. We find no statistical evidence of additional information content related to accretion histories in stellar halo metallicity profiles beyond what is contained in surface-brightness profiles. Age gradients in the stellar halo do not appear to be sensitive to galactic accretion histories, and none of the stellar population gradients studied are strongly correlated with the mean merger mass-ratio. Future observations that reach large radii outside galaxies will have the best potential to constrain galactic accretion histories.

We test whether halo age and galaxy age are correlated at fixed halo and galaxy mass. The formation histories, and thus ages, of dark matter halos correlate with their large-scale density $\rho$, an effect known as assembly bias. We test whether this correlation extends to galaxies by measuring the dependence of galaxy stellar age on $\rho$. To clarify the comparison between theory and observation, and to remove the strong environmental effects on satellites, we use galaxy group catalogs to identify central galaxies and measure their quenched fraction, $f_Q$, as a function of large-scale environment. Models that match halo age to central galaxy age predict a strong positive correlation between $f_Q$ and $\rho$. However, we show that the amplitude of this effect depends on the definition of halo age: assembly bias is significantly reduced when removing the effects of splashback halos---those halos that are central but have passed through a larger halo or experienced strong tidal encounters. Defining age using halo mass at its peak value rather than current mass removes these effects. In SDSS data, at M$_{\rm gal}\gtrsim 10^{10.0}$ M_sol/h$^2$, there is a $\sim 5\%$ increase in $f_Q$ from low to high densities, which is in agreement with predictions of dark matter halos using peak halo mass. At lower stellar mass there is little to no correlation of $f_Q$ with $\rho$. For these galaxies, age-matching is inconsistent with the data across the wide range the halo formation metrics that we tested. This implies that halo formation history has a small but statistically significant impact on quenching of star formation at high masses, while the quenching process in low-mass central galaxies is uncorrelated with halo formation history.

Recently a population of large, very low surface brightness, spheroidal galaxies was identified in the Coma cluster. The apparent survival of these Ultra Diffuse Galaxies (UDGs) in a rich cluster suggests that they have very high masses. Here we present the stellar kinematics of Dragonfly 44, one of the largest Coma UDGs, using a 33.5 hr integration with DEIMOS on the Keck II telescope. We find a velocity dispersion of 47 km/s, which implies a dynamical mass of M_dyn=0.7x10^10 M_sun within its deprojected half-light radius of r_1/2=4.6 kpc. The mass-to-light ratio is M/L=48 M_sun/L_sun, and the dark matter fraction is 98 percent within the half-light radius. The high mass of Dragonfly 44 is accompanied by a large globular cluster population. From deep Gemini imaging taken in 0.4" seeing we infer that Dragonfly 44 has 94 globular clusters, similar to the counts for other galaxies in this mass range. Our results add to other recent evidence that many UDGs are "failed" galaxies, with the sizes, dark matter content, and globular cluster systems of much more luminous objects. We estimate the total dark halo mass of Dragonfly 44 by comparing the amount of dark matter within r=4.6 kpc to enclosed mass profiles of NFW halos. The enclosed mass suggests a total mass of ~10^12 M_sun, similar to the mass of the Milky Way. The existence of nearly-dark objects with this mass is unexpected, as galaxy formation is thought to be maximally-efficient in this regime.

Stellar spectra are often modeled and fit by interpolating within a rectilinear grid of synthetic spectra to derive the stars' labels: stellar parameters and elemental abundances. However, the number of synthetic spectra needed for a rectilinear grid grows exponentially with the label space dimensions, precluding the simultaneous and self-consistent fitting of more than a few elemental abundances. Shortcuts such as fitting subsets of labels separately can introduce unknown systematics and do not produce correct error covariances in the derived labels. In this paper we present a new approach -- CHAT (Convex Hull Adaptive Tessellation) -- which includes several new ideas for inexpensively generating a sufficient stellar synthetic library, using linear algebra and the concept of an adaptive, data-driven grid. A convex hull approximates the region where the data lie in the label space. A variety of tests with mock datasets demonstrate that CHAT can reduce the number of required synthetic model calculations by three orders of magnitude in an 8D label space. The reduction will be even larger for higher-dimensional label spaces. In CHAT the computational effort increases only linearly with the number of labels that are fit simultaneously. Around each of these grid points in label space an approximate synthetic spectrum can be generated through linear expansion using a set of "gradient spectra" that represent flux derivatives at every wavelength point with respect to all labels. These techniques provide new opportunities to fit the full stellar spectra from large surveys with 15-30 labels simultaneously.

Stellar spectra depend on the stellar parameters and on dozens of photospheric elemental abundances. Simultaneous fitting of these $\mathcal{N}\sim 10-40$ model labels to observed spectra has been deemed unfeasible, because the number of ab initio spectral model grid calculations scales exponentially with $\mathcal{N}$. We suggest instead the construction of a polynomial spectral model (PSM) of order $\mathcal{O}$ for the model flux at each wavelength. Building this approximation requires a minimum of only ${\mathcal{N}+\mathcal{O}\choose\mathcal{O}}$ calculations: e.g. a quadratic spectral model ($\mathcal{O}=2$) to fit $\mathcal{N}=20$ labels simultaneously, can be constructed from as few as $231$ ab initio spectral model calculations; in practice, a somewhat larger number ($\sim 300-1000$) of randomly chosen models lead to a better performing PSM. Such a PSM can be a good approximation only over a portion of label space, which will vary case by case. Yet, taking the APOGEE survey as an example, a single quadratic PSM provides a remarkably good approximation to the exact ab initio spectral models across much of this survey: for random labels within that survey the PSM approximates the flux to within $10^{-3}$, and recovers the abundances to within $\sim 0.02$ dex rms of the exact models. This enormous speed-up enables the simultaneous many-label fitting of spectra with computationally expensive ab initio models for stellar spectra, such as non-LTE models. A PSM also enables the simultaneous fitting of observational parameters, such as the spectrum's continuum or line-spread function.

This is the first of a series of papers presenting the Modules for Experiments in Stellar Astrophysics (MESA) Isochrones and Stellar Tracks (MIST) project, a new comprehensive set of stellar evolutionary tracks and isochrones computed using MESA, a state-of-the-art open-source 1D stellar evolution package. In this work, we present models with solar-scaled abundance ratios covering a wide range of ages ($5 \leq \rm \log(Age)\;[yr] \leq 10.3$), masses ($0.1 \leq M/M_{\odot} \leq 300$), and metallicities ($-2.0 \leq \rm [Z/H] \leq 0.5$). The models are self-consistently and continuously evolved from the pre-main sequence to the end of hydrogen burning, the white dwarf cooling sequence, or the end of carbon burning, depending on the initial mass. We also provide a grid of models evolved from the pre-main sequence to the end of core helium burning for $-4.0 \leq \rm [Z/H] < -2.0$. We showcase extensive comparisons with observational constraints as well as with some of the most widely used existing models in the literature. The evolutionary tracks and isochrones can be downloaded from the project website at http://waps.cfa.harvard.edu/MIST/.

We present the first direct comparison between Balmer line and panchromatic SED-based SFRs for z~2 galaxies. For this comparison we used 17 star-forming galaxies selected from the MOSFIRE Deep Evolution Field (MOSDEF) survey, with $3\sigma$ detections for H$\alpha$ and at least two IR bands (Spitzer/MIPS 24$\mu$m and Herschel/PACS 100 and 160$\mu$m, and in some cases Herschel/SPIRE 250, 350, and 500$\mu$m). The galaxies have total IR (8-1000$\mu$m) luminosities of $\sim10^{11.4}-10^{12.4}\,\textrm{L}_\odot$ and star-formation rates (SFRs) of $\sim30-250\,\textrm{M}_\odot\,\mathrm{yr^{-1}}$. We fit the UV-to-far-IR SEDs with flexible stellar population synthesis (FSPS) models - which include both stellar and dust emission - and compare the inferred SFRs with the SFR(H$\alpha$,H$\beta$) values corrected for dust attenuation using Balmer decrements. The two SFRs agree with a scatter of 0.17 dex. Our results imply that the Balmer decrement accurately predicts the obscuration of the nebular lines and can be used to robustly calculate SFRs for star-forming galaxies at z~2 with SFRs up to $\sim200\,\textrm{M}_\odot\,\mathrm{yr^{-1}}$. We also use our data to assess SFR indicators based on modeling the UV-to-mid-IR SEDs or by adding SFR(UV) and SFR(IR), for which the latter is based on the mid-IR only or on the full IR SED. All these SFRs show a poorer agreement with SFR(H$\alpha$,H$\beta$) and in some cases large systematic biases are observed. Finally, we show that the SFR and dust attenuation derived from the UV-to-near-IR SED alone are unbiased when assuming a delayed exponentially declining star-formation history.

The analysis of stellar populations has, by and large, been developed for two limiting cases: spatially-resolved stellar populations in the color-magnitude diagram, and integrated light observations of distant systems. In between these two extremes lies the semi-resolved regime, which encompasses a rich and relatively unexplored realm of observational phenomena. Here we develop the concept of pixel color magnitude diagrams (pCMDs) as a powerful technique for analyzing stellar populations in the semi-resolved regime. pCMDs show the distribution of imaging data in the plane of pixel luminosity vs. pixel color. A key feature of pCMDs is that they are sensitive to all stars, including both the evolved giants and the unevolved main sequence stars. An important variable in this regime is the mean number of stars per pixel, $N_{\rm pix}$. Simulated pCMDs demonstrate a strong sensitivity to the star formation history (SFH) and allow one to break degeneracies between age, metallicity and dust based on two filter data for values of $N_{\rm pix}$ up to at least $10^4$. We extract pCMDs from Hubble Space Telescope (HST) optical imaging of M31 and derive non-parametric SFHs from $10^6$ yr to $10^{10}$ yr for both the crowded disk and bulge regions (where $N_{\rm pix}\approx30-10^3$). From analyzing a small region of the disk we find a non-parametric SFH that is smooth and consistent with an exponential decay timescale of 4 Gyr. The bulge SFH is also smooth and consistent with a 2 Gyr decay timescale. pCMDs will likely play an important role in maximizing the science returns from next generation ground and space-based facilities.

Stars born from the same molecular cloud should be nearly homogeneous in their element abundances. The concept of chemical tagging is to identify members of disrupted clusters by their clustering in element abundance space. Chemical tagging requires large samples of stars with precise abundances for many individual elements. With uncertainties of $\sigma_{[X/{\rm Fe}]}$ and $\sigma_{\rm [Fe/H]} \simeq 0.05$ for 10 elements measured for $>10^4$ stars, the APOGEE DR12 spectra may be the first well-suited data set to put this idea into practice. We find that even APOGEE data offer only $\sim 500$ independent volume elements in the 10-dimensional abundance space, when we focus on the $\alpha$-enhanced Galactic disk. We develop and apply a new algorithm to search for chemically homogeneous sets of stars against a dominant background. By injecting star clusters into the APOGEE data set we show that chemically homogeneous clusters with masses $\gtrsim 3 \times 10^7 \, {\rm M}_\odot$ would be easily detectable and yet no such signal is seen in the data. By generalizing this approach, we put a first abundance-based constraint on the cluster mass function for the old disk stars in the Milky Way.

Long period variable stars arise in the final stages of the asymptotic giant branch phase of stellar evolution. They have periods of up to ~1000d and amplitudes that can exceed a factor of three in the I-band flux. These stars pulsate predominantly in their fundamental mode, which is a function of mass and radius, and so the pulsation periods are sensitive to the age of the underlying stellar population. The overall number of long period variables in a population is directly related to their lifetime, which is difficult to predict from first principles because of uncertainties associated with stellar mass-loss and convective mixing. The time variability of these stars has not been previously taken into account when modeling the spectral energy distributions of galaxies. Here we construct time-dependent stellar population models that include the effects of long period variable stars, and report the ubiquitous detection of this expected `pixel shimmer' in the massive metal-rich galaxy M87. The pixel light curves display a variety of behaviors, including linearly rising and falling curves, semi-periodic curves, and sudden increases or decreases in the flux level. The observed variation of 0.1-1% is very well matched to the predictions of our models. The data provide a strong and novel constraint on the properties of variable stars in an old and metal-rich stellar population, and we infer that the lifetime of long period variables in M87 is shorter by approximately 30% compared to predictions from the latest stellar evolution models.

The influence of both major and minor mergers is expected to significantly affect gradients of stellar ages and metallicities in the outskirts of galaxies. Measurements of observed gradients are beginning to reach large radii in galaxies, but a theoretical framework for connecting the findings to a picture of galactic build-up is still in its infancy. We analyze stellar populations of a statistically representative sample of quiescent galaxies over a wide mass range from the Illustris simulation. We measure metallicity and age profiles in the stellar halos of quiescent Illustris galaxies ranging in stellar mass from $10^{10}$ to $10^{12} M_\odot$, accounting for observational projection and luminosity-weighting effects. We find wide variance in stellar population gradients between galaxies of similar mass, with typical gradients agreeing with observed galaxies. We show that, at fixed mass, the fraction of stars born in-situ within galaxies is correlated with the metallicity gradient in the halo, confirming that stellar halos contain unique information about the build-up and merger histories of galaxies.

We present J-band spectroscopy of passive galaxies focusing on the Na I doublet at 1.14 {\mu}m. Like the Na I 0.82 {\mu}m doublet, this feature is strong in low-mass stars and hence may provide a useful probe of the initial mass function (IMF). From high signal-to-noise composite spectra, we find that Na I 1.14 {\mu}m increases steeply with increasing velocity dispersion, {\sigma}, and for the most massive galaxies (\sigma > 300 km/s) is much stronger than predicted from synthetic spectra with Milky-Way-like IMFs and solar abundances. Reproducing Na I 1.14 {\mu}m at high {\sigma} likely requires either a very high [Na/H], or a bottom-heavy IMF, or a combination of both. Using the Na D line to break the degeneracy between IMF and abundance, we infer [Na/H] $\approx$ +0.5 and a steep IMF (single-slope-equivalent x $\approx$ 3.2, where x = 2.35 for Salpeter), for the high-\sigma galaxies. At lower mass ({\sigma} = 50-100 km/s), the line strengths are compatible with MW-like IMFs and near-solar [Na/H]. We highlight two galaxies in our sample where strong gravitational lensing masses favour MW-like IMFs. Like the high-{\sigma} sample on average, these galaxies have strong Na I 1.14 \mu m; taken in isolation their sodium indices imply bottom-heavy IMFs which are hard to reconcile with the lensing masses. An alternative full-spectrum-fitting approach, applied to the high-\sigma sample, recovers an IMF less heavy than Salpeter, but under-predicts the Na I 1.14 \mu m line at the 5{\sigma} level. We conclude that current models struggle to reproduce this feature in the most massive galaxies without breaking other constraints, and caution against over-reliance on the sodium lines in spectroscopic IMF studies.

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

For the first time in history, humans have reached the point where it is possible to construct a revolutionary space-based observatory that has the capability to find dozens of Earth-like worlds, and possibly some with signs of life. This same telescope, designed as a long-lived facility, would also produce transformational scientific advances in every area of astronomy and astrophysics from black hole physics to galaxy formation, from star and planet formation to the origins of the Solar System. The Association of Universities for Research in Astronomy (AURA) commissioned a study on a next-generation UVOIR space observatory with the highest possible scientific impact in the era following JWST. This community-based study focuses on the future space-based options for UV and optical astronomy that significantly advance our understanding of the origin and evolution of the cosmos and the life within it. The committee concludes that a space telescope equipped with a 12-meter class primary mirror can find and characterize dozens of Earth-like planets and make fundamental advances across nearly all fields of astrophysics. The concept is called the High Definition Space Telescope (HDST). The telescope would be located at the Sun-Earth L2 point and would cover a spectral range that, at a minimum, runs from 0.1 to 2 microns. Unlike JWST, HDST will not need to operate at cryogenic temperatures. HDST can be made to be serviceable on orbit but does not require servicing to complete its primary scientific objectives. We present the scientific and technical requirements for HDST and show that it could allow us to determine whether or not life is common outside the Solar System. We do not propose a specific design for such a telescope, but show that designing, building and funding such a facility is feasible beginning in the next decade - if the necessary strategic investments in technology begin now.

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.

It is now well-established that the elemental abundance patterns of stars holds key clues not only to their formation but also to the assembly histories of galaxies. One of the most exciting possibilities is the use of stellar abundance patterns as "chemical tags" to identify stars that were born in the same molecular cloud. In this paper we assess the prospects of chemical tagging as a function of several key underlying parameters. We show that in the fiducial case of $10^4$ distinct cells in chemical space and $10^5-10^6$ stars in the survey, one can expect to detect $\sim 10^2-10^3$ groups that are $\ge 5\sigma$ overdensities in the chemical space. However, we find that even very large overdensities in chemical space do not guarantee that the overdensity is due to a single set of stars from a common birth cloud. In fact, for our fiducial model parameters, the typical $5\sigma$ overdensity is comprised of stars from a wide range of clusters with the most dominant cluster contributing only 25% of the stars. The most important factors limiting the identification of disrupted clusters via chemical tagging are the number of chemical cells in the chemical space and the survey sampling rate of the underlying stellar population. Both of these factors can be improved through strategic observational plans. While recovering individual clusters through chemical tagging may prove challenging, we show, in agreement with previous work, that different CMFs imprint different degrees of clumpiness in chemical space. These differences provide the opportunity to statistically reconstruct the slope and high mass cutoff of CMF and its evolution through cosmic time.

Selecting centrally quiescent galaxies from the Sloan Digital Sky Survey (SDSS) to create high signal-to-noise (>100) stacked spectra with minimal emission line contamination, we accurately and precisely model the central stellar populations of barred and unbarred quiescent disk galaxies. By splitting our sample by redshift, we can use the fixed size of the SDSS fiber to model the stellar populations at different radii within galaxies. At 0.02<z<0.04, the SDSS fiber radius corresponds to ~1 kpc, which is the typical half-light radii of both classical bulges and disky pseudobulges. Assuming that the SDSS fiber primarily covers the bulges at these redshifts, our analysis shows that there are no significant differences in the stellar populations, i.e., stellar age, [Fe/H], [Mg/Fe], and [N/Fe], of the bulges of barred vs. unbarred quiescent disk galaxies. Modeling the stellar populations at different redshift intervals from z=0.020 to z=0.085 at fixed stellar masses produces an estimate of the stellar population gradients out to about half the typical effective radius of our sample, assuming null evolution over this ~1 Gyr epoch. We find that there are no noticeable differences in the slopes of the azimuthally averaged gradients of barred vs. unbarred quiescent disk galaxies. These results suggest that bars are not a strong influence on the chemical evolution of quiescent disk galaxies.

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

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

We 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 results from a blind survey to identify strong gravitational lenses among the population of low-redshift early-type galaxies. The SINFONI Nearby Elliptical Lens Locator Survey (SNELLS) uses integral-field infrared spectroscopy to search for lensed emission line sources behind massive lens candidates at $z$<0.055. From 27 galaxies observed, we have recovered one previously-known lens (ESO325-G004) at $z$=0.034, and discovered two new systems, at $z$=0.031 and $z$=0.052. All three lens galaxies have high velocity dispersions (\sigma>300 km/s) and \alpha-element abundances ([Mg/Fe]>0.3). From the lensing configurations we derive total J-band mass-to-light ratios of 1.8$\pm$0.1, 2.1$\pm$0.1 and 1.9$\pm$0.2 within the $\sim$2 kpc Einstein radius. Correcting for estimated dark-matter contributions, and comparing to stellar population models with a Milky Way (Kroupa) initial mass function (IMF), we determine the "mass excess factor", \alpha. Assuming the lens galaxies have "old" stellar populations (10$\pm$1 Gyr), the average IMF mass factor is $\langle\alpha\rangle$=1.10$\pm$0.08$\pm$0.10, where the first error is random and the second is systematic. If we instead fit the stellar populations from 6dF optical survey spectra, all three galaxies are consistent with being old, but the age errors are 3-4 Gyr, due to limited signal-to-noise. The IMF constraints are therefore looser in this case, with $\langle\alpha\rangle$ = $1.23^{+0.16}_{-0.13}\pm{0.10}$. Our results are thus consistent with a Kroupa IMF (\alpha=1.00) on average, and strongly reject very heavy IMFs with \alpha>2. A Salpeter IMF (\alpha=1.55) is inconsistent at the 3.5$\sigma$ level if the galaxies are old, but cannot be excluded using age constraints derived from the currently-available optical spectra.

We revisit previous suggestions that the heating provided by the winds of dying low-mass stars plays an important role in preventing star formation in quiescent galaxies. At the end of their asymptotic giant branch phase, intermediate and low-mass stars eject their envelopes rapidly in a super-wind phase, usually giving rise to planetary nebulae. In spheroidal galaxies with high stellar velocity dispersions, the interaction of these ejected envelopes with the ambient diffuse gas can lead to significant, isotropic and steady-state heating that scales as $\dot{M}_\ast\sigma_\ast^2$. We show that cooling of the central regions of the hot diffuse halo gas can be delayed for a Hubble time for halos more massive than $\sim10^{12.5}M_{\odot}$ at $0<z<2$, although stellar heating alone is unlikely to forestall cooling in the most massive clusters at $z=0$. This mechanism provides a natural explanation for the strong trend of galaxy quiescence with stellar surface density and velocity dispersion. In addition, since the ejected material will thermalize to $kT\sim\sigma_\ast^2$, this mechanism provides an explanation for the observed similarity between the central temperature of the hot diffuse gas and $\sigma_\ast^2$, a result which is not trivial in light of the short inferred cooling times of the hot gas. The main uncertainty in this analysis is the ultimate fate of the stellar ejecta. Preventing accumulation of the ejecta in the central regions may require energy input from another source, such as Type Ia supernovae. Detailed simulations of the interaction of the stellar wind with the ambient gas are required to better quantify the net effect of AGB heating.

We report the discovery of 47 low surface brightness objects in deep images of a 3 x 3 degree field centered on the Coma cluster, obtained with the Dragonfly Telephoto Array. The objects have central surface brightness mu(g,0) ranging from 24 - 26 mag/arcsec^2 and effective radii r_e = 3"-10", as measured from archival Canada France Hawaii Telescope images. From their spatial distribution we infer that most or all of the objects are galaxies in the Coma cluster. This relatively large distance is surprising as it implies that the galaxies are very large: with r_e = 1.5 - 4.6 kpc their sizes are similar to those of L* galaxies even though their median stellar mass is only ~6 x 10^7 Solar masses. The galaxies are relatively red and round, with <g-i> = 0.8 and <b/a> = 0.74. One of the 47 galaxies is fortuitously covered by a deep Hubble Space Telescope ACS observation. The ACS imaging shows a large spheroidal object with a central surface brightness mu(g,0) = 25.8 mag/arcsec^2, a Sersic index n=0.6, and an effective radius of 7", corresponding to 3.4 kpc at the distance of Coma. The galaxy is not resolved into stars, consistent with expectations for a Coma cluster object. We speculate that UDGs may have lost their gas supply at early times, possibly resulting in very high dark matter fractions.

We present an overview of a new integral field spectroscopic survey called MaNGA (Mapping Nearby Galaxies at Apache Point Observatory), one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) that began on 2014 July 1. MaNGA will investigate the internal kinematic structure and composition of gas and stars in an unprecedented sample of 10,000 nearby galaxies. We summarize essential characteristics of the instrument and survey design in the context of MaNGA's key science goals and present prototype observations to demonstrate MaNGA's scientific potential. MaNGA employs dithered observations with 17 fiber-bundle integral field units that vary in diameter from 12" (19 fibers) to 32" (127 fibers). Two dual-channel spectrographs provide simultaneous wavelength coverage over 3600-10300 A at R~2000. With a typical integration time of 3 hr, MaNGA reaches a target r-band signal-to-noise ratio of 4-8 (per A, per 2" fiber) at 23 AB mag per sq. arcsec, which is typical for the outskirts of MaNGA galaxies. Targets are selected with stellar mass greater than 1e9 Msun using SDSS-I redshifts and i-band luminosity to achieve uniform radial coverage in terms of the effective radius, an approximately flat distribution in stellar mass, and a sample spanning a wide range of environments. Analysis of our prototype observations demonstrates MaNGA's ability to probe gas ionization, shed light on recent star formation and quenching, enable dynamical modeling, decompose constituent components, and map the composition of stellar populations. MaNGA's spatially resolved spectra will enable an unprecedented study of the astrophysics of nearby galaxies in the coming 6 yr.

We present stellar atmosphere and evolution models of main sequence stars in two stellar populations of the Galactic globular cluster NGC 6752. These populations represent the two extremes of light-element abundance variations in the cluster. NGC 6752 is a benchmark cluster in the study of multiple stellar populations because of the rich array of spectroscopic abundances and panchromatic Hubble Space Telescope photometry. The spectroscopic abundances are used to compute stellar atmosphere and evolution models. The synthetic spectra for the two populations show significant differences in the ultraviolet and, for the coolest temperatures, in the near-infrared. The stellar evolution models exhibit insignificant differences in the H-R diagram except on the lower main sequence. The appearance of multiple sequences in the colour-magnitude diagrams (CMDs) of NGC 6752 is almost exclusively due to spectral effects caused by the abundance variations. The models reproduce the observed splitting and/or broadening of sequences in a range of CMDs. The ultraviolet CMDs are sensitive to variations in carbon, nitrogen, and oxygen but the models are not reliable enough to directly estimate abundance variations from photometry. On the other hand, the widening of the lower main sequence in the near-infrared CMD, driven by oxygen-variation via the water molecule, is well-described by the models and can be used to estimate the range of oxygen present in a cluster from photometry. We confirm that it is possible to use multiband photometry to estimate helium variations among the different populations, with the caveat that the estimated amount of helium-enhancement is model-dependent.

We introduce a new method to determine the relative contributions of different types of stars to the integrated light of nearby early-type galaxies. As is well known, the surface brightness of these galaxies shows pixel-to-pixel fluctuations due to Poisson variations in the number of giant stars. Differential spectroscopy of pixels as a function of fluctuation strength ("fluctuation spectroscopy") effectively measures the spectral variation of stars as a function of their luminosity, information that is otherwise difficult to obtain for individual stars outside of the Local Group. We apply this technique to the elliptical galaxy NGC 4472, using HST/ACS imaging in six narrow-band ramp filters tuned to spectral features in the range 0.8-1.0 micron. Pixels with +- 5% broad-band variations show differential color variations of 0.1% - 1.0% in the narrow-band filters. These variations are primarily due to the systematic increase in TiO absorption strength with increasing luminosity on the upper giant branch. The data are very well reproduced by the same Conroy & van Dokkum (2012) stellar population synthesis model that is the best fit to the integrated light, with residuals in the range 0.03% - 0.09%. Models with ages or metallicities that are significantly different from the integrated-light values do not yield good fits. We can also rule out several modifications to the underlying model, including the presence of a significant (>3% of the light) population of late M giants. The current observations constitute a powerful test of the expected luminosities and temperatures of metal-rich giants in massive early-type galaxies. Studies of pixels with much larger (negative) fluctuations will provide unique information on main sequence stars and the stellar initial mass function.

As part of the SLUGGS survey, we stack 1137 Keck DEIMOS spectra of globular clusters from 10 galaxies to study their stellar populations in detail. The stacked spectra have median signal to noise ratios of $\sim 90$ \AA$^{-1}$. Besides the calcium triplet, we study weaker sodium, magnesium, titanium and iron lines as well as the H$\alpha$ and higher order Paschen hydrogen lines. In general, the stacked spectra are consistent with old ages and a Milky Way-like initial mass function. However, we see different metal line index strengths at fixed colour and magnitude, and differences in the calcium triplet--colour relation from galaxy to galaxy. We interpret this as strong evidence for variations in the globular cluster colour--metallicity relation between galaxies. Two possible explanations for the colour--metallicity relation variations are that the average ages of globular clusters vary from galaxy to galaxy or that the average abundances of light elements (i.e. He, C, N and O) differ between galaxies. Stacking spectra by magnitude, we see that the colours become redder and metal line indices stronger with brighter magnitudes. These trends are consistent with the previously reported `blue tilts' being mass--metallicity relations.

We present a new technique to estimate the evolution of the very faint end of the UV luminosity function (LF) out to $z\sim5$. Measured star formation histories (SFHs) from the fossil record of Local Group galaxies are used to reconstruct the LF down to M$_{UV}\sim-5$ at $z\sim5$ and M$_{UV}\sim-1.5$ at $z<1$. Such faint limits are well beyond the current observational limits and are likely to remain beyond the limits of next generation facilities. The reconstructed LFs, when combined with direct measurements of the LFs at higher luminosity, are well-fit by a standard Schechter function with no evidence of a break to the faintest limits probed by this technique. The derived faint end slope, $\alpha$, steepens from $\approx-1.2$ at $z<1$ to $\approx-1.6$ at $4<z<5$. We test the effects of burstiness in the SFHs and find the recovered LFs to be only modestly affected. Incompleteness corrections for the faintest Local Group galaxies and the (unlikely) possibility of significant luminosity-dependent destruction of dwarf galaxies between high redshift and the present epoch are important uncertainties. These and other uncertainties can be mitigated with more detailed modeling and future observations. The reconstructed faint end LF from the fossil record can therefore be a powerful and complementary probe of the high redshift faint galaxies believed to play a key role in the reionization of the Universe.

This paper explores several simple model variations for the connections among quasars, galaxies, and dark matter halos for redshifts 1 < z < 6. A key component of these models is that we enforce a self-consistent black hole (BH) history by tracking both BH mass and BH growth rate at all redshifts. We connect objects across redshift with a simple constant-number-density procedure, and choose a fiducial model with a relationship between BH and galaxy growth rates that is linear and evolves in a simple way with redshift. Within this fiducial model, we find the quasar luminosity function (QLF) by calculating an "intrinsic" luminosity based on either the BH mass or BH growth rate, and then choosing a model of quasar variability with either a lognormal or truncated power-law distribution of instantaneous luminosities. This gives four model variations, which we fit to the observed QLF at each redshift. With the best-fit models in hand, we undertake a detailed comparison of the four fiducial models, and explore changes to our fiducial model of the BH-galaxy relationship. Each model variation can successfully fit the observed QLF, the shape of which is generally set by the "intrinsic" luminosity at the faint end and by the scatter due to variability at the bright end. We focus on accounting for the reasons that physically different models can make such similar predictions, and on identifying what observational data or physical arguments are most essential in breaking the degeneracies among models.

We present results from modeling the optical spectra of a large sample of quiescent galaxies between 0.1 < z < 0.7 from the Sloan Digital Sky Survey (SDSS) and the AGN and Galaxy Evolution Survey (AGES). We examine how the stellar ages and abundance patterns of galaxies evolve over time as a function of stellar mass from 10^{9.6}-10^{11.8} Msun. Galaxy spectra are stacked in bins of mass and redshift, and modeled over a wavelength range from 4000 A to 5500 A. Full spectrum stellar population synthesis modeling provides estimates of the age and the abundances of the elements Fe, Mg, C, N, and Ca. We find negligible evolution in elemental abundances at fixed stellar mass over roughly 7 Gyr of cosmic time. In addition, the increase in stellar ages with time for massive galaxies is consistent with passive evolution since z = 0.7. Taken together, these results favor a scenario in which the inner ~ 0.3-3 R_e of massive quiescent galaxies have been passively evolving over the last half of cosmic time. Interestingly, the derived stellar ages are considerably younger than the age of the universe at all epochs, consistent with an equivalent single-burst star formation epoch of z < 1.5. These young stellar population ages coupled with the existence of massive quiescent galaxies at z > 1 indicate the inhomogeneous nature of the z < 0.7 quiescent population. The data also permit the addition of newly-quenched galaxies at masses below ~10^{10.5} Msun at z < 0.7. Additionally, we analyze very deep Keck DEIMOS spectra of the two brightest quiescent galaxies in a cluster at z = 0.83. There is tentative evidence that these galaxies are older than their counterparts in low-density environments. In the Appendix, we demonstrate that our full spectrum modeling technique allows for accurate and reliable modeling of galaxy spectra to low S/N (~20 A^{-1}) and/or low spectral resolution (R ~ 500).

The dense interiors of massive galaxies are among the most intriguing environments in the Universe. In this paper we ask when these dense cores were formed and determine how galaxies gradually assembled around them. We select galaxies that have a stellar mass >3x10^10 Msun inside r=1 kpc out to z=2.5, using the 3D-HST survey and data at low redshift. Remarkably, the number density of galaxies with dense cores appears to have decreased from z=2.5 to the present. This decrease is probably mostly due to stellar mass loss and the resulting adiabatic expansion, with some contribution from merging. We infer that dense cores were mostly formed at z>2.5, consistent with their largely quiescent stellar populations. While the cores appear to form early, the galaxies in which they reside show strong evolution: their total masses increase by a factor of 2-3 from z=2.5 to z=0 and their effective radii increase by a factor of 5-6. As a result, the contribution of dense cores to the total mass of the galaxies in which they reside decreases from ~50% at z=2.5 to ~15% at z=0. Because of their early formation, the contribution of dense cores to the total stellar mass budget of the Universe is a strong function of redshift. The stars in cores with M_1kpc>3x10^10 Msun make up ~0.1% of the stellar mass density of the Universe today but 10%-20% at z~2, depending on their IMF. The formation of these cores required the conversion of ~10^11 Msun of gas into stars within ~1 kpc, while preventing significant star formation at larger radii.

On dwarf galaxy scales, the different shapes of the galaxy stellar mass function and the dark halo mass function require a star-formation efficiency (SFE) in these systems that is currently more than 1 dex lower than that of Milky Way-size halos. Here, we argue that this trend may actually be reversed at high redshift. Specifically, by combining the resolved star-formation histories of nearby isolated dwarfs with the simulated mass-growth rates of dark matter halos, we show that the assembly of these systems occurs in two phases: (1) an early, fast halo accretion phase with a rapidly deepening potential well, characterized by a high SFE; and (2) a late slow halo accretion phase where, perhaps as a consequence of reionization, the SFE is low. Nearby dwarfs have more old stars than predicted by assuming a constant or decreasing SFE with redshift, a behavior that appears to deviate qualitatively from the trends seen amongst more massive systems. Taken at face value, the data suggest that, at sufficiently early epochs, dwarf galaxy halos above the atomic cooling mass limit can be among the most efficient sites of star formation in the universe.

We present results from a fully cosmological, very high-resolution, LCDM "zoom-in" simulation of a group of seven field dwarf galaxies with present-day virial masses in the range M_vir=4.4e8-3.6e10 Msun. The simulation includes a blastwave scheme for supernova feedback, a star formation recipe based on a high gas density threshold, metal-dependent radiative cooling, a scheme for the turbulent diffusion of metals and thermal energy, and a uniform UV background that modifies the ionization and excitation state of the gas. The properties of the simulated dwarfs are strongly modulated by the depth of the gravitational potential well. All three halos with M_vir < 1e9 Msun are devoid of stars, as they never reach the density threshold for star formation of 100 atoms/cc. The other four, M_vir > 1e9 Msun dwarfs have blue colors, low star formation efficiencies, high cold gas to stellar mass ratios, and low stellar metallicities. Their bursty star formation histories are characterized by peak specific star formation rates in excess of 50-100 1/Gyr, far outside the realm of normal, more massive galaxies, and in agreement with observations of extreme emission-line starbursting dwarfs by the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey. Metal-enriched galactic outflows produce sub-solar effective yields and pollute with heavy elements a Mpc-size region of the intergalactic medium, but are not sufficient to completely quench star formation activity and are not ubiquitous in our dwarfs. Within the limited size of the sample, our simulations appear to simultaneously reproduce the observed stellar mass and cold gas content, resolved star formation histories, stellar kinematics, and metallicities of field dwarfs in the Local Volume.

The nature of dust in distant galaxies is not well understood, and until recently few direct dust measurements have been possible. We investigate dust in distant star-forming galaxies using near-infrared grism spectra of the 3D-HST survey combined with archival multi-wavelength photometry. These data allow us to make a direct comparison between dust around star-forming regions ($A_{V,\mathrm{HII}}$) and the integrated dust content ($A_{V,\mathrm{star}}$). We select a sample of 163 galaxies between $1.36\le{}z\le1.5$ with H$\alpha$ signal-to-noise ratio $\ge5$ and measure Balmer decrements from stacked spectra to calculate $A_{V,\mathrm{HII}}$. First, we stack spectra in bins of $A_{V,\mathrm{star}}$, and find that $A_{V,\mathrm{HII}}=1.86\,A_{V,\mathrm{star}}$, with a significance of $\sigma=1.7$. Our result is consistent with the two-component dust model, in which galaxies contain both diffuse and stellar birth cloud dust. Next, we stack spectra in bins of specific star formation rate ($\log\,\mathrm{SSFR}$), star formation rate ($\log\,\mathrm{SFR}$), and stellar mass ($\log{}M_*$). We find that on average $A_{V,\mathrm{HII}}$ increases with SFR and mass, but decreases with increasing SSFR. Interestingly, the data hint that the amount of extra attenuation decreases with increasing SSFR. This trend is expected from the two-component model, as the extra attenuation will increase once older stars outside the star-forming regions become more dominant in the galaxy spectrum. Finally, using Balmer decrements we derive dust-corrected H$\alpha$ SFRs, and find that stellar population modeling produces incorrect SFRs if rapidly declining star formation histories are included in the explored parameter space.

We examine the star formation histories (SFHs) of galaxies in smoothed particle hydrodynamics (SPH) simulations, compare them to parametric models that are commonly used in fitting observed galaxy spectral energy distributions, and examine the efficacy of these parametric models as practical tools for recovering the physical parameters of galaxies. The commonly used tau-model, with SFR ~ exp(-t/tau), provides a poor match to the SFH of our SPH galaxies, with a mismatch between early and late star formation that leads to systematic errors in predicting colours and stellar mass-to-light ratios. A one-parameter lin-exp model, with SFR ~ t*exp(-t/tau), is much more successful on average, but it fails to match the late-time behavior of the bluest, most actively star-forming galaxies and the passive, "red and dead" galaxies. We introduce a 4-parameter model, which transitions from lin-exp to a linear ramp after a transition time, which describes our simulated galaxies very well. We test the ability of these parametrised models to recover (at z=0, 0.5, and 1) the stellar mass-to-light ratios, specific star formation rates, and stellar population ages from the galaxy colours, computed from the full SPH star formation histories using the FSPS code of Conroy et al. (2009). Fits with tau-models systematically overestimate M/L by ~ 0.2 dex, overestimate population ages by ~ 1-2 Gyr, and underestimate sSFR by ~ 0.05 dex. Fits with lin-exp are less biased on average, but the 4-parameter model yields the best results for the full range of galaxies. Marginalizing over the free parameters of the 4-parameter model leads to slightly larger statistical errors than 1-parameter fits but essentially removes all systematic biases, so this is our recommended procedure for fitting real galaxies.

Galaxies that are several virial radii beyond groups/clusters show preferentially quiescent star formation rates. Using a galaxy group/cluster catalog from the Sloan Digital Sky Survey, together with a cosmological N-body simulation, we examine the origin of this environmental quenching beyond the virial radius. Accounting for the clustering of groups/clusters, we show that central galaxies show enhanced SFR quenching out to 2.5 virial radii beyond groups/clusters, and we demonstrate that this extended environmental enhancement can be explained simply by 'ejected' satellite galaxies that orbit beyond their host halo's virial radius. We show that ejected satellites typically orbit for several Gyr beyond the virial radius before falling back in, and thus they compose up to 40% of all central galaxies near groups/clusters. We show that a model in which ejected satellites experience the same SFR quenching as satellites within a host halo can explain essentially all environmental dependence of galaxy quenching. Furthermore, ejected satellites (continue to) lose significant halo mass, an effect that is potentially observable via gravitational lensing. The SFRs/colors and stellar-to-halo masses of ejected satellites highlight the importance of environmental history and present challenges to models of galaxy occupation that ignore such history.

We have used near-ultraviolet (NUV) to mid-infrared (MIR) composite spectral energy distributions (SEDs) to simultaneously model the attenuated stellar and dust emission of 0.5 < z < 2.0 galaxies. These composite SEDs were previously constructed from the photometric catalogs of the NEWFIRM Medium-Band Survey, by stacking the observed photometry of galaxies that have similar rest-frame NUV-to-NIR SEDs. In this work, we include a stacked MIPS 24 micron measurement for each SED type to extend the SEDs to rest-frame MIR wavelengths. Consistent with previous studies, the observed MIR emission for most SED types is higher than expected from only the attenuated stellar emission. We fit the NUV-to-MIR composite SEDs by the Flexible Stellar Population Synthesis (SPS) models, which include both stellar and dust emission. We compare the best-fit star formation rates (SFRs) to the SFRs based on simple UV+IR estimators. Interestingly, the UV and IR luminosities overestimate SFRs - compared to the model SFRs - by more than ~ 1 dex for quiescent galaxies, while for the highest star-forming galaxies in our sample the two SFRs are broadly consistent. The difference in specific SFRs also shows a gradually increasing trend with declining specific SFR, implying that quiescent galaxies have even lower specific SFRs than previously found. Contributions from evolved stellar populations to both the UV and the MIR SEDs most likely explain the discrepancy. Based on this work, we conclude that SFRs should be determined from modeling the attenuated stellar and dust emission simultaneously, instead of employing simple UV+IR-based SFR estimators.

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

We investigate the optical spectral region of spectra of 1000 stars searching for IMF-sensitive features to constrain the low-mass end of the initial mass function (IMF) slope in elliptical galaxies. We use the MILES stellar library in the wavelength range [3500-7500] {\AA} to select indices that are sensitive to cool dwarf stars and that only weakly depend on age and metallicity. In this wavelength range, the response of a change in the effective temperature of the cool red giant population is similar to the response of a change in the number of dwarf stars in the galaxy. We investigate the degeneracy between IMF variation and Teff,RGB and show that it is possible to break it with the new IMF indicators defined here. In particular, we define a CaH1 index around {\lambda}6380 {\AA} that arises purely from cool dwarfs, does not strongly depend on age and is anti-correlated with [a/Fe]. This index allows the determination of the the IMF slope when combined with TiO lines and age- and metallicity-dependent features. The use of several indicators is crucial to break degeneracies between IMF, age, abundance pattern and Teff of the RGB population. We measure line-index strengths of the optical IMF indicators in the SSP models and compare these with index strengths of the same features in stacked SSDSS early-type galaxy spectra with varying velocity dispersions. Using 10 different indicators, we find a clear trend of a steepening IMF with increasing velocity dispersion described by the equation x=(2.3+-0.1)*log{\sigma}200+(2.13+-0.15), where x is the IMF slope and {\sigma}200 is the central stellar velocity dispersion measured in units of 200km/s. We compute the 'IMF mismatch' parameter as the ratio of stellar mass-to-light ratio predicted from the x-{\sigma}200 relation to that inferred from SSP models assuming a Salpeter IMF and find good agreement with independent published results.

The stellar populations of galaxies hold vital clues to their formation histories. In this paper we present results based on modeling stacked spectra of early-type galaxies drawn from the Sloan Digital Sky Survey (SDSS) as a function of velocity dispersion, sigma, from 90 km/s to 300 km/s. The spectra are of extremely high quality, with typical S/N of 1000/A, and a wavelength coverage of 4000A-8800A. Our population synthesis model includes variation in 16 elements from C to Ba, the shift in effective temperature, Delta(Teff), of the stars with respect to a solar metallicity isochrone, amongst other parameters. In our approach we fit the full optical spectra rather than a select number of spectral indices and are able to, for the first time, measure the abundances of the elements V, Cr, Mn, Co, and Ni from the integrated light of distant galaxies. Our main results are as follows: 1) light-weighted stellar ages range from 6-12 Gyr from low to high sigma; 2) [Fe/H] varies by less than 0.1 dex across the entire sample; 3) Mg closely tracks O, and both increase from ~0.0 at low sigma to ~0.25 at high sigma; Si and Ti show a shallower rise with sigma, and Ca tracks Fe rather than O; 4) the iron peak elements V, Cr, Mn, and Ni track Fe, while Co tracks O, suggesting that Co forms primarily in massive stars; 5) C and N track O over the full sample and [C/Fe] and [N/Fe] exceed 0.2 at high sigma; and 6) the variation in Delta(Teff) with total metallicity follows theoretical predictions based on stellar evolution theory. Our derived [Mg/Fe] and [O/Fe] abundance ratios are 0.05-0.1 dex lower than most previous determinations. Under the conventional interpretation that the variation in these ratios is due to star formation timescale variations, our results suggest longer star formation timescales for massive early-type galaxies than previous studies. (ABRIDGED)

We present Keck/MOSFIRE K-band spectroscopy of the first mass-selected sample of galaxies at $z\sim2.3$. Targets are selected from the 3D-HST Treasury survey. The six detected galaxies have a mean [NII]$\lambda$6584/H$\alpha$ ratio of $0.27\pm0.01$, with a small standard deviation of 0.05. This mean value is similar to that of UV-selected galaxies of the same mass. The mean gas-phase oxygen abundance inferred from the [NII]/H$\alpha$ ratios depends on the calibration method, and ranges from 12+log(O/H)$_{gas}=8.57$ for the {Pettini} & {Pagel} (2004) calibration to 12+log(O/H)$_{gas}= 8.87$ for the {Maiolino} {et~al.} (2008) calibration. Measurements of the stellar oxygen abundance in nearby quiescent galaxies with the same number density indicate 12+log(O/H)$_{stars}= 8.95$, similar to the gas-phase abundances of the $z\sim2.3$ galaxies if the {Maiolino} {et~al.} (2008) calibration is used. This suggests that these high-redshift star forming galaxies may be progenitors of today's massive early-type galaxies. The main uncertainties are the absolute calibration of the gas-phase oxygen abundance and the incompleteness of the $z\sim2.3$ sample: the galaxies with detected H$\alpha$ tend to be larger and have higher star formation rates than the galaxies without detected H$\alpha$, and we may still be missing the most dust-obscured progenitors.

Several independent lines of evidence suggest that the stellar initial mass function (IMF) in early-type galaxies becomes increasingly `bottom-heavy' with increasing galaxy mass and/or velocity dispersion, sigma. Here we consider evidence for IMF variation in a sample of relatively compact early-type galaxies drawn from the Sloan Digital Sky Survey. These galaxies are of sufficiently high stellar density that a dark halo likely makes a minor contribution to the total dynamical mass, Mdyn, within one effective radius. We fit our detailed stellar population synthesis models to the stacked absorption line spectra of these galaxies in bins of sigma and find evidence from IMF-sensitive spectral features for a bottom-heavy IMF at high sigma. We also apply simple `mass-follows-light' dynamical models to the same data and find that Mdyn is significantly higher than what would be expected if these galaxies were stellar dominated and had a universal Milky Way IMF. Adopting Mdyn M_* therefore implies that the IMF is `heavier' at high sigma. Most importantly, the quantitative amount of inferred IMF variation is very similar between the two techniques, agreeing to within 0.1 dex in mass. The agreement between two independent techniques, when applied to the same data, provides compelling evidence for systematic variation in the IMF as a function of early-type galaxy velocity dispersion. Any alternative explanations must reproduce both the results from dynamical and stellar population-based techniques.

We investigate the use of the halo mass-gap statistic --- defined as the logarithmic difference in mass between the host halo and its most massive satellite subhalo --- as a probe of halo age and concentration. A cosmological N-body simulation is used to study N ~25, 000 group/cluster sized halos in the mass range 10^12.5 < M_halo/M_sol < 10^14.5. In agreement with previous work, we find that halo mass-gap is related to halo formation time and concentration. On average, older and more highly concentrated halos have larger halo mass-gaps, and this trend is stronger than the mass-concentration relation over a similar dynamic range. However, there is a large amount of scatter owing to the transitory nature of the satellite subhalo population, which limits the use of the halo mass-gap statistic on an object-by-object basis. For example, we find that 20% of very large halo mass-gap systems (akin to "fossil groups") are young, and have likely experienced a recent merger between a massive satellite subhalo and the central subhalo. We relate halo mass-gap to the observable stellar mass-gap via abundance matching. Using a galaxy group catalog constructed from the Sloan Digital Sky Survey Data Release 7, we find that the star formation and structural properties of galaxies at fixed mass show no trend with stellar mass-gap. This is despite a variation in halo-age of ~2.5 Gyr over ~1.2 dex in stellar mass-gap. Thus, we find no evidence to suggest that the halo formation history significantly affects galaxy properties.

In the past decade, the notion that globular clusters (GCs) are composed of coeval stars with homogeneous initial chemical compositions has been challenged by growing evidence that they host an intricate stellar population mix, likely indicative of a complex history of star formation and chemical enrichment. Several models have been proposed to explain the existence of multiple stellar populations in GCs, but no single model provides a fully satisfactory match to existing data. Correlations between chemistry and global parameters such as cluster mass or luminosity are fundamental clues to the physics of GC formation. In this Letter, we present an analysis of the mean abundances of Fe, Mg, C, N, and Ca for 72 old GCs from the Andromeda galaxy. We show for the first time that there is a correlation between the masses of GCs and the mean stellar abundances of nitrogen, spanning almost two decades in mass. This result sheds new light on the formation of GCs, providing important constraints on their internal chemical evolution and mass loss history.

This letter utilizes composite spectral energy distributions (SEDs) constructed from NEWFIRM Medium-Band Survey photometry to constrain the dust attenuation curve in 0.5<z<2.0 galaxies. Based on similarities between the full SED shapes (0.3-8 micron), we have divided galaxies in 32 different spectral classes and stacked their photometry. As each class contains galaxies over a range in redshift, the resulting rest-frame SEDs are well-sampled in wavelength and show various spectral features including Halpha and the UV dust bump at 2175 Angstrom. We fit all composite SEDs with flexible stellar population synthesis models, while exploring attenuation curves with varying slopes and UV bump strengths. The Milky Way and Calzetti law provide poor fits at UV wavelengths for nearly all SEDs. Consistent with previous studies, we find that the best-fit attenuation law varies with spectral type. There is a strong correlation between the best-fit dust slope and UV bump strength, with steeper laws having stronger bumps. Moreover, the attenuation curve correlates with specific star formation rate (SFR), with more active galaxies having shallower dust curves and weaker bumps. There is also a weak correlation with inclination. The observed trends can be explained by differences in the dust-to-star geometry, a varying grain size distribution, or a combination of both. Our results have several implications for galaxy evolution studies. First, the assumption of a universal dust model leads to biases in derived galaxy properties. Second, the presence of a dust bump may result in underestimated values for the UV slope, used to correct SFRs of distant galaxies.

We report the discovery of a remarkable ultra-compact dwarf galaxy around the massive Virgo elliptical galaxy NGC 4649 (M60), which we term M60-UCD1. With a dynamical mass of 2.0 x 10^8 M_sun but a half-light radius of only ~ 24 pc, M60-UCD1 is more massive than any ultra-compact dwarfs of comparable size, and is arguably the densest galaxy known in the local universe. It has a two-component structure well-fit by a sum of Sersic functions, with an elliptical, compact (r_h=14 pc; n ~ 3.3) inner component and a round, exponential, extended (r_h=49 pc) outer component. Chandra data reveal a variable central X-ray source with L_X ~ 10^38 erg/s that could be an active galactic nucleus associated with a massive black hole or a low-mass X-ray binary. Analysis of optical spectroscopy shows the object to be old (~> 10 Gyr) and of solar metallicity, with elevated [Mg/Fe] and strongly enhanced [N/Fe] that indicates light element self-enrichment; such self-enrichment may be generically present in dense stellar systems. The velocity dispersion (~ 70 km/s) and resulting dynamical mass-to-light ratio (M/L_V=4.9 +/- 0.7) are consistent with---but slightly higher than---expectations for an old, metal-rich stellar population with a Kroupa initial mass function. The presence of a massive black hole or a mild increase in low-mass stars or stellar remnants is therefore also consistent with this M/L_V. The stellar density of the galaxy is so high that no dynamical signature of dark matter is expected. However, the properties of M60-UCD1 suggest an origin in the tidal stripping of a nucleated galaxy with M_B ~ -18 to -19.

Galaxies with the mass of the Milky Way dominate the stellar mass density of the Universe but it is uncertain how and when they were assembled. Here we study progenitors of these galaxies out to z=2.5, using data from the 3D-HST and CANDELS Treasury surveys. We find that galaxies with present-day stellar masses of log(M)~10.7 built ~90% of their stellar mass since z=2.5, with most of the star formation occurring before z=1. In marked contrast to the assembly history of massive elliptical galaxies, mass growth is not limited to large radii: the mass in the central 2 kpc of the galaxies increased by a factor of 3.2+-0.8 between z=2.5 and z=1. We therefore rule out simple models in which bulges were fully assembled at high redshift and disks gradually formed around them. Instead, bulges (and black holes) likely formed in lockstep with disks, through bar instabilities, migration, or other processes. We find that after z=1 the growth in the central regions gradually stopped and the disk continued to build, consistent with recent studies of the gas distributions in z~1 galaxies and the properties of many spiral galaxies today.

Satellite galaxies in rich clusters are subject to numerous physical processes that can significantly influence their evolution. However, the typical L* satellite galaxy resides in much lower mass galaxy groups, where the processes capable of altering their evolution are generally weaker and have had less time to operate. To investigate the extent to which satellite and central galaxy evolution differs, we separately model the stellar mass - halo mass (M* -Mh) relation for these two populations over the redshift interval 0 < z < 1. This relation for central galaxies is constrained by the galaxy stellar mass function while the relation for satellite galaxies is constrained against recent measurements of the galaxy two-point correlation function (2PCF). At z ~ 0 the satellites, on average, have ~10% larger stellar masses at fixed peak subhalo mass compared to central galaxies of the same halo mass. This is required in order to reproduce the observed stellar mass-dependent 2PCF and satellite fractions. At low masses our model slightly under-predicts the correlation function at ~1 Mpc scales. At z ~ 1 the satellite and central galaxy M*-Mh relations are consistent within the errors, and the model provides an excellent fit to the clustering data. At present, the errors on the clustering data at z ~ 2 are too large to constrain the satellite model. A simple model in which satellite and central galaxies share the same M*-Mh relation is able to reproduce the extant z ~ 2 clustering data. We speculate that the striking similarity between the satellite and central galaxy M*-Mh relations since z ~ 2 arises because the central galaxy relation evolves very weakly with time and because the stellar mass of the typical satellite galaxy has not changed significantly since it was accreted. [Abridged]

We present a robust method to constrain average galaxy star formation rates, star formation histories, and the intracluster light as a function of halo mass. Our results are consistent with observed galaxy stellar mass functions, specific star formation rates, and cosmic star formation rates from z=0 to z=8. We consider the effects of a wide range of uncertainties on our results, including those affecting stellar masses, star formation rates, and the halo mass function at the heart of our analysis. As they are relevant to our method, we also present new calibrations of the dark matter halo mass function, halo mass accretion histories, and halo-subhalo merger rates out to z=8. We also provide new compilations of cosmic and specific star formation rates; more recent measurements are now consistent with the buildup of the cosmic stellar mass density at all redshifts. Implications of our work include: halos near 10^12 Msun are the most efficient at forming stars at all redshifts, the baryon conversion efficiency of massive halos drops markedly after z ~ 2.5 (consistent with theories of cold-mode accretion), the ICL for massive galaxies is expected to be significant out to at least z ~ 1-1.5, and dwarf galaxies at low redshifts have higher stellar mass to halo mass ratios than previous expectations and form later than in most theoretical models. Finally, we provide new fitting formulae for star formation histories that are more accurate than the standard declining tau model. Our approach places a wide variety of observations relating to the star formation history of galaxies into a self-consistent framework based on the modern understanding of structure formation in LCDM. Constraints on the stellar mass-halo mass relationship and star formation rates are available for download at http://www.peterbehroozi.com/data.html .

Satellite galaxies in groups and clusters are more likely to have low star formation rates (SFR) and lie on the red-sequence than central (field) galaxies. Using galaxy group/cluster catalogs from SDSS DR7, together with a cosmological N-body simulation to track satellite orbits, we examine the star formation histories and quenching timescales of satellites of M_star > 5 x 10^9 M_sun at z=0. We first explore satellite infall histories: group preprocessing and ejected orbits are critical aspects of satellite evolution, and properly accounting for these, satellite infall typically occurred at z~0.5, or ~5 Gyr ago. To obtain accurate initial conditions for the SFRs of satellites at their time of first infall, we construct an empirical parametrization for the evolution of central galaxy SFRs and quiescent fractions. With this, we constrain the importance and efficiency of satellite quenching as a function of satellite and host halo mass, finding that satellite quenching is the dominant process for building up all quiescent galaxies at M_star < 10^10 M_sun. We then constrain satellite star formation histories, finding a 'delayed-then-rapid' quenching scenario: satellite SFRs evolve unaffected for 2-4 Gyr after infall, after which star formation quenches rapidly, with an e-folding time of < 0.8 Gyr. These quenching timescales are shorter for more massive satellites but do not depend on host halo mass: the observed increase in satellite quiescent fraction with halo mass arises simply because of satellites quenching in a lower mass group prior to infall (group preprocessing), which is responsible for up to half of quenched satellites in massive clusters. Because of the long time delay before quenching starts, satellites experience significant stellar mass growth after infall, nearly identical to central galaxies. This fact provides key physical insight into the subhalo abundance matching method.

The spectral energy distributions (SEDs) of galaxies are shaped by nearly every physical property of the system, including the star formation history, metal content, abundance pattern, dust mass, grain size distribution, star-dust geometry, and interstellar radiation field. The principal goal of stellar population synthesis (SPS) is to extract these variables from observed SEDs. In this review I provide an overview of the SPS technique and discuss what can be reliably measured from galaxy SEDs. Topics include stellar masses, star formation rates and histories, metallicities and abundance patterns, dust properties, and the stellar initial mass function.

The detailed abundance patterns of the stars within galaxies provide a unique window into the history of star formation (SF) at early times. Two widely used `chronometers' include the alpha and iron-peak elements, which are created on short and long timescales, respectively. These two clocks have been widely used to estimate SF timescales from moderate-resolution spectra of early-type galaxies. Elements formed via s-process neutron captures (e.g., Sr and Ba) comprise a third type of chronometer, as the site of the main s-process is believed to be intermediate and low-mass asymptotic giant branch stars. The [alpha/Ba] ratio in particular should provide a powerful new constraint on the SF histories of galaxies, in part because it is insensitive to the uncertain distribution of type Ia SNe detonation times and the overall Ia rate. Here we present new measurements of the abundance of Sr and Ba in nearby early-type galaxies by applying stellar population synthesis tools to high S/N optical spectra. We find a strong anti-correlation between [Mg/Fe] and [Ba/Fe], and a strong positive correlation between [Mg/Ba] and galaxy velocity dispersion. These trends are consistent with the idea that more massive galaxies formed their stars on shorter timescales compared to less massive galaxies, and rule out several other proposed explanations for the observed super-solar [Mg/Fe] values in massive galaxies. In contrast, [Sr/Fe]~0, with no strong variation across the sample. It is difficult to interpret the Sr trends without detailed chemical evolution models owing to the multiplicity of proposed nucleosynthetic sites for Sr.

Using reconstructed galaxy star formation histories, we calculate the instantaneous efficiency of galaxy star formation (i.e., the star formation rate divided by the baryon accretion rate) from $z=8$ to the present day. This efficiency exhibits a clear peak near a characteristic halo mass of 10^11.7 Msun, which coincides with longstanding theoretical predictions for the mass scale relevant to virial shock heating of accreted gas. Above the characteristic halo mass, the efficiency falls off as the mass to the minus four-thirds power; below the characteristic mass, the efficiency falls off at an average scaling of mass to the two-thirds power. By comparison, the shape and normalization of the efficiency change very little since z=4. We show that a time-independent star formation efficiency simply explains the shape of the cosmic star formation rate since z=4 in terms of dark matter accretion rates. The rise in the cosmic star formation from early times until z=2 is especially sensitive to galaxy formation efficiency. The mass dependence of the efficiency strongly limits where most star formation occurs, with the result that two-thirds of all star formation has occurred inside halos within a factor of three of the characteristic mass, a range that includes the mass of the Milky Way.

We present a simple model for the relationship between quasars, galaxies, and dark matter halos from 0.5<z<6. In the model, black hole (BH) mass is linearly related to galaxy mass, and galaxies are connected to dark matter halos via empirically constrained relations. A simple "scattered" light bulb model for quasars is adopted, wherein BHs shine at a fixed fraction of the Eddington luminosity during accretion episodes, and Eddington ratios are drawn from a lognormal distribution that is redshift-independent. This model has two free, physically meaningful parameters at each redshift: the normalization of the Mbh-Mgal relation and the quasar duty cycle; these parameters are fit to the observed quasar luminosity function (LF) over the interval 0.5<z<6. This simple model provides an excellent fit to the LF at all epochs, and also successfully predicts the observed projected two-point correlation of quasars from 0.5<z<2.5. It is significant that a single quasar duty cycle at each redshift is capable of reproducing the extant observations. The data are therefore consistent with a scenario wherein quasars are equally likely to exist in galaxies, and therefore dark matter halos, over a wide range in masses. The knee in the quasar LF is a reflection of the knee in the stellar mass-halo mass relation. Future constraints on the quasar LF and quasar clustering at high redshift will provide strong constraints on the model. In the model, the autocorrelation function of quasars becomes a strong function of luminosity only at the very highest luminosities, and will be difficult to observe because such quasars are so rare. Cross-correlation techniques may provide useful constraints on the bias of such rare objects.

The spectral absorption lines in early-type galaxies contain a wealth of information regarding the detailed abundance pattern, star formation history, and stellar initial mass function (IMF) of the underlying stellar population. Using our new population synthesis model that accounts for the effect of variable abundance ratios of 11 elements, we analyze very high quality absorption line spectra of 38 early-type galaxies and the nuclear bulge of M31. These data extend to 1um and they therefore include the IMF-sensitive spectral features NaI, CaII, and FeH at 0.82um, 0.86um and 0.99um, respectively. The models fit the data well, with typical rms residuals ~1%. Strong constraints on the IMF and therefore the stellar mass-to-light ratio, (M/L)_stars, are derived for individual galaxies. We find that the IMF becomes increasingly bottom-heavy with increasing velocity dispersion and [Mg/Fe]. At the lowest dispersions and [Mg/Fe] values the derived IMF is consistent with the Milky Way IMF, while at the highest dispersions and [Mg/Fe] values the derived IMF contains more low-mass stars (is more bottom-heavy) than even a Salpeter IMF. Our best-fit (M/L)_stars values do not exceed dynamically-based M/L values. We also apply our models to stacked spectra of four metal-rich globular clusters in M31 and find an (M/L)_stars that implies fewer low-mass stars than a Milky Way IMF, again agreeing with dynamical constraints. We discuss other possible explanations for the observed trends and conclude that variation in the IMF is the simplest and most plausible.

Due to its proximity (9 Mpc) and the strongly bimodal color distribution of its spectroscopically well-sampled globular cluster (GC) system, the early-type galaxy NGC 3115 provides one of the best available tests of whether the color bimodality widely observed in GC systems generally reflects a true metallicity bimodality. Color bimodality has alternatively been attributed to a strongly nonlinear color--metallicity relation reflecting the influence of hot horizontal branch stars. Here we couple Subaru Suprime-Cam gi photometry with Keck/DEIMOS spectroscopy to accurately measure GC colors and a CaT index that measures the CaII triplet. We find the NGC 3115 GC system to be unambiguously bimodal in both color and the CaT index. Using simple stellar population models, we show that the CaT index is essentially unaffected by variations in horizontal branch morphology over the range of metallicities relevant to GC systems (and is thus a robust indicator of metallicity) and confirm bimodality in the metallicity distribution. We assess the existing evidence for and against multiple metallicity subpopulations in early and late-type galaxies and conclude that metallicity bi/multimodality is common. We briefly discuss how this fundamental characteristic links directly to the star formation and assembly histories of galaxies.

All globular clusters (GCs) studied to date show evidence for internal variation in their light element abundances. These variations have been interpreted as evidence for multiple star formation episodes within GCs, with secondary episodes fueled, at least in part, by the ejecta of AGB stars from a first generation of stars. A major puzzle emerging from this otherwise plausible scenario is that the fraction of stars associated with the second episode of star formation is observed to be much larger than expected for a standard IMF. The present work investigates this tension by modeling the observed anti-correlation between [Na/Fe] and [O/Fe] for 20 Galactic GCs. If the abundance pattern of the retained AGB ejecta does not depend on GC mass at fixed [Fe/H], then a strong correlation is found between the fraction of current GC stellar mass comprised of pure AGB ejecta, f_p, and GC mass. This fraction varies from 0.20 at low masses (10^4.5 Msun) to 0.45 at high masses (10^6.5 Msun). The fraction of mass associated with pure AGB ejecta is directly related to the total mass of the cluster at birth; the ratio between the initial and present mass in stars can therefore be derived. Assuming a star formation efficiency of 50%, the observed Na-O anti-correlations imply that GCs were at least 10-20 times more massive at birth. These factors are lower limits because any mass-loss mechanism that removes first and second generation stars equally will leave f_p unchanged. The mass-dependence of f_p probably arises because lower mass GCs are unable to retain all of the AGB ejecta from the first stellar generation. Recent observations of elemental abundances in intermediate-age LMC clusters are re-interpreted and shown to be consistent with this basic scenario. A convincing explanation of these trends is currently lacking.

We present a new and innovative near-infrared multi-band ultraprecise spectroimager (NIMBUS) for SOFIA. This design is capable of characterizing a large sample of extrasolar planet atmospheres by measuring elemental and molecular abundances during primary transit and occultation. This wide-field spectroimager would also provide new insights into Trans-Neptunian Objects (TNO), Solar System occultations, brown dwarf atmospheres, carbon chemistry in globular clusters, chemical gradients in nearby galaxies, and galaxy photometric redshifts. NIMBUS would be the premier ultraprecise spectroimager by taking advantage of the SOFIA observatory and state of the art infrared technologies. This optical design splits the beam into eight separate spectral bandpasses, centered around key molecular bands from 1 to 4 microns. Each spectral channel has a wide field of view for simultaneous observations of a reference star that can decorrelate time-variable atmospheric and optical assembly effects, allowing the instrument to achieve ultraprecise calibration for imaging and photometry for a wide variety of astrophysical sources. NIMBUS produces the same data products as a low-resolution integral field spectrograph over a large spectral bandpass, but this design obviates many of the problems that preclude high-precision measurements with traditional slit and integral field spectrographs. This instrument concept is currently not funded for development.

The strength of gravity-sensitive absorption lines in the integrated light of old stellar populations is one of the few direct probes of the stellar initial mass function (IMF) outside of the Milky Way. Owing to the advent of fully depleted CCDs with little or no fringing it has recently become possible to obtain accurate measurements of these features. Here we present spectra covering the wavelength ranges 0.35 - 0.55 micron and 0.72 - 1.03 micron for the bulge of M31 and 34 early-type galaxies from the SAURON sample, obtained with the Low Resolution Imaging Spectrometer on Keck. The signal-to-noise ratio is >200 per Angstrom out to 1 micron, which is sufficient to measure gravity-sensitive features for individual galaxies and to determine how they depend on other properties of the galaxies. Combining the new data with previously obtained spectra for globular clusters in M31 and the most massive elliptical galaxies in the Virgo cluster we find that the dwarf-sensitive Na I doublet and the FeH Wing-Ford band increase systematically with velocity dispersion, while the giant-sensitive Ca II triplet decreases with dispersion. These trends are consistent with a varying IMF, such that galaxies with deeper potential wells have more dwarf-enriched mass functions. In a companion paper (Conroy & van Dokkum 2012) we use a comprehensive stellar population synthesis model to demonstrate that IMF effects can be separated from age and abundance variations and quantify the IMF variation among early-type galaxies.

Using galaxy group/cluster catalogs created from the Sloan Digital Sky Survey Data Release 7, we examine in detail the specific star formation rate (SSFR) distribution of satellite galaxies and its dependence on stellar mass, host halo mass, and halo-centric radius. All galaxies, regardless of central-satellite designation, exhibit a similar bimodal SSFR distribution, with a strong break at SSFR ~ 10^-11 yr^-1 and the same high SSFR peak; in no regime is there ever an excess of galaxies in the `green valley'. Satellite galaxies are simply more likely to lie on the quenched (`red sequence') side of the SSFR distribution. Furthermore, the satellite quenched fraction excess above the field galaxy value is nearly independent of galaxy stellar mass. An enhanced quenched fraction for satellites persists in groups with halo masses down to 3 x 10^11 Msol and increases strongly with halo mass and toward halo center. We find no detectable quenching enhancement for galaxies beyond ~2R_vir around massive clusters once these galaxies have been decomposed into centrals and satellites. These trends imply that (1) galaxies experience no significant environmental effects until they cross within ~R_vir of a more massive host halo, (2) after this, star formation in active satellites continues to evolve in the same manner as active central galaxies for several Gyrs, and (3) once begun, satellite star formation quenching occurs rapidly. These results place strong constraints on satellite-specific quenching mechanisms, as we will discuss further in companion papers.

Approximately 30% of all massive stars in the Galaxy are runaways with velocities exceeding 30 km/s. Their high speeds allow them to travel ~0.1-1 kpc away from their birth place before they explode at the end of their several Myr lifetimes. At high redshift, when galaxies were much smaller than in the local universe, runaways could venture far from the dense inner regions of their host galaxies. From these large radii, and therefore low column densities, much of their ionizing radiation is able to escape into the intergalactic medium. Runaways may therefore significantly enhance the overall escape fraction of ionizing radiation, fesc, from small galaxies at high redshift. We present simple models of the high-redshift runaway population and its impact on fesc as a function of halo mass, size, and redshift. We find that the inclusion of runaways enhances fesc by factors of ~1.1-8, depending on halo mass, galaxy geometry, and the mechanism of runaway production, implying that runaways may contribute 50-90% of the total ionizing radiation escaping from high-redshift galaxies. Runaways may therefore play an important role in reionizing the universe.

We present a comprehensive study of 250,000 galaxies targeted by the Baryon Oscillation Spectroscopic Survey (BOSS) up to z ~ 0.7 with the specific goal of identifying and characterising a population of galaxies that has evolved without significant merging. We compute a likelihood that each BOSS galaxy is a progenitor of the Luminous Red Galaxies (LRGs) sample, targeted by SDSS-I/II up z ~ 0.5, by using the fossil record of LRGs and their inferred star-formation histories, metallicity histories and dust content. We determine merger rates, luminosity growth rates and the evolution of the large-scale clustering between the two surveys, and we investigate the effect of using different stellar population synthesis models in our conclusions. We demonstrate that our sample is slowly evolving (of the order of 2 +/- 1.5% per Gyr by merging). Our conclusions refer to the bright and massive end of the galaxy population, with Mi0.55 < -22, and M* > 1E11.2 Msolar, corresponding roughly to 95% and 40% of the LRGs and BOSS galaxy populations, respectively. Our analysis further shows that any possible excess of flux in BOSS galaxies, when compared to LRGs, from potentially unresolved targets at z ~ 0.55 must be less than 1% in the r0.55-band (approximately equivalent to the g-band in the rest-frame of galaxies at z=0.55). When weighting the BOSS galaxies based on the predicted properties of the LRGs, and restricting the analysis to the reddest BOSS galaxies, we find an evolution of the large-scale clustering that is consistent with dynamical passive evolution, assuming a standard cosmology. We conclude that our likelihoods give a weighted sample that is as clean and as close to passive evolution (in dynamical terms, i.e. no or negligible merging) as possible, and that is optimal for cosmological studies.

We describe the design and data sample from the DEEP2 Galaxy Redshift Survey, the densest and largest precision-redshift survey of galaxies at z ~ 1 completed to date. The survey has conducted a comprehensive census of massive galaxies, their properties, environments, and large-scale structure down to absolute magnitude M_B = -20 at z ~ 1 via ~90 nights of observation on the DEIMOS spectrograph at Keck Observatory. DEEP2 covers an area of 2.8 deg^2 divided into four separate fields, observed to a limiting apparent magnitude of R_AB=24.1. Objects with z < 0.7 are rejected based on BRI photometry in three of the four DEEP2 fields, allowing galaxies with z > 0.7 to be targeted ~2.5 times more efficiently than in a purely magnitude-limited sample. Approximately sixty percent of eligible targets are chosen for spectroscopy, yielding nearly 53,000 spectra and more than 38,000 reliable redshift measurements. Most of the targets which fail to yield secure redshifts are blue objects that lie beyond z ~ 1.45. The DEIMOS 1200-line/mm grating used for the survey delivers high spectral resolution (R~6000), accurate and secure redshifts, and unique internal kinematic information. Extensive ancillary data are available in the DEEP2 fields, particularly in the Extended Groth Strip, which has evolved into one of the richest multiwavelength regions on the sky. DEEP2 surpasses other deep precision-redshift surveys at z ~ 1 in terms of galaxy numbers, redshift accuracy, sample number density, and amount of spectral information. We also provide an overview of the scientific highlights of the DEEP2 survey thus far. This paper is intended as a handbook for users of the DEEP2 Data Release 4, which includes all DEEP2 spectra and redshifts, as well as for the publicly-available DEEP2 DEIMOS data reduction pipelines. [Abridged]

Low-mass stars (M<0.4Msun) are thought to comprise the bulk of the stellar mass of galaxies but they constitute only of order a percent of the bolometric luminosity of an old stellar population. Directly estimating the number of low-mass stars from integrated flux measurements of old stellar systems is therefore possible but very challenging given the numerous variables that can affect the light at the percent level. Here we present a new population synthesis model created specifically for the purpose of measuring the low-mass initial mass function (IMF) down to ~0.1Msun for metal-rich stellar populations with ages in the range 3-13.5 Gyr. Our fiducial model is based on the synthesis of three separate isochrones and a combination of optical and near-IR empirical stellar libraries in order to produce integrated light spectra over the wavelength interval 0.35mu<lambda<2.4mu at a resolving power of R~2000. New synthetic stellar atmospheres and spectra have been computed in order to model the spectral variations due to changes in individual elemental abundances including C, N, Na, Mg, Si, Ca, Ti, Fe, and generic alpha elements. We demonstrate the power of combining blue spectral features with surface gravity-sensitive near-IR features in order to simultaneously constrain the low-mass IMF, stellar population age, metallicity, and abundance pattern from integrated light measurements. Finally, we show that the shape of the low-mass IMF can also be directly constrained by employing a suite of surface gravity-sensitive spectral features, each of which is most sensitive to a particular mass interval.

We investigate the evolution of the MBH-{\sigma} relation by examining the relationship between the intrinsic scatter in the MBH-{\sigma} relation and galaxy bolometric nuclear luminosity, the latter being a probe of the accretion rate of the black hole (BH). Our sample is composed of galaxies with classical bulges when possible, of which 38 have dynamically measured BHs masses, and 17 have BHs masses measured by reverberation mapping. In order to obtain the bolometric nuclear luminosity for galaxies with low nuclear luminosity, we convert the X-ray nuclear luminosity measured by Chandra to bolometric luminosity. We find that the scatter in the MBH-{\sigma} relation is uncorrelated with nuclear luminosity over seven orders of magnitude in luminosity, with the high luminosity end approaching the Eddington luminosity. This suggests that at the present epoch galaxies evolve along the MBH-{\sigma} relation. This conclusion is consistent with the standard paradigm that BHs grow contemporaneously with their host stellar spheroids.

Measurements of the total amount of stars locked up in galaxies as a function of host halo mass contain key clues about the efficiency of processes that regulate star formation. We derive the total stellar mass fraction f_star as a function of halo mass M500c from z=0.2 to z=1 using two complementary methods. First, we derive f_star using a statistical Halo Occupation Distribution model jointly constrained by data from lensing, clustering, and the stellar mass function. This method enables us to probe f_star over a much wider halo mass range than with group or cluster catalogs. Second, we derive f_star at group scales using a COSMOS X-ray group catalog and we show that the two methods agree to within 30%. We quantify the systematic uncertainty on f_star using abundance matching methods and we show that the statistical uncertainty on f_star (~10%) is dwarfed by systematic uncertainties associated with stellar mass measurements (~45% excluding IMF uncertainties). Assuming a Chabrier IMF, we find 0.012<f_star<0.025 at M500c=10^13 Msun and 0.0057<f_star<0.015 at M500c=10^14 Msun. These values are significantly lower than previously published estimates. We investigate the cause of this difference and find that previous work has overestimated f_star due to a combination of inaccurate stellar mass estimators and/or because they have assumed that all galaxies in groups are early type galaxies with a constant mass-to-light ratio. Contrary to previous claims, our results suggest that the mean value of f_star is always significantly lower than f_gas for halos above 10^13 Msun. Combining our results with recently published gas mas fractions, we find a shortfall in f_star+f_gas at R500c compared to the cosmic mean. This shortfall varies with halo mass and becomes larger towards lower halos masses.

We recently found that massive cluster elliptical galaxies have strong Na I 8183,8195 and FeH 9916 Wing-Ford band absorption, indicating the presence of a very large population of stars with masses <~ 0.3 M_sun. Here we test this result by comparing the elliptical galaxy spectra to those of luminous globular clusters associated with M31. These globular clusters have similar metallicities, abundance ratios and ages as massive elliptical galaxies but their low dynamical mass-to-light ratios rule out steep stellar initial mass functions (IMFs). From high quality Keck spectra we find that the dwarf-sensitive absorption lines in globular clusters are significantly weaker than in elliptical galaxies, and consistent with normal IMFs. The differences in the Na I and Wing-Ford indices are 0.027 +- 0.007 mag and 0.017 +- 0.006 mag respectively. We directly compare the two classes of objects by subtracting the averaged globular cluster spectrum from the averaged elliptical galaxy spectrum. The difference spectrum is well fit by the difference between a stellar population synthesis model with a bottom-heavy IMF and one with a bottom-light IMF. We speculate that the slope of the IMF may vary with velocity dispersion, although it is not yet clear what physical mechanism would be responsible for such a relation.

The conjecture that the ancient globular clusters (GCs) formed at the center of their own dark matter halos was first proposed by Peebles (1984), and has recently been revived to explain the puzzling abundance patterns observed within many GCs. In this paper we demonstrate that the outer stellar density profile of isolated GCs is very sensitive to the presence of an extended dark halo. The GCs NGC 2419, located at 90 kpc from the center of our Galaxy, and MGC1, located at ~200 kpc from the center of M31, are ideal laboratories for testing the scenario that GCs formed at the centers of massive dark halos. Comparing analytic models to observations of these GCs, we conclude that these GCs cannot be embedded within dark halos with a virial mass greater than 10^6 Msun, or, equivalently, the dark matter halo mass-to-stellar mass ratio must be Mdm/M_*<1. If these GCs have indeed orbited within weak tidal fields throughout their lifetimes, then these limits imply that these GCs did not form within their own dark halos. Recent observations of an extended stellar halo in the GC NGC 1851 are also interpreted in the context of our analytic models. Implications of these results for the formation of GCs are briefly discussed.

We improve estimates of stellar mass and mass-weighted average age of Type Ia supernova (SN Ia) host galaxies by combining UV and near-IR photometry with optical photometry in our analysis. Using 206 SNe Ia drawn from the full three-year SDSS-II Supernova Survey (median redshift of z {\approx} 0.2) and multi-wavelength host-galaxy photometry from SDSS, GALEX, and UKIDSS, we present evidence of a correlation (1.9{\sigma} confidence level) between the residuals of SNe Ia about the best-fit Hubble relation and the mass-weighted average age of their host galaxies. The trend is such that older galaxies host SNe Ia that are brighter than average after standard light-curve corrections are made. We also confirm, at the 3.0{\sigma} level, the trend seen by previous studies that more massive galaxies often host brighter SNe Ia after light-curve correction.

The properties of dark matter halos, including mass growth, correlate with larger scale environment at fixed mass, an effect known as assembly bias. However, whether this environmental dependence manifests itself in galaxy properties remains unclear. We apply a group-finding algorithm to DR7 of the SDSS to estimate the halo mass of each galaxy and to decompose galaxies into those that exist at the centers of distinct halos and those that orbit as satellites within larger halos. Using the 4000-A break as a measure of star formation history, we examine the correlation between the quenched fraction of galaxies, f_q, and large-scale environment, rho. At all galaxy magnitudes, there is a positive, monotonic relationship between f_q and rho. We use the group catalog to decompose this correlation into the contribution from central and satellite galaxies as a function of halo mass. Because satellites are more likely to be quenched than central galaxies, the observed f_q-rho correlation is primarily due to variations of the halo mass function with environment, which causes a larger fraction of satellite galaxies at high rho. For low-mass central galaxies (Mgal <~ 10^10.0 Msol/h^2), there is no correlation between f_q and rho. These results are inconsistent with the strong assembly bias of dark matter halos seen in this mass regime if recent galaxy growth at all correlates with recent halo growth, as we demonstrate through a high resolution N-body simulation. We also find that the mean stellar age of quenched central galaxies is independent of rho at fixed Mgal, while the formation times of low mass halos vary significantly. We conclude that the processes that halt the star formation of low mass central galaxies are not correlated to the formation histories of their host halos, and old galaxies do not reside preferentially in old halos. (Abridged)

We investigate the origin of the relations between stellar mass and optical circular velocity for early-type (ETG) and late-type (LTG) galaxies --- the Faber-Jackson (FJ) and Tully-Fisher (TF) relations. We combine measurements of dark halo masses (from satellite kinematics and weak lensing), and the distribution of baryons in galaxies (from a new compilation of galaxy scaling relations), with constraints on dark halo structure from cosmological simulations. The principle unknowns are the halo response to galaxy formation and the stellar initial mass function (IMF). The slopes of the TF and FJ relations are naturally reproduced for a wide range of halo response and IMFs. However, models with a universal IMF and universal halo response cannot simultaneously reproduce the zero points of both the TF and FJ relations. For a model with a universal Chabrier IMF, LTGs require halo expansion, while ETGs require halo contraction. A Salpeter IMF is permitted for high mass (sigma > 180 km/s) ETGs, but is inconsistent for intermediate masses, unless V_circ(R_e)/sigma_e > 1.6. If the IMF is universal and close to Chabrier, we speculate that the presence of a major merger may be responsible for the contraction in ETGs while clumpy accreting streams and/or feedback leads to expansion in LTGs. Alternatively, a recently proposed variation in the IMF disfavors halo contraction in both types of galaxies. Finally we show that our models naturally reproduce flat and featureless circular velocity profiles within the optical regions of galaxies without fine-tuning.

Stellar population synthesis (SPS) models are a key ingredient of many galaxy evolution studies. Unfortunately, the models are still poorly calibrated for certain stellar evolution stages. Of particular concern is the treatment of the thermally-pulsing asymptotic giant branch (TP-AGB) phase, as different implementations lead to systematic differences in derived galaxy properties. Post-starburst galaxies are a promising calibration sample, as TP-AGB stars are thought to be most prominently visible during this phase. Here, we use post-starburst galaxies in the NEWFIRM medium-band survey to assess different SPS models. The available photometry allows the selection of a homogeneous and well-defined sample of 62 post-starburst galaxies at 0.7<z<2.0, from which we construct a well-sampled composite spectral energy distribution (SED) over the range 1200-40 000 Angstrom. The SED is well-fit by the Bruzual & Charlot SPS models, while the Maraston models do not reproduce the rest-frame optical and near-infrared parts of the SED simultaneously. When the fitting is restricted to lambda < 6000 Angstrom, the Maraston models overpredict the near-infrared luminosity, implying that these models give too much weight to TP-AGB stars. Using the flexible SPS models by Conroy et al, and assuming solar metallicity, we find that the contribution of TP-AGB stars to the integrated SED is a factor of ~3 lower than predicted by the latest Padova TP-AGB models. Whether this is due to lower bolometric luminosities, shorter lifetimes, and/or heavy dust obscuration of TP-AGB stars remains to be addressed. Altogether, our data demand a low contribution from TP-AGB stars to the SED of post-starburst galaxies.

Nearly all globular clusters (GCs) studied to date show evidence for multiple stellar populations, in stark contrast to the conventional view that GCs are a mono-metallic, coeval population of stars. Building on earlier work, we propose a simple physical model for the early evolution (several 10^8 yr) of GCs. We consider the effects of stellar mass-loss, type II and prompt type Ia supernovae, ram pressure, and accretion from the ambient ISM on the development of a young GC's own gas reservoir. In our model, type II SNe from a first generation of star formation clears the GC of its initial gas reservoir. Over the next several 10^8 yr, mass lost from AGB stars and matter accreted from the ambient ISM collect at the center of the GC. This material must remain quite cool (T~10^2K), but does not catastrophically cool on a crossing time because of the high Lyman-Werner flux density in young GCs. The collection of gas within the GC must compete with ram pressure from the ambient ISM. After several 10^8 yr, the Lyman-Werner photon flux density drops by more than three orders of magnitude, allowing molecular hydrogen and then stars to form. After this second generation of star formation, type II SNe from the second generation and then prompt type Ia SNe associated with the first generation maintain a gas-free GC, thereby ending the cycle of star formation events. Our model makes clear predictions for the presence or absence of multiple stellar populations within GCs as a function of GC mass and formation environment. Analyzing intermediate-age LMC clusters, we find evidence for a mass threshold of ~10^4 Msun below which LMC clusters appear to be truly coeval. This threshold mass is consistent with our predictions for the mass at which ram pressure is capable of clearing gas from clusters in the LMC at the present epoch. (ABRIDGED)

The stellar initial mass function (IMF) describes the mass distribution of stars at the time of their formation and is of fundamental importance for many areas of astrophysics. The IMF is reasonably well constrained in the disk of the Milky Way but we have very little direct information on the form of the IMF in other galaxies and at earlier cosmic epochs. Here we investigate the stellar mass function in elliptical galaxies by measuring the strength of the Na I doublet and the Wing-Ford molecular FeH band in their spectra. These lines are strong in stars with masses <0.3 Msun and weak or absent in all other types of stars. We unambiguously detect both signatures, consistent with previous studies that were based on data of lower signal-to-noise ratio. The direct detection of the light of low mass stars implies that they are very abundant in elliptical galaxies, making up >80% of the total number of stars and contributing >60% of the total stellar mass. We infer that the IMF in massive star-forming galaxies in the early Universe produced many more low mass stars than the IMF in the Milky Way disk, and was probably slightly steeper than the Salpeter form in the mass range 0.1 - 1 Msun.

The attenuation of starlight by interstellar dust is investigated in a sample of low redshift, disk-dominated star-forming galaxies using photometry from GALEX and SDSS. By considering broadband colors as a function of galaxy inclination we are able to confidently separate trends arising from increasing dust opacity from possible differences in stellar populations, since stellar populations do not correlate with inclination. All commonly employed dust attenuation curves (such as the Calzetti curve for starbursts, or a power-law curve) provide poor fits to the ultraviolet colors for moderately and highly inclined galaxies. This conclusion rests on the fact that the average FUV-NUV color varies little from face-on to edge-on galaxies, while other colors such as NUV-u and u-r vary strongly with inclination. After considering a number of model variations, we are led to speculate that the presence of the strong dust extinction feature at 2175A seen in the Milky Way (MW) extinction curve is responsible for the observed trends. Independent of our interpretation, these results imply that the modeling of dust attenuation in the ultraviolet is significantly more complicated than traditionally assumed. These results also imply a very weak dependence of the FUV-NUV color on total FUV attenuation, and we conclude from this that it is extremely difficult to use only the observed UV spectral slope to infer the total UV dust attenuation, as is commonly done. We propose several simple tests that might finally identify the grain population responsible for the 2175A feature.

We conduct a comprehensive analysis of the relationship between central galaxies and their host dark matter halos, as characterized by the stellar mass-halo mass (SM-HM) relation, with rigorous consideration of uncertainties. Our analysis focuses on results from the abundance matching technique, which assumes that every dark matter halo or subhalo above a specific mass threshold hosts one galaxy. We discuss the quantitative effects of uncertainties in observed galaxy stellar mass functions (GSMFs) (including stellar mass estimates and counting uncertainties), halo mass functions (including cosmology and uncertainties from substructure), and the abundance matching technique used to link galaxies to halos (including scatter in this connection). Our analysis results in a robust estimate of the SM-HM relation and its evolution from z=0 to z=4. The shape and evolution are well constrained for z < 1. The largest uncertainties at these redshifts are due to stellar mass estimates; however, failure to account for scatter in stellar masses at fixed halo mass can lead to errors of similar magnitude in the SM-HM relation for central galaxies in massive halos. We also investigate the SM-HM relation to z=4, although the shape of the relation at higher redshifts remains fairly unconstrained when uncertainties are taken into account. These results will provide a powerful tool to inform galaxy evolution models. [Abridged]

Using estimates of dark halo masses from satellite kinematics, weak gravitational lensing, and halo abundance matching, combined with the Tully-Fisher and Faber-Jackson relations, we derive the mean relation between the optical, V_opt, and virial, V_200, circular velocities of early- and late-type galaxies at redshift z~0. For late-type galaxies V_opt ~ V_200 over the velocity range V_opt=90-260 km/s, and is consistent with V_opt = V_maxh (the maximum circular velocity of NFW dark matter haloes in the concordance LCDM cosmology). However, for early-type galaxies V_opt \ne V_200, with the exception of early-type galaxies with V_opt simeq 350 km/s. This is inconsistent with early-type galaxies being, in general, globally isothermal. For low mass (V_opt < 250 km/s) early-types V_opt > V_maxh, indicating that baryons have modified the potential well, while high mass (V_opt > 400 km/s) early-types have V_opt < V_maxh. Folding in measurements of the black hole mass - velocity dispersion relation, our results imply that the supermassive black hole - halo mass relation has a logarithmic slope which varies from ~1.4 at halo masses of ~10^{12} Msun/h to ~0.65 at halo masses of 10^{13.5} Msun/h. The values of V_opt/V_200 we infer for the Milky Way and M31 are lower than the values currently favored by direct observations and dynamical models. This offset is due to the fact that the Milky Way and M31 have higher V_opt and lower V_200 compared to typical late-type galaxies of the same stellar masses. We show that current high resolution cosmological hydrodynamical simulations are unable to form galaxies which simultaneously reproduce both the V_opt/V_200 ratio and the V_opt-M_star (Tully-Fisher/Faber-Jackson) relation.

A novel technique is employed for estimating attenuation curves in galaxies where only photometry and spectroscopic redshifts are available. This technique provides a powerful measure of particular extinction features such as the UV bump at 2175\A, which has been observed in environments ranging from the Milky Way to high-redshift star-forming galaxies. Knowledge of the typical strength of the UV bump as a function of environment and redshift is crucial for converting restframe UV flux into star formation rates. The UV bump will impart a unique signature as it moves through various filters due to redshifting; its presence can therefore be disentangled from other stellar population effects. The utility of this technique is demonstrated with a large sample of galaxies drawn from the DEEP2 Galaxy Redshift Survey. The observed B-R color of star-forming galaxies at 0.6<z<1.4 disfavors the presence of a UV bump as strong as observed in the Milky Way, and instead favors restframe UV (1800A<lambda<3000A) attenuation curves similar to the Milky Way without a UV bump or a power-law with index delta=-0.7. Stronger constraints on the strength of the UV bump in galaxies can be achieved if independent constraints on the V-band optical depth are available.

Models for the formation and evolution of galaxies readily predict physical properties such as the star formation rates, metal enrichment histories, and, increasingly, gas and dust content of synthetic galaxies. Such predictions are frequently compared to the spectral energy distributions of observed galaxies via the stellar population synthesis (SPS) technique. Substantial uncertainties in SPS exist, and yet their relevance to the task of comparing galaxy evolution models to observations has received little attention. In the present work we begin to address this issue by investigating the importance of uncertainties in stellar evolution, the initial stellar mass function (IMF), and dust and interstellar medium (ISM) properties on the translation from models to observations. We demonstrate that these uncertainties translate into substantial uncertainties in the ultraviolet, optical, and near-infrared colors of synthetic galaxies. Aspects that carry significant uncertainties include the logarithmic slope of the IMF above 1 Msun, dust attenuation law, molecular cloud disruption timescale, clumpiness of the ISM, fraction of unobscured starlight, and treatment of advanced stages of stellar evolution including blue stragglers, the horizontal branch, and the thermally-pulsating asymptotic giant branch. The interpretation of the resulting uncertainties in the derived colors is highly non-trivial because many of the uncertainties are likely systematic, and possibly correlated with the physical properties of galaxies. We therefore urge caution when comparing models to observations.

Stellar population synthesis (SPS) provides the link between the stellar and dust content of galaxies and their observed spectral energy distributions. In the present work we perform a comprehensive calibration of our own flexible SPS (FSPS) model against a suite of data. Several public SPS models are intercompared, including the models of Bruzual & Charlot (BC03), Maraston (M05) and FSPS. The relative strengths and weaknesses of these models are evaluated, with the following conclusions: 1) The FSPS and BC03 models compare favorably with MC data at all ages, whereas M05 colors are too red and the age-dependence is incorrect; 2) All models yield similar optical and near-IR colors for old metal-poor systems, and yet they all provide poor fits to the integrated J-K and V-K colors of both MW and M31 star clusters; 4) All models predict ugr colors too red, D4000 strengths too strong and Hdelta strengths too weak compared to massive red sequence galaxies, under the assumption that such galaxies are composed solely of old metal-rich stars; 5) FSPS and, to a lesser extent, BC03 can reproduce the optical and near-IR colors of post-starburst galaxies, while M05 cannot. Reasons for these discrepancies are explored. The failure at predicting the ugr colors, D4000, and Hdelta strengths can be explained by some combination of a minority population of metal-poor stars, young stars, blue straggler and/or blue horizontal branch stars, but not by appealing to inadequacies in either theoretical stellar atmospheres or canonical evolutionary phases (e.g., the main sequence turn-off). We emphasize that due to a lack of calibrating star cluster data in regions of the metallicity-age plane relevant for galaxies, all of these models continue to suffer from serious uncertainties that are difficult to quantify. (ABRIDGED)

A simple, observationally-motivated model is presented for understanding how halo masses, galaxy stellar masses, and star formation rates are related, and how these relations evolve with time. The relation between halo mass and galaxy stellar mass is determined by matching the observed spatial abundance of galaxies to the expected spatial abundance of halos at multiple epochs -- i.e. more massive galaxies are assigned to more massive halos at each epoch. Halos at different epochs are connected by halo mass accretion histories estimated from N-body simulations. The halo--galaxy connection at fixed epochs in conjunction with the connection between halos across time provides a connection between observed galaxies across time. With approximations for the impact of merging and accretion on the growth of galaxies, one can then directly infer the star formation histories of galaxies as a function of stellar and halo mass. This model is tuned to match both the observed evolution of the stellar mass function and the normalization of the observed star formation rate -- stellar mass relation to z~1. The data demands, for example, that the star formation rate density is dominated by galaxies with Mstar~10^{10.0-10.5} Msun from 0<z<1, and that such galaxies over these epochs reside in halos with M~10^{11.5-12.5} Msun. The star formation rate -- halo mass relation is approximately Gaussian over the range 0<z<1 with a mildly evolving mean and normalization. This model is then used to shed light on a number of issues, including 1) a clarification of ``downsizing'', 2) the lack of a sharp characteristic halo mass at which star formation is truncated, and 3) the dominance of star formation over merging to the stellar build-up of galaxies with Mstar<10^11 Msun at z<1.

The stellar masses, mean ages, metallicities, and star formation histories of galaxies are now commonly estimated via stellar population synthesis (SPS) techniques. SPS relies on stellar evolution calculations from the main sequence to stellar death, stellar spectral libraries, phenomenological dust models, and stellar initial mass functions (IMFs). The present work is the first in a series that explores the impact of uncertainties in key phases of stellar evolution and the IMF on the derived physical properties of galaxies and the expected luminosity evolution for a passively evolving set of stars. A Monte-Carlo Markov-Chain approach is taken to fit near-UV through near-IR photometry of a representative sample of low- and high-redshift galaxies with this new SPS model. Significant results include the following: 1) including uncertainties in stellar evolution, stellar masses at z~0 carry errors of ~0.3 dex at 95% CL with little dependence on luminosity or color, while at z~2, the masses of bright red galaxies are uncertain at the ~0.6 dex level; 2) either current stellar evolution models, current observational stellar libraries, or both, do not adequately characterize the metallicity-dependence of the thermally-pulsating asymptotic giant branch phase; 3) conservative estimates on the uncertainty of the slope of the IMF in the solar neighborhood imply that luminosity evolution per unit redshift is uncertain at the ~0.4 mag level in the K-band, which is a substantial source of uncertainty for interpreting the evolution of galaxy populations across time; 4) The more plausible assumption of a distribution of stellar metallicities, rather than a fixed value as is usually assumed, can have significant effects on the interpretation of colors blueward of the V-band. (ABRIDGED)

We present a simple formalism to interpret two galaxy statistics, the UV luminosity function and two-point correlation functions for star-forming galaxies at z~4, 5, 6 in the context of LCDM cosmology. Both statistics are the result of how star formation takes place in DM halos, and thus are used to constrain how UV light depends on halo properties such as mass. The two measures were taken from the GOODS data, thus ideal for joint analysis. The two physical quantities we explore are the SF duty cycle, and the range of L_UV that a halo of mass M can have (mean and variance). The former addresses the typical duration of SF activity in halos while the latter addresses the averaged SF history and regularity of gas inflow into these systems. We explore various physical models consistent with data, and find the following: 1) the typical duration of SF observed in the data is <0.4 Gyr (1 sig), 2) the inferred scaling law between L_UV and halo mass M from the observed slope of the LFs is roughly linear at all redshifts, and 3) L_UV for a fixed halo mass decreases with time, implying that the SF efficiency (after dust extinction) is higher at earlier times. We explore several physical scenarios relating star formation to halo mass, but find that these scenarios are indistinguishable due to the limited range of halo mass probed by our data. In order to discriminate between different scenarios, we discuss constraining the bright-faint galaxy cross-correlation functions and luminosity-dependence of galaxy bias. (Abridged)

Groups and clusters of galaxies occupy a special position in the hierarchy of large-scale cosmic structures because they are the largest and the most massive (from ~10^13 Msun to over 10^15 Msun) objects in the universe that have had time to undergo gravitational collapse. The large masses of clusters imply that their contents have been accreted from regions of ~8-40 comoving Mpc in size and should thus be representative of the mean matter content of the universe. During the next decade sensitive multi-wavelength observations should be able to map the radial distributions of all main mass components (stars, cold, warm, and hot gas and total mass) at z<~ 1 out to the virial radius. At the same time, comparative studies of real and simulated cluster samples sould allow us to use clusters as veritable astrophysical laboratories for studying galaxy formation, as well as testing our theoretical models of structure formation and underlying assumptions about fundamental physics governing the universe.

We use high-resolution N-body simulations, combined with a halo occupation model of galaxy bias, to investigate voids in the galaxy distribution. Our goal is to address the 'void phenomenon' of Peebles (2001), which presents the observed dearth of faint galaxies in voids as a challenge to the current cosmology. In our model, galaxy luminosity is determined only as a function of dark matter halo mass. With this simple assumption, we demonstrate that large, empty voids of ~15 Mpc/h in diameter are expected even for galaxies seven magnitudes fainter than L*. The predictions of our model are in excellent agreement with several statistical measures; (i) the luminosity function of galaxies in underdense regions, (ii) nearest neighbor statistics of dwarf galaxies, (iii) the void probability function of faint galaxies. In the transition between filaments and voids in the dark matter, the halo mass function changes abruptly, causing the maximum galaxy luminosity to decrease by ~5 magnitudes over a range of ~1 Mpc/h. Thus the boundary between filaments and voids in the galaxy distribution is nearly as sharp for dwarfs as for ~L* objects. These results support a picture in which galaxy formation is driven predominantly by the mass of the host dark matter halo, and is nearly independent of the larger-scale halo environment. Further, they demonstrate that LCDM, combined with a straightforward bias model, naturally explains the existence of the void phenomenon.

A variety of physical heating mechanisms are combined with radiative cooling to explore, via one dimensional hydrodynamic simulations, the expected thermal properties of the intracluster medium (ICM) in the context of the cooling flow problem. Energy injection from type Ia supernovae, thermal conduction, and dynamical friction (DF) from orbiting satellite galaxies are considered. The novel feature of this work is the exploration of a wide range of efficiencies of each heating process. While the latter two can provide a substantial amount of energy, neither mechanism operating alone can produce nor maintain an ICM in thermal balance over cosmological timescales, in stark contrast with observations. For simulated clusters with initially isothermal temperature profiles, both mechanisms acting in combination result in long-term thermal balance for a range of ICM temperatures and for central electron densities less than n_e~0.02 cm^-3; at greater densities catastrophic cooling invariably occurs. Furthermore, these heating mechanisms can neither produce nor maintain clusters with a declining temperature profile in the central regions, implying that the observed "cooling-core'' clusters, which have such declining temperature profiles, cannot be maintained with these mechanisms alone. Thus, while there appears to be an abundant supply of energy capable of heating the ICM in clusters, it is extremely difficult for the energy deposition to occur in such a way that the ICM remains in thermal balance over cosmological time-scales. These results strongly suggest that a more dynamic heating process such as feedback from a central black hole is required to generate the properties of observed intracluster media. (ABRIDGED)

Star-forming galaxies constitute the majority of galaxies with stellar masses >10^10 M_Sun/h^2 at z~2 and dominate the star-formation rate density of the Universe at this early epoch. It is thus critical to understand their origins, evolution, and connection to the underlying dark matter distribution. To this end, we identify the dark matter halos (including subhalos) that are likely to contain star-forming galaxies at z~2 (z2SFGs) within a large dissipationless cosmological simulation and then use halo merger histories to follow the evolution of z2SFG descendants to z~1 and z~0. The evolved halos at these epochs are then confronted with an array of observational data in order to uncover the likely descendants of z2SFGs. Though the evolved halos have clustering strengths comparable to red galaxies at z~1 and z~0, we find that the bulk of z2SFGs do not evolve into red galaxies, at either epoch. This conclusion is based primarily on the fact that the space density of z2SFGs is much higher than that of lower redshift red galaxies, even when accounting for the merging of z2SFG descendants, which decreases the number density of z2SFG descendants by at most a factor of two by z~0. Of the ~50% of z2SFGs that survive to z~0, ~70% reside at the center of z~0 dark matter halos with M>10^12 M_Sun/h. Halo occupation modeling of z~0 galaxies suggests that such halos are occupied by galaxies with M_r<-20.5, implying that these z2SFGs evolve into ``typical'' ~ L* galaxies today, including our own Galaxy. The remaining ~30% become satellite galaxies by z~0, and comparison to halo occupation modeling suggests that they are rather faint, with M_r<-19.5. (ABRIDGED)

We use measurements of the projected galaxy correlation function w_p and galaxy void statistics to test whether the galaxy content of halos of fixed mass is systematically different in low density environments. We present new measurements of the void probability function (VPF) and underdensity probability function (UPF) from Data Release Four of the Sloan Digital Sky Survey, as well as new measurements of the VPF from the full data release of the Two-Degree Field Galaxy Redshift Survey. We compare these measurements to predictions calculated from models of the Halo Occupation Distribution (HOD) that are constrained to match both w_p and the space density of galaxies. The standard implementation of the HOD assumes that galaxy occupation depends on halo mass only, and is independent of local environment. For luminosity-defined samples, we find that the standard HOD prediction is a good match to the observations, and the data exclude models in which galaxy formation efficiency is reduced in low-density environments. More remarkably, we find that the void statistics of red and blue galaxies (at L ~ 0.4L_*) are perfectly predicted by standard HOD models matched to the correlation function of these samples, ruling out "assembly bias" models in which galaxy color is correlated with large-scale environment at fixed halo mass. We conclude that the luminosity and color of field galaxies are determined predominantly by the mass of the halo in which they reside and have little direct dependence on the environment in which the host halo formed. In broader terms, our results show that the sizes and emptiness of voids found in the distribution of L > 0.2L_* galaxies are in excellent agreement with the predictions of a standard cosmological model with a simple connection between galaxies and dark matter halos. (abridged)

In the LCDM cosmology dark matter halos grow primarily through the accretion of smaller halos. Much of the mass in a halo of 10^14Msun comes in through accretion of ~10^13Msun halos. If each such halo hosted one luminous red galaxy (LRG) then the accretion of so many halos is at odds with the observed number of LRGs in clusters unless these accreted LRGs merge or disrupt on relatively short timescales (~2Gyr). These timescales are consistent with classical dynamical friction arguments, and imply that 2-3 LRGs have merged or disrupted within each halo more massive than 10^14Msun by z=0. The total amount of stellar mass brought into these massive halos by z=0 is consistent with observations once the intracluster light (ICL) is included. If disrupted LRGs build up the ICL, then the hierarchical growth of massive halos implies that a substantial amount of ICL should also surround satellite LRGs, as suggested by recent observations of the Virgo cluster. Finally, we point out that these results are entirely consistent with a non-evolving clustering strength and halo occupation distribution, and note that observations of the latter in fact support the hypothesis that merging/disruption of massive galaxies does indeed take place at late times.

We use a set of simulation-based models for the dissipationless evolution of galaxies since z=1 to constrain the fate of accreted satellites embedded in dark matter subhalos. These models assign stellar mass to dark matter halos at z=1 by relating the observed galaxy stellar mass function (GSMF) to the halo+subhalo mass function monotonically. The evolution of the stellar mass content is then followed using halo merger trees extracted from N-body simulations. Our models are differentiated only in the fate assigned to satellite galaxies once subhalos, within which satellites are embedded, disrupt. These models are confronted with the observed evolution in the massive end of the GSMF, the z~0 brightest cluster galaxy (BCG)-cluster mass relation, and the combined BCG and intracluster light (ICL) luminosity distribution -- all observables expected to evolve approximately dissipationlessly since z=1. The combined observational constraints favor a model in which the vast majority (>80%) of satellite stars from disrupted subhalos go into the ICL. Conversely, models that leave behind a significant population of satellite galaxies once the subhalo has disrupted are strongly disfavored, as are models that put a significant fraction of satellite stars into the BCG. Our results show that observations of the ICL provide useful and unique constraints on models of galaxy merging and the dissipationless evolution of galaxies in groups and clusters.

We measure the evolution in the virial mass-to-light ratio (M_{200}/L_B) and virial-to-stellar mass ratio (M_{200}/M_\ast) for isolated ~ L* galaxies between z~1 and z~0 by combining data from the DEEP2 Galaxy Redshift Survey and the Sloan Digital Sky Survey. Utilizing the motions of satellite galaxies around isolated galaxies, we measure line-of-sight velocity dispersions and derive dark matter halo virial masses for these host galaxies. At both epochs the velocity dispersion of satellites correlates with host galaxy stellar mass, \sigma\propto M_\ast^{0.4+/-0.1}, while the relation between satellite velocity dispersion and host galaxy B-band luminosity may grow somewhat shallower from \sigma\propto L_B^{0.6+/-0.1} at z~1 to \sigma\propto L_B^{0.4+/-0.1} at z~0. The evolution in M_200/M_\ast from z~1 to z~0 displays a bimodality insofar as host galaxies with stellar mass below M_\ast ~10^{11} M_Sun/h maintain a constant ratio (the intrinsic increase is constrained to a factor of 1.1+/-0.7) while host galaxies above M_\ast ~10^{11} M_Sun/h experience a factor of 4+/-3 increase in their virial-to-stellar mass ratio. This result can be easily understood if galaxies below this stellar mass scale continue to form stars while star formation in galaxies above this scale is quenched and the dark matter halos of galaxies both above and below this scale grow in accordance with LCDM cosmological simulations. Host galaxies that are red in U-B color have larger satellite dispersions and hence reside on average in more massive halos than blue galaxies at both z~1 and z~0. The redshift and host galaxy stellar mass dependence of M_200/M_\ast agrees qualitatively with the Millennium Run semi-analytic model of galaxy formation. (ABRIDGED)

The evolution of the B-band galaxy luminosity function is measured using a sample of more than 11,000 galaxies with spectroscopic redshifts from the DEEP2 Redshift Survey. The rest-frame M_B versus U-B color-magnitude diagram of DEEP2 galaxies shows that the color-magnitude bi-modality seen in galaxies locally is still present at redshifts z > 1. Dividing the sample at the trough of this color bimodality into predominantly red and blue galaxies, we find that the luminosity function of each galaxy color type evolves differently. Blue counts tend to shift to brighter magnitudes at constant number density, while the red counts remain largely constant at a fixed absolute magnitude. Using Schechter functions with fixed faint-end slopes we find that M*_B for blue galaxies brightens by ~ 1.3 magnitudes per unit redshift, with no significant evolution in number density. For red galaxies M*_B brightens somewhat less with redshift, while the formal value of phi* declines. When the population of blue galaxies is subdivided into two halves using the rest-frame color as the criterion, the measured evolution of both blue subpopulations is very similar.

We employ high-resolution dissipationless simulations of the concordance LCDM cosmology to model the observed luminosity dependence and evolution of galaxy clustering through most of the age of the universe, from z~5 to z~0. We use a simple, non-parametric model which monotonically relates galaxy luminosities to the maximum circular velocity of dark matter halos (V_max) by preserving the observed galaxy luminosity function in order to match the halos in simulations with observed galaxies. The novel feature of the model is the use of the maximum circular velocity at the time of accretion, V_max,acc, for subhalos, the halos located within virial regions of larger halos. We argue that for subhalos in dissipationless simulations, V_max,acc reflects the luminosity and stellar mass of the associated galaxies better than the circular velocity at the epoch of observation, V_max,now. The simulations and our model L-V_max relation predict the shape, amplitude, and luminosity dependence of the two-point correlation function in excellent agreement with the observed galaxy clustering in the SDSS data at z~0 and in the DEEP2 samples at z~1 over the entire probed range of projected separations, 0.1<r_p/(Mpc/h)<10.0. In particular, the small-scale upturn of the correlation function from the power-law form in the SDSS and DEEP2 luminosity-selected samples is reproduced very well. At z~3-5, our predictions also match the observed shape and amplitude of the angular two-point correlation function of Lyman-break galaxies (LBGs) on both large and small scales, including the small-scale upturn.

We present measurements of the void probability function (VPF) at z~1 using data from the DEEP2 Redshift Survey and its evolution to z~0 using data from the Sloan Digital Sky Survey (SDSS). We measure the VPF as a function of galaxy color and luminosity in both surveys and find that it mimics trends displayed in the two-point correlation function, $\xi$; namely that samples of brighter, red galaxies have larger voids (i.e. are more strongly clustered) than fainter, blue galaxies. We also clearly detect evolution in the VPF with cosmic time, with voids being larger in comoving units at z~0. We find that the reduced VPF matches the predictions of a `negative binomial' model for galaxies of all colors, luminosities, and redshifts studied. This model lacks a physical motivation, but produces a simple analytic prediction for sources of any number density and integrated two-point correlation function, \bar{\xi}. This implies that differences in the VPF across different galaxy populations are consistent with being due entirely to differences in the population number density and \bar{\xi}. The robust result that all galaxy populations follow the negative binomial model appears to be due to primarily to the clustering of dark matter halos. The reduced VPF is insensitive to changes in the parameters of the halo occupation distribution, in the sense that halo models with the same \bar{\xi} will produce the same VPF. For the wide range of galaxies studied, the VPF therefore does not appear to provide useful constraints on galaxy evolution models that cannot be gleaned from studies of \bar{\xi} alone. (abridged)

Using the first 25% of DEEP2 Redshift Survey data, we probe the line-of-sight velocity dispersion profile for isolated galaxies with absolute B-band magnitude -22<M_B-5log(h)<-21 at z=0.7-1.0, using satellite galaxies as luminous tracers of the underlying velocity distribution. Measuring the velocity dispersion beyond a galactocentric radius of ~200 kpc/h (physical) permits us to determine the total mass, including dark matter, around these bright galaxies. We find a line-of-sight velocity dispersion (sigma_los) of 162^{+44}_{-30} km/s at ~110 kpc/h, 136^{+26}_{-20} km/s at ~230 kpc/h, and 150^{+55}_{-38} km/s at ~320 kpc/h. Assuming an NFW model for the dark matter density profile, this corresponds to a mass within r_{200} of M_200=5.5^{+2.5}_{-2.0} x 10^12 M_Sun/h for our sample of satellite hosts with mean luminosity ~2.5L*. Roughly $~60% of these host galaxies have early-type spectra and are red in restframe (U-B) color, consistent with the overall DEEP2 sample in the same luminosity and redshift range. The halo mass determined for DEEP2 host galaxies is consistent with that measured in the Sloan Digital Sky Survey for host galaxies within a similar luminosity range relative to M*_B. This comparison is insensitive to the assumed halo mass profile, and implies an increase in the dynamical mass-to-light ratio (M_200/L_B) of isolated galaxies which host satellites by a factor of ~2.5 from z ~ 1 to z ~ 0. Our results are consistent with scenarios in which galaxies populate dark matter halos similarly from z ~ 0 to z ~ 1, except for ~1 magnitude of evolution in the luminosity of all galaxies.

The DEEP2 and COMBO-17 surveys are used to study the evolution of the luminosity function of red and blue galaxies to $z \sim 1$. Schechter function fits show that, since $z = 1$, $M^*_B$ dims by $\sim$ 1.3 mag per unit redshift for both color classes, $\phi^*$ of blue galaxies shows little change, while $\phi^*$ for red galaxies has formally nearly quadrupled. At face value, the number density of blue galaxies has remained roughly constant since $ z = 1$, whereas that of red galaxies has been rising. Luminosity densities support both conclusions, but we note that most red-galaxy evolution occurs between our data and local surveys and in our highest redshift bin, where the data are weakest. We discuss the implications of having most red galaxies emerge after $z = 1$ from precursors among the blue population, taking into account the properties of local and distant E/S0s. We suggest a ``mixed'' scenario in which some blue galaxies have their star-formation quenched in gas-rich mergers, migrate to the red sequence with a variety of masses, and merge further on the red sequence in one or more purely stellar mergers. E/S0s of a given mass today will have formed via different routes, in a manner that may help to explain the fundamental plane and other local scaling laws.

We use the first 25% of the DEEP2 Galaxy Redshift Survey spectroscopic data to identify groups and clusters of galaxies in redshift space. The data set contains 8370 galaxies with confirmed redshifts in the range 0.7 <= z <= 1.4, over one square degree on the sky. Groups are identified using an algorithm (the Voronoi-Delaunay Method) that has been shown to accurately reproduce the statistics of groups in simulated DEEP2-like samples. We optimize this algorithm for the DEEP2 survey by applying it to realistic mock galaxy catalogs and assessing the results using a stringent set of criteria for measuring group-finding success, which we develop and describe in detail here. We find in particular that the group-finder can successfully identify ~78% of real groups and that ~79% of the galaxies that are true members of groups can be identified as such. Conversely, we estimate that ~55% of the groups we find can be definitively identified with real groups and that ~46% of the galaxies we place into groups are interloper field galaxies. Most importantly, we find that it is possible to measure the distribution of groups in redshift and velocity dispersion, n(sigma,z), to an accuracy limited by cosmic variance, for dispersions greater than 350km/s. We anticipate that such measurements will allow strong constraints to be placed on the equation of state of the dark energy in the future. Finally, we present the first DEEP2 group catalog, which assigns 32% of the galaxies to 899 distinct groups with two or more members, 153 of which have velocity dispersions above 350km/s. We provide locations, redshifts and properties for this high-dispersion subsample. This catalog represents the largest sample to date of spectroscopically detected groups at z~1.