Measuring the chemical composition of galaxies is crucial to our understanding of galaxy formation and evolution models. However, such measurements are extremely challenging for quiescent galaxies at high redshifts, which have faint stellar continua and compact sizes, making it difficult to detect absorption lines and nearly impossible to spatially resolve them. Gravitational lensing offers the opportunity to study these galaxies with detailed spectroscopy that can be spatially resolved. In this work, we analyze deep spectra of MRG-M0138, a lensed quiescent galaxy at z = 1.98 which is the brightest of its kind, with an H-band magnitude of 17.1. Taking advantage of full spectral fitting, we measure $[{\rm Mg/Fe}]=0.51\pm0.05$, $[\rm{Fe/H}]=0.26\pm0.04$, and, for the first time, the stellar abundances of 6 other elements in this galaxy. We further constrained, also for the first time in a $z\sim2$ galaxy, radial gradients in stellar age, [Fe/H], and [Mg/Fe]. We detect no gradient in age or [Mg/Fe] and a slightly negative gradient in [Fe/H], which has a slope comparable to that seen in local early-type galaxies. Our measurements show that not only is MRG-M0138 very Mg-enhanced compared to the centers of local massive early-type galaxies, it is also very iron rich. These dissimilar abundances suggest that even the inner regions of massive galaxies have experienced significant mixing of stars in mergers, in contrast to a purely inside-out growth model. The abundance pattern observed in MRG-M0138 challenges simple galactic chemical evolution models that vary only the star formation timescale and shows the need for more elaborate models.

We introduce LATIS, the Ly$\alpha$ Tomography IMACS Survey, a spectroscopic survey at Magellan designed to map the z=2.2-2.8 intergalactic medium (IGM) in three dimensions by observing the Ly$\alpha$ forest in the spectra of galaxies and QSOs. Within an area of 1.7 deg${}^2$, we will observe approximately half of $\gtrsim L^*$ galaxies at z=2.2-3.2 for typically 12 hours, providing a dense network of sightlines piercing the IGM with an average transverse separation of 2.5 $h^{-1}$ comoving Mpc (1 physical Mpc). At these scales, the opacity of the IGM is expected to be closely related to the dark matter density, and LATIS will therefore map the density field in the $z \sim 2.5$ universe at $\sim$Mpc resolution over the largest volume to date. Ultimately LATIS will produce approximately 3800 spectra of z=2.2-3.2 galaxies that probe the IGM within a volume of $4 \times 10^6 h^{-3}$ Mpc${}^3$, large enough to contain a representative sample of structures from protoclusters to large voids. Observations are already complete over one-third of the survey area. In this paper, we describe the survey design and execution. We present the largest IGM tomographic maps at comparable resolution yet made. We show that the recovered matter overdensities are broadly consistent with cosmological expectations based on realistic mock surveys, that they correspond to galaxy overdensities, and that we can recover structures identified using other tracers. LATIS is conducted in Canada-France-Hawaii Telescope Legacy Survey fields, including COSMOS. Coupling the LATIS tomographic maps with the rich data sets collected in these fields will enable novel studies of environment-dependent galaxy evolution and the galaxy-IGM connection at cosmic noon.

Dark matter-only simulations predict that dark matter halos have steep, cuspy inner density profiles, while observations of dwarf galaxies find a range of inner slopes that are often much shallower. There is debate whether this discrepancy can be explained by baryonic feedback or if it may require modified dark matter models. In Paper 1 of this series, we obtained high-resolution integral field H$\alpha$ observations for 26 dwarf galaxies with $M_*=10^{8.1}-10^{9.7}\textrm{M}_\odot$. We derived rotation curves from our observations, which we use here to construct mass models. We model the total mass distribution as the sum of a generalized Navarro-Frenk-White (NFW) dark matter halo and the stellar and gaseous components. Our analysis of the slope of the dark matter density profile focuses on the inner 300-800 pc, chosen based on the resolution of our data and the region resolved by modern hydrodynamical simulations. The inner slope measured using ionized and molecular gas tracers is consistent, and it is additionally robust to the choice of stellar mass-to-light ratio. We find a range of dark matter profiles, including both cored and cuspy slopes, with an average of $\rho_{\rm DM}\sim r^{-0.74\pm 0.07}$, shallower than the NFW profile, but steeper than those typically observed for lower-mass galaxies with $M_*\sim 10^{7.5}\textrm{M}_\odot$. Simulations that reproduce the observed slopes in those lower-mass galaxies also produce slopes that are too shallow for galaxies in our mass range. We therefore conclude that supernova feedback models do not yet provide a fully satisfactory explanation for the observed trend in dark matter slopes.

We present the discovery and early evolution of ASASSN-19bt, a tidal disruption event (TDE) discovered by the All-Sky Automated Survey for Supernovae (ASAS-SN) at a distance of $d\simeq115$ Mpc and the first TDE to be detected by TESS. As the TDE is located in the TESS Continuous Viewing Zone, our dataset includes 30-minute cadence observations starting on 2018 July 25, and we precisely measure that the TDE begins to brighten $\sim8.3$ days before its discovery. Our dataset also includes 18 epochs of Swift UVOT and XRT observations, 2 epochs of XMM-Newton observations, 13 spectroscopic observations, and ground data from the Las Cumbres Observatory telescope network, spanning from 32 days before peak through 37 days after peak. ASASSN-19bt thus has the most detailed pre-peak dataset for any TDE. The TESS light curve indicates that the transient began to brighten on 2019 January 21.6 and that for the first 15 days its rise was consistent with a flux $\propto t^2$ power-law model. The optical/UV emission is well-fit by a blackbody SED, and ASASSN-19bt exhibits an early spike in its luminosity and temperature roughly 32 rest-frame days before peak and spanning up to 14 days that has not been seen in other TDEs, possibly because UV observations were not triggered early enough to detect it. It peaked on 2019 March 04.9 at a luminosity of $L\simeq1.3\times10^{44}$ ergs s$^{-1}$ and radiated $E\simeq3.2\times10^{50}$ ergs during the 41-day rise to peak. X-ray observations after peak indicate a softening of the hard X-ray emission prior to peak, reminiscent of the hard/soft states in X-ray binaries.

Dark-matter-only simulations predict that dark matter halos have cusp-like inner density profiles, while observations of low-mass galaxies have found a range of inner slopes that are typically much shallower. It is still not well established whether this discrepancy can be explained by baryonic feedback or if it may require modified dark matter models. To better understand the diversity of dark matter profiles in dwarf galaxies, we undertook a survey of 26 low-mass galaxies ($\log M_*/\textrm{M}_\odot = 8.4-9.8$, $v_{\rm max} = 50-140$ km s$^{-1}$) within 30 Mpc using the Palomar Cosmic Web Imager, which is among the largest integral field spectroscopic surveys of its type. In this paper, we derive H$\alpha$ velocity fields for the full sample with a typical spatial resolution of $\sim$160 pc. We extract rotation curves and verify their robustness to several choices in the analysis. We present a method for improving the velocity precision obtained from image slicing spectrographs using narrowband H$\alpha$ images. For 11 galaxies, we compare the H$\alpha$ velocity fields to CO kinematics measured using CARMA, finding the maps to be in good agreement. The standard deviation of the difference is typically $\sim$7 km s$^{-1}$, comparable to the level of turbulence in the interstellar medium, showing that the two tracers have substantially the same bulk kinematics. In a companion paper, we will use the rotation curves produced here to construct mass models of the galaxies and determine their dark matter density profiles.

We investigate the stellar populations for a sample of 24 quiescent galaxies at 1.5 < z < 2.5 using deep rest-frame optical spectra obtained with Keck MOSFIRE. By fitting templates simultaneously to the spectroscopic and photometric data, and exploring a variety of star formation histories, we obtain robust measurements of median stellar ages and residual levels of star formation. After subtracting the stellar templates, the stacked spectrum reveals the Halpha and [NII] emission lines, providing an upper limit on the ongoing star formation rate of 0.9 +/- 0.1 Msun/yr. By combining the MOSFIRE data to our sample of Keck LRIS spectra at lower redshift, we analyze in a consistent manner the quiescent population at 1 < z < 2.5. We find a tight relation (with a scatter of 0.13 dex) between the stellar age and the rest-frame U-V and V-J colors, which can be used to estimate the age of quiescent galaxies given their colors. Applying this age--color relation to large, photometric samples, we are able to model the number density evolution for quiescent galaxies of various ages. We find evidence for two distinct quenching paths: a fast quenching that produces compact post-starburst systems, and a slow quenching of larger galaxies. Fast quenching accounts for about a fifth of the growth of the red sequence at z~1.4, and half at z~2.2. We conclude that fast quenching is triggered by dramatic events such as gas-rich mergers, while slow quenching is likely caused by a different physical mechanism.

Stellar kinematics provide insights into the masses and formation histories of galaxies. At high redshifts, spatially resolving the stellar kinematics of quiescent galaxies is challenging due to their compact sizes. Using deep near-infrared spectroscopy, we have measured the resolved stellar kinematics of four quiescent galaxies at z=1.95-2.64, introduced in Paper I, that are gravitationally lensed by galaxy clusters. Analyses of two of these have previously been reported individually by Newman et al. and Toft et al., and for the latter we present new observations. All four galaxies show significant rotation and can be classified as "fast rotators." In the three systems for which the lensing constraints permit a reconstruction of the source, we find that all are likely to be highly flattened (intrinsic ellipticities of $\approx0.75-0.85$) disk-dominated galaxies with rapid rotation speeds of $V_{\rm max}=290-352$ km/s and predominantly rotational support, as indicated by the ratio $(V/\sigma)_{R_e}=1.7-2.3$. Compared to coeval star-forming galaxies of similar mass, the quiescent galaxies have smaller $V/\sigma$. Given their high masses $M_{\rm dyn} \gtrsim 2\times10^{11} M_{\odot}$, we argue that these galaxies are likely to evolve into "slow rotator" elliptical galaxies whose specific angular momentum is reduced by a factor of 5-10. This provides strong evidence for merger-driven evolution of massive galaxies after quenching. Consistent with indirect evidence from earlier morphological studies, our small but unique sample suggests that the kinematic transformations that produced round, dispersion-supported elliptical galaxies were not generally coincident with quenching. Such galaxies probably emerged later via mergers that increased their masses and sizes while also eroding their rotational support.

Quiescent galaxies at $z \gtrsim 2$ are compact and have weak or absent emission lines, making it difficult to spatially resolve their kinematics and stellar populations using ground-based spectroscopy. Gravitationally lensed examples provide a promising route forward, but such objects are very rare. We describe a search in the fields of 232 galaxy clusters that has uncovered five bright ($H_{\rm AB} < 20$) lensed galaxies with red near-infrared colors. These include MRG-M0138, the brightest lensed galaxy known in the near-infrared. Analysis of near-infrared spectra and multiband photometry confirms that all are quiescent galaxies at z=1.95-2.64 with stellar ages of 0.5-1.4 Gyr (corresponding to formation epochs $z_{\rm form} \simeq 3-4$) and stellar masses of $10^{11.6-12.8} \mu^{-1} M_{\odot}$, where $\mu$ is the magnification. In three cases we derive lens models and reconstruct the source structure; these galaxies are massive ($M_* \gtrsim 10^{11.0} M_{\odot}$) and follow the mass--size relation defined by unlensed samples. In two of these three galaxies, the main structural component is an inclined disk. Weak emission lines are detected in four of five galaxies with high ratios [NII]/H$\alpha$~2-6 that are inconsistent with a star formation origin. Based on the line ratios, the H$\alpha$ equivalent widths, and the distribution and kinematics of the gas, we infer that shocks are likely present in at least two galaxies and could be present in all of the line emitters. We speculate that these could be analogs of local galaxies in which AGN jet-driven outflows are thought to heat the interstellar medium and suppress star formation. In further papers we will present spatially resolved measurements of the stellar populations and kinematics of this unique sample.

Modern data empower observers to describe galaxies as the spatially and biographically complex objects they are. We illustrate this through case studies of four, $z\sim1.3$ systems based on deep, spatially resolved, 17-band + G102 + G141 Hubble Space Telescope grism spectrophotometry. Using full spectrum rest-UV/-optical continuum fitting, we characterize these galaxies' observed $\sim$kpc-scale structures and star formation rates (SFRs) and reconstruct their history over the age of the universe. The sample's diversity---passive to vigorously starforming; stellar masses $\log M_*/M_\odot=10.5$ to $11.2$---enables us to draw spatio-temporal inferences relevant to key areas of parameter space (Milky Way- to super-Andromeda-mass progenitors). Specifically, we find signs that bulge mass-fractions ($B/T$) and SF history shapes/spatial uniformity are linked, such that higher $B/T$s correlate with "inside-out growth" and central specific SFRs that peaked above the global average for all starforming galaxies at that epoch. Conversely, the system with the lowest $B/T$ had a flat, spatially uniform SFH with normal peak activity. Both findings are consistent with models positing a feedback-driven connection between bulge formation and the switch from rising to falling SFRs ("quenching"). While sample size forces this conclusion to remain tentative, this work provides a proof-of-concept for future efforts to refine or refute it: JWST, WFIRST, and the 30-m class telescopes will routinely produce data amenable to this and more sophisticated analyses. These samples---spanning representative mass, redshift, SFR, and environmental regimes---will be ripe for converting into thousands of sub-galactic-scale empirical windows on what individual systems actually looked like in the past, ushering in a new dialog between observation and theory.

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 report the discovery of a massive $\log(M/M_\odot)=10.74^{+0.18}_{-0.16}$ galaxy at the same redshift as a carbon-monoxide-bearing sub-damped Lyman $\alpha$ absorber (sub-DLA) seen in the spectrum of QSO J1439+1117. The galaxy, J1439B, is located 4\farcs7 from the QSO sightline, a projected distance of 38 physical kpc at $z=2.4189$, and exhibits broad optical emission lines ($\sigma_{\rm{[O III]}}=303 \pm 12$~\kms) with ratios characteristic of excitation by an active galactic nucleus (AGN). The galaxy has a factor of $\sim$9 lower star formation than is typical of star-forming galaxies of the same mass and redshift. The nearby sub-DLA is highly enriched, suggesting its galactic counterpart must be massive if it follows the $z\sim2$ mass-metallicity relationship. Metallic absorption within the circumgalactic medium of the sub-DLA and J1439B is spread over a velocity range $\Delta v > 1000$ \kms, suggesting an energetic origin. We explore the possibility that a different galaxy could be responsible for the rare absorber, and conclude that it is unlikely based on imaging, integral-field spectroscopy, and high-$z$ massive galaxy pair statistics. We argue that the gas seen in absorption against the QSO was likely ejected from the galaxy J1439B and therefore provides a unique observational probe of AGN feedback in the distant universe.

We present deep near-infrared spectra for a sample of 24 quiescent galaxies in the redshift range 1.5 < z < 2.5 obtained with the MOSFIRE spectrograph at the W. M. Keck Observatory. In conjunction with a similar dataset we obtained in the range 1 < z < 1.5 with the LRIS spectrograph, we analyze the kinematic and structural properties for 80 quiescent galaxies, the largest homogeneously-selected sample to date spanning 3 Gyr of early cosmic history. Analysis of our Keck spectra together with measurements derived from associated HST images reveals increasingly larger stellar velocity dispersions and smaller sizes to redshifts beyond z~2. By classifying our sample according to Sersic indices, we find that among disk-like systems the flatter ones show a higher dynamical to stellar mass ratio compared to their rounder counterparts which we interpret as evidence for a significant contribution of rotational motion. For this subset of disk-like systems, we estimate that V/sigma, the ratio of the circular velocity to the intrinsic velocity dispersion, is a factor of two larger than for present-day disky quiescent galaxies. We use the velocity dispersion measurements also to explore the redshift evolution of the dynamical to stellar mass ratio, and to measure for the first time the physical size growth rate of individual systems over two distinct redshift ranges, finding a faster evolution at earlier times. We discuss the physical origin of this time-dependent growth in size in the context of the associated reduction of the systematic rotation.

We present the Magellan/FIRE detection of highly-ionized CIV 1550 and OIII] 1666 in a deep infrared spectrum of the z=6.11 gravitationally lensed low-mass galaxy RXC J2248.7-4431-ID3, which has previously-known Lyman-alpha. No corresponding emission is detected at the expected location of HeII 1640. The upper limit on HeII paired with detection of OIII] and CIV constrains possible ionization scenarios. Production of CIV and OIII] requires ionizing photons of 2.5-3.5 Ryd, but once in that state their multiplet emission is powered by collisional excitation at lower energies (~0.5 Ryd). As a pure recombination line, HeII emission is powered by 4 Ryd ionizing photons. The data therefore require a spectrum with significant power at 3.5 Ryd but a rapid drop toward 4.0 Ryd. This hard spectrum with a steep drop is characteristic of low-metallicity stellar populations, and less consistent with soft AGN excitation, which features more 4 Ryd photons and hence higher HeII flux. The conclusions based on ratios of metal line detections to Helium non-detection are strengthened if the gas metallicity is low. RXJ2248-ID3 adds to the growing handful of reionization-era galaxies with UV emission line ratios distinct from the general z=2-3 population, in a way that suggests hard ionizing spectra that do not necessarily originate in AGN.

Observations of strong gravitational lensing, stellar kinematics, and larger-scale tracers enable accurate measures of the distribution of dark matter (DM) and baryons in massive early-type galaxies (ETGs). While such techniques have been applied to galaxy-scale and cluster-scale lenses, the paucity of intermediate-mass systems with high-quality data has precluded a uniform analysis of mass-dependent trends. With the aim of bridging this gap, we present new observations and analyses of 10 group-scale lenses at <z>=0.36 characterized by Einstein radii theta_Ein=2.5"-5.1" and a mean halo mass of M_200=10^14.0 Msol. We measure a mean concentration c_200=5.0+-0.8 consistent with unmodified cold dark matter halos. By combining our data with other lens samples, we analyze the mass structure of ETGs in 10^13 Msol-10^15 Msol halos using homogeneous techniques. We show that the slope of the total density profile gamma_tot within the effective radius depends on the stellar surface density, as demonstrated previously, but also on the halo mass. We analyze these trends using halo occupation models and resolved stellar kinematics with the goal of testing the universality of the DM profile. Whereas the central galaxies of clusters require a shallow inner DM density profile, group-scale lenses are consistent with a Navarro-Frenk-White profile or one that is slightly contracted. The largest uncertainties arise from the sample size and likely radial gradients in stellar populations. We conclude that the net effect of baryons on the DM distribution may not be universal, but more likely varies with halo mass due to underlying trends in star formation efficiency and assembly history.

We report the discovery of RG1M0150, a massive, recently quenched galaxy at z=2.636 that is multiply imaged by the cluster MACSJ0150.3-1005. We derive a stellar mass of log M_*=11.49+0.10-0.16 and a half-light radius of R_e,maj =1.8+-0.4 kpc. Taking advantage of the lensing magnification, we are able to spatially resolve a remarkably massive yet compact quiescent galaxy at z>2 in ground-based near-infrared spectroscopic observations using Magellan/FIRE and Keck/MOSFIRE. We find no gradient in the strength of the Balmer absorption lines over 0.6 R_e - 1.6 R_e, which are consistent with an age of 760 Myr. Gas emission in [NII] broadly traces the spatial distribution of the stars and is coupled with weak Halpha emission (log [NII]/Halpha = 0.6+-0.2), indicating that OB stars are not the primary ionizing source. The velocity dispersion within the effective radius is sigma_e = 271+-41 km/s. We detect rotation in the stellar absorption lines for the first time beyond z~1. Using a two-integral Jeans model that accounts for observational effects, we measure a dynamical mass of log M_dyn =11.24+-0.14 and V/sigma=0.70+-0.21. This is a high degree of rotation considering the modest observed ellipticity of 0.12+-0.08, but it is consistent with predictions from dissipational merger simulations that produce compact remnants. The mass of RG1M0150 implies that it is likely to become a slowly rotating elliptical. If it is typical, this suggests that the progenitors of massive ellipticals retain significant net angular momentum after quenching which later declines, perhaps through accretion of satellites.

We analyze the stellar populations of a sample of 62 massive (log Mstar/Msun > 10.7) galaxies in the redshift range 1 < z < 1.6, with the main goal of investigating the role of recent quenching in the size growth of quiescent galaxies. We demonstrate that our sample is not biased toward bright, compact, or young galaxies, and thus is representative of the overall quiescent population. Our high signal-to-noise ratio Keck LRIS spectra probe the rest-frame Balmer break region which contains important absorption line diagnostics of recent star formation activity. We obtain improved measures of the various stellar population parameters, including the star-formation timescale tau, age and dust extinction, by fitting templates jointly to both our spectroscopic and broad-band photometric data. We identify which quiescent galaxies were recently quenched and backtrack their individual evolving trajectories on the UVJ color-color plane finding evidence for two distinct quenching routes. By using sizes measured in the previous paper of this series, we confirm that the largest galaxies are indeed among the youngest at a given redshift. This is consistent with some contribution to the apparent growth from recent arrivals, an effect often called progenitor bias. However, we calculate that recently-quenched objects can only be responsible for about half the increase in average size of quiescent galaxies over a 1.5 Gyr period, corresponding to the redshift interval 1.25 < z < 2. The remainder of the observed size evolution arises from a genuine growth of long-standing quiescent galaxies.

Using the MOSFIRE near-infrared multi-slit spectrograph on the Keck 1 Telescope, we have secured high signal-to-noise ratio absorption line spectra for six massive galaxies with redshift 2 < z < 2.5. Five of these galaxies lie on the red sequence and show signatures of passive stellar populations in their rest-frame optical spectra. By fitting broadened spectral templates we have determined stellar velocity dispersions and, with broad-band HST and Spitzer photometry and imaging, stellar masses and effective radii. Using this enlarged sample of galaxies we confirm earlier suggestions that quiescent galaxies at z > 2 have small sizes and large velocity dispersions compared to local galaxies of similar stellar mass. The dynamical masses are in very good agreement with stellar masses (log Mstar/Mdyn = -0.02 +/- 0.03), although the average stellar-to-dynamical mass ratio is larger than that found at lower redshift (-0.23 +/- 0.05). By assuming evolution at fixed velocity dispersion, not only do we confirm a surprisingly rapid rate of size growth but we also consider the necessary evolutionary track on the mass-size plane and find a slope alpha = dlogR / dlogM > ~2 inconsistent with most numerical simulations of minor mergers. Both results suggest an additional mechanism may be required to explain the size growth of early galaxies.

We present Hubble Space Telescope imaging and grism spectroscopy in the field of the distant galaxy cluster JKCS 041 using the Wide Field Camera 3. We confirm that JKCS 041 is a rich cluster and derive a redshift z=1.80 via the spectroscopic identification of 19 member galaxies, of which 15 are quiescent. These are centered upon diffuse X-ray emission seen by the Chandra observatory. As JKCS 041 is the most distant known cluster with such a large, spectroscopically confirmed quiescent population, it provides a unique opportunity to study the effect of the environment on galaxy properties at early epochs. We construct high-quality composite spectra of the quiescent cluster members that reveal prominent Balmer and metallic absorption lines. Using these, we measure the mean stellar ages in two bins of stellar mass. The quiescent cluster members' ages agree remarkably closely with that inferred by Whitaker et al. for similarly selected samples in the field, supporting the idea that the cluster environment is more efficient at truncating star formation while not having a strong effect on the mean epoch of quenching. We find some evidence (90% confidence) for a lower fraction of disk-like quiescent systems in JKCS 041 compared to a sample of coeval field galaxies drawn from the CANDELS survey. Taking this into account, we do not detect a significant difference between the mass-radius relations of the quiescent JKCS 041 members and our z~1.8 field sample. Finally, we demonstrate how differences in the morphological mixture of quenched systems can complicate measures of the environmental dependence of size growth.

Using deep two-color near-infrared HST imaging and unbiased grism spectroscopy we present a detailed study of the z=1.803 JKCS041 cluster. Uniquely, for a high redshift cluster, we confirm a mass of $\log M=14.2$ in solar units using four different techniques based on the X-ray temperature, the X-ray luminosity and the cluster richness. JKCS041 is thus a progenitor of a local system like the Coma cluster. Our rich dataset and the abundant population of 14 spectroscopically-confirmed red sequence galaxies allows us to explore the past star formation history of this system in unprecedented detail. Remarkably, we find a prominent red sequence down to stellar masses as low as $\log M=9.8$, corresponding to a mass range of 2 dex. These quiescent galaxies are concentrated around the cluster center with a core radius of 330 kpc. Blue members are few and avoid the cluster center. In JKCS041 quenching was therefore largely completed by a look-back time of 10 Gyr and we can constrain the epoch at which this occurred via spectroscopic age-dating of the individual galaxies. Most galaxies were quenched about 1.1 Gyr prior to the epoch of observation. The less massive quiescent galaxies are somewhat younger, corresponding to a decrease in age of 650 Myr per mass dex, but the scatter in age at fixed mass is only 380 Myr (at $\log M=11$). The size-mass relation of quiescent galaxies in JKCS041 is consistent with that observed for local clusters within our uncertainties. Comparing our data on JKCS041 with 41 clusters at lower redshift, we find that the form of the mass function of red sequence galaxies has hardly evolved in the past 10 Gyr, both in terms of its faint end slope and characteristic mass. Despite observing JKCS041 soon after its quenching and the three fold expected increase in mass in the next 10 Gyr, it is already remarkably similar to present-day clusters.

We present Keck LRIS spectroscopy for a sample of 103 massive galaxies with redshifts 0.9 < z < 1.6. Of these, 56 are quiescent with high signal-to-noise absorption line spectra, enabling us to determine robust stellar velocity dispersions for the largest sample yet available beyond a redshift of 1. Together with effective radii measured from deep Hubble Space Telescope images, we calculate dynamical masses and address key questions relating to the puzzling size growth of quiescent galaxies over 0 < z < 2. We examine the relationship between stellar and dynamical masses at high redshift, finding that it closely follows that determined locally. We also confirm the utility of the locally-established empirical calibration which enables high-redshift velocity dispersions to be estimated photometrically, and we determine its accuracy to be 35%. To address recent suggestions that progenitor bias - the continued arrival of recently-quenched larger galaxies - can largely explain the size evolution of quiescent galaxies, we examine the growth at fixed velocity dispersion assuming this quantity is largely unaffected by the merger history. We demonstrate that significant size and mass growth have clearly occurred in individual systems. Parameterizing the relation between mass and size growth over 0 < z < 1.6 as R \propto M^alpha, we find alpha = 1.6 +- 0.3, in agreement with theoretical expectations from simulations of minor mergers. Relaxing the assumption that the velocity dispersion is unchanging, we examine growth assuming a constant ranking in galaxy velocity dispersion. This approach is applicable only to the large-dispersion tail of the distribution, but yields a consistent growth rate of alpha = 1.4 +- 0.2. Both methods confirm that progenitor bias alone is insufficient to explain our new observations and that quiescent galaxies have grown in both size and stellar mass over 0 < z < 1.6.

A number of recent challenges to the standard Lambda-CDM paradigm relate to discrepancies that arise in comparing the abundance and kinematics of local dwarf galaxies with the predictions of numerical simulations. Such arguments rely heavily on the assumption that the local dwarf and satellite galaxies form a representative distribution in terms of their stellar-to-halo mass ratios. To address this question, we present new, deep spectroscopy using DEIMOS on Keck for 82 low mass (10^7-10^9 solar masses) star-forming galaxies at intermediate redshift (z=0.2-1). For 50 percent of these we are able to determine resolved rotation curves using nebular emission lines and thereby construct the stellar mass Tully-Fisher relation to masses as low as 10^7 solar masses. Using scaling relations determined from weak lensing data, we convert this to a stellar-to-halo mass (SHM) relation for comparison with abundance matching predictions. We find a discrepancy between the propagated predictions from simulations compared to our observations, and suggest possible reasons for this as well as future tests that will be more effective.

We provide a new observational test for a key prediction of the \Lambda CDM cosmological model: the contributions of mergers with different halo-to-main-cluster mass ratios to cluster-sized halo growth. We perform this test by dynamically analyzing seven galaxy clusters, spanning the redshift range $0.13 < z_c < 0.45$ and caustic mass range $0.4-1.5$ $10^{15} h_{0.73}^{-1}$ M$_{\odot}$, with an average of 293 spectroscopically-confirmed bound galaxies to each cluster. The large radial coverage (a few virial radii), which covers the whole infall region, with a high number of spectroscopically identified galaxies enables this new study. For each cluster, we identify bound galaxies. Out of these galaxies, we identify infalling and accreted halos and estimate their masses and their dynamical states. Using the estimated masses, we derive the contribution of different mass ratios to cluster-sized halo growth. For mass ratios between ~0.2 and ~0.7, we find a ~1 $\sigma$ agreement with \Lambda CDM expectations based on the Millennium simulations I and II. At low mass ratios, $\lesssim 0.2$, our derived contribution is underestimated since the detection efficiency decreases at low masses, $\sim 2 \times 10^{14}$ $h_{0.73}^{-1}$ M$_{\odot}$. At large mass ratios, $\gtrsim 0.7$, we do not detect halos probably because our sample, which was chosen to be quite X-ray relaxed, is biased against large mass ratios. Therefore, at large mass ratios, the derived contribution is also underestimated.

[abridged] The mass-size relation of early-type galaxies (ETGs) has been largely studied in the last years to probe the mass assembly of the most massive objects in the Universe. In this paper, we focus on the mass-size relation of quiescent massive ETGs (Mstar/Msol > 3*10^10) living in massive clusters (M200 ~ 10^14 Mstar) at 0.8< z <1.5, as compared to those living in the field at the same epoch. Our sample contains ~ 400 ETGs in clusters and the same number in the field. Therefore, our sample is approximately an order of magnitude larger than previous studies in the same redshift range for galaxy clusters. We find that ETGs living in clusters are between ~30-50% larger than galaxies with the same stellar mass residing in the field. We parametrize the size using the mass-normalized size, gamma=Re/Mstar^0.57. The gamma distributions in both environments peak at the same position but the distributions in clusters are more skewed towards larger sizes. Since this size difference is not observed in the local Universe, the size evolution at fixed stellar mass from z~1.5 of cluster galaxies is less steep ((1+z)-0.53pm0.04) than the evolution of field galaxies ((1+z)-0.92pm0.04). The size difference seems to be essentially driven by the galaxies residing in the clusters cores (R<0.5*R200). If part of the size evolution is due to mergers, the difference we see between cluster and field galaxies could be due to higher merger rates in clusters at higher redshift, probably during the formation phase of the clusters when velocity dispersions are lower. We cannot exclude however that the difference is driven by newly quenched galaxies which are quenched more efficiently in clusters. The implications of these results for the hierarchical growth of ETGs will be discussed in a companion paper.

Clusters of galaxies are excellent locations to probe the distribution of baryons and dark matter (DM) over a wide range of scales. We study a sample of seven massive, relaxed galaxy clusters with centrally-located brightest cluster galaxies (BCGs) at z=0.2-0.3. Using the observational tools of strong and weak gravitational lensing, combined with resolved stellar kinematics within the BCG, we measure the total radial density profile, comprising both dark and baryonic matter, over scales of ~3-3000 kpc. Lensing-derived mass profiles typically agree with independent X-ray estimates within ~15%, suggesting that departures from hydrostatic equilibrium are small and that the clusters in our sample (except A383) are not strongly elongated along the line of sight. The inner logarithmic slope gamma_tot of the total density profile measured over r/r200=0.003-0.03, where rho_tot ~ r^(-gamma_tot), is found to be nearly universal, with a mean <gamma_tot> = 1.16 +- 0.05 (random) +0.05-0.07 (systematic) and an intrinsic scatter of < 0.13 (68% confidence). This is further supported by the very homogeneous shape of the observed velocity dispersion profiles, obtained via Keck spectroscopy, which are mutually consistent after a simple scaling. Remarkably, this slope agrees closely with numerical simulations that contain only dark matter, despite the significant contribution of stellar mass on the scales we probe. The Navarro-Frenk-White profile characteristic of collisionless cold dark matter is a better description of the total mass density at radii >~ 5-10 kpc than that of dark matter alone. Hydrodynamical simulations that include baryons, cooling, and feedback currently provide a poorer match. We discuss the significance of our findings for understanding the assembly of BCGs and cluster cores, particularly the influence of baryons on the inner DM halo. [abridged]

We present stellar and dark matter (DM) density profiles for a sample of seven massive, relaxed galaxy clusters derived from strong and weak gravitational lensing and resolved stellar kinematic observations within the centrally-located brightest cluster galaxies (BCGs). In Paper I of the series, we demonstrated that the total density profile derived from these data, which span 3 decades in radius, is consistent with numerical DM-only simulations at radii >~ 5-10 kpc, despite the significant contribution of stellar material in the core. Here we decompose the inner mass profiles of these clusters into stellar and dark components. Parametrizing the DM density profile as a power law rho_DM ~ r^{-\beta} on small scales, we find a mean slope <\beta> = 0.50 +- 0.10 (random) +0.14-0.13 (systematic). Alternatively, cored Navarro-Frenk-White (NFW) profiles with <log r_core/kpc> = 1.14 +- 0.13 (random) +0.14-0.22 (systematic) provide an equally good description. These density profiles are significantly shallower than canonical NFW models at radii <~ 30 kpc, comparable to the effective radii of the BCGs. The inner DM profile is correlated with the distribution of stars in the BCG, suggesting a connection between the inner halo and the assembly of stars in the central galaxy. The stellar mass-to-light ratio inferred from lensing and stellar dynamics is consistent with that inferred using stellar population synthesis models if a Salpeter initial mass function is adopted. We compare these results to theories describing the interaction between baryons and DM in cluster cores, including adiabatic contraction models and the possible effects of galaxy mergers and active galactic nucleus feedback, and evaluate possible signatures of alternative DM candidates.

Through extended integrations using the recently-installed deep depletion CCD on the red arm of the Keck I Low Resolution Imaging Spectrograph, we present new measurements of the resolved spectra of 70 morphologically-selected star-forming galaxies with i_AB<24.1 in the redshift range 1<z<1.7. Using the formalism introduced in Paper I of this series and available HST ACS images, we successfully recover rotation curves using the extended emission line distribution of [O II] 3727 A to 2.2 times the disk scale radius for a sample of 42 galaxies. Combining these measures with stellar masses derived from HST and ground-based near-infrared photometry enables us to construct the stellar mass Tully-Fisher relation in the time interval between the well-constructed relation defined at z~1 in Paper I and the growing body of resolved dynamics probed with integral field unit spectrographs at z>2. Remarkably, we find a well-defined Tully-Fisher relation with up to 60% increase in scatter and stellar mass zero-point shift constraint of 0.02+/-0.02 dex since z~1.7, compared to the local relation. Although our sample is incomplete in terms of either a fixed stellar mass or star formation rate limit, we discuss the implications that typical star-forming disk galaxies evolve to arrive on a well-defined Tully-Fisher relation within a surprisingly short period of cosmic history.

Early-type galaxies (ETGs) are observed to be more compact at z>2 than in the local Universe. Remarkably, much of this size evolution appears to take place in a short (1.8 Gyr) time span between z=2.2 and z=1.3, which poses a serious challenge to hierarchical galaxy formation models where mergers occurring on a similar timescale are the main mechanism for galaxy growth. We compute the merger-driven redshift evolution of stellar mass Mstar\propto(1+z)^aM, half-mass radius Re\propto(1+z)^aR and velocity-dispersion sigma0\propto(1+z)^asigma predicted by concordance Lambda cold dark matter for a typical massive ETG in the redshift range z=1.3-2.2. Neglecting dissipative processes, and thus maximizing evolution in surface density, we find -1.5<aM<-0.6, -1.9<aR<-0.7 and 0.06<asigma<0.22, under the assumption that the accreted satellites are spheroids. It follows that the predicted z=2.2 progenitors of z=1.3 ETGs are significantly less compact (on average a factor of 2 larger Re at given Mstar) than the quiescent galaxies observed at z>2. Furthermore, we find that the scatter introduced in the size-mass correlation by the predicted merger-driven growth is difficult to reconcile with the tightness of the observed scaling law. We conclude that - barring unknown systematics or selection biases in the current measurements - minor and major mergers with spheroids are not sufficient to explain the observed size growth of ETGs within the standard model.

The presence of extremely compact galaxies at z~2 and their subsequent growth in physical size has been the cause of much puzzlement. We revisit the question using deep infrared Wide Field Camera 3 data to probe the rest-frame optical structure of 935 host galaxies selected with 0.4<z<2.5 and stellar masses M* > 10^10.7 Msol using optical and near-infrared photometry in the UKIRT Ultra Deep Survey and GOODS-South fields of the CANDELS survey. At each redshift, the most compact sources are those with little or no star formation, and we find that the mean size of these systems grows by a factor of 3.5 +- 0.3 over this redshift interval. The new data are sufficiently deep to enable us to identify companions to these hosts whose stellar masses are ten times smaller, while still yielding suitably accurate photometric redshifts to define a likely physical association. By searching for faint companions around 404 quiescent hosts within a projected physical annulus 10 < R < 30 kpc/h, we estimate the minor merger rate over the redshift range 0.4 < z < 2. After correcting for contamination from projected pairs, we find that 13-18% of quiescent hosts have likely physical companions with stellar mass ratios of 0.1 or greater. Mergers of these companions will typically increase the host mass by 6+-2% per merger timescale. We estimate the minimum growth rate necessary to explain the declining abundance of compact galaxies. Using a simple model of merging motivated by recent numerical simulations, we then assess whether mergers of the faint companions with their hosts are sufficient to explain this minimal rate. We find that mergers with mass ratios > 0.1 may explain most of the size evolution observed at z >~ 1 if a relatively short merger timescale is assumed, but the rapid growth seen at higher redshift likely requires additional physical processes.

We extend our analyses of the dark matter (DM) distribution in relaxed clusters to the case of Abell 383, a luminous X-ray cluster at z=0.189 with a dominant central galaxy and numerous strongly-lensed features. Following our earlier papers, we combine strong and weak lensing constraints secured with Hubble Space Telescope and Subaru imaging with the radial profile of the stellar velocity dispersion of the central galaxy, essential for separating the baryonic mass distribution in the cluster core. Hydrostatic mass estimates from Chandra X-ray observations further constrain the solution. These combined datasets provide nearly continuous constraints extending from 2 kpc to 1.5 Mpc in radius, allowing stringent tests of results from recent numerical simulations. Two key improvements in our data and its analysis make this the most robust case yet for a shallow slope \beta of the DM density profile \rho_DM ~ r^-\beta on small scales. First, following deep Keck spectroscopy, we have secured the stellar velocity dispersion profile to a radius of 26 kpc for the first time in a lensing cluster. Secondly, we improve our previous analysis by adopting a triaxial DM distribution and axisymmetric dynamical models. We demonstate that in this remarkably well-constrained system, the logarithmic slope of the DM density at small radii is \beta < 1.0 (95% confidence). An improved treatment of baryonic physics is necessary, but possibly insufficient, to reconcile our observations with the recent results of high-resolution simulations.

We present deep Keck spectroscopy for 17 morphologically-selected field spheroidals in the redshift range 1.05<z<1.60 in order to investigate the continuity in physical properties between the claimed massive compact red galaxies ("nuggets") at z~2 and well-established data for massive spheroidal galaxies below z~1. By combining Keck-based stellar velocity dispersions with HST-based sizes, we find that the most massive systems (Mdyn > 10^11 Msol) grew in size over 0<z<1.6 as (1+z)^(-0.75 +- 0.10) (i.e., x2 since z=1.5) whereas intermediate mass systems (10^11 Msol > Mdyn > 10^10 Msol) did not grow significantly. These trends are consistent with a picture in which more massive spheroidals formed at higher redshift via "wetter" mergers involving greater dissipation. To examine growth under the favored "dry" merger hypothesis, we also examine size growth at a fixed velocity dispersion. This test, uniquely possible with our dynamical data, allows us to consider the effects of "progenitor bias." Above our completeness limit (sigma > 200 km/s), we find size growth consistent with that inferred for the mass-selected sample, thus ruling out strong progenitor bias. To maintain continuity in the growth of massive galaxies over the past 10 Gyr, our new results imply that size evolution over 1.3<z<2.3, a period of 1.9 Gyr, must have been even more dramatic than hitherto claimed if the red sources at z>2 are truly massive and compact.

We present a detailed analysis of the baryonic and dark matter distribution in the lensing cluster Abell 611 (z=0.288), with the goal of determining the dark matter profile over an unprecedented range of cluster-centric distance. By combining three complementary probes of the mass distribution, weak lensing from deep multi-color imaging, strong lensing constraints based on the identification of multiply-imaged sources, and resolved stellar velocity dispersion measures for the brightest cluster galaxy (BCG), we extend the methodology for separating the dark and baryonic mass components introduced by Sand et al. (2008). Our resulting dark matter profile samples the cluster from ~3 kpc to 3.25 Mpc, thereby providing an excellent basis for comparisons with recent numerical models. We demonstrate that only by combining our three observational techniques can degeneracies in constraining the form of the dark matter profile be broken on scales crucial for detailed comparisons with numerical simulations. Our analysis reveals that a simple Navarro, Frenk, and White (NFW) profile is an unacceptable fit to our data. We confirm earlier claims that the inner profile of the dark matter profile deviates significantly from the NFW form and find a inner logarithmic slope \beta flatter than 0.3 (68%; where rho_DM ~ r^{-\beta} at small radii). In order to reconcile our data with cluster formation in a LambdaCDM cosmology, we speculate that it may be necessary to revise our understanding of the nature of baryon--dark matter interactions in cluster cores. Comprehensive weak and strong lensing data, when coupled with kinematic information on the brightest cluster galaxy, can readily be applied to a larger sample of clusters to test the universality of these results.

We report the successful identification of the type of the supernova responsible for the supernova remnant SNR 0509-675 in the Large Magellanic Cloud (LMC) using Gemini spectra of surrounding light echoes. The ability to classify outbursts associated with centuries-old remnants provides a new window into several aspects of supernova research and is likely to be successful in providing new constraints on additional LMC supernovae as well as their historical counterparts in the Milky Way Galaxy (MWG). The combined spectrum of echo light from SNR 0509-675 shows broad emission and absorption lines consistent with a supernova (SN) spectrum. We create a spectral library consisting of 26 SNe Ia and 6 SN Ib/c that are time-integrated, dust-scattered by LMC dust, and reddened by the LMC and MWG. We fit these SN templates to the observed light echo spectrum using $\chi^2$ minimization as well as correlation techniques, and we find that overluminous 91T-like SNe Ia with $\dm15<0.9$ match the observed spectrum best.

Light echoes are interesting because of the wealth of information they offer about their progenitors and the reflecting dust. Due to their faint surface brightnesses, difference imaging is necessary to separate most light echoes from the sky background. However, difference images reveal only the relative fluxes between two epochs. Obtaining absolute fluxes for individual epochs has traditionally relied on a single template image that is free of light echoes. Since such an observation is normally unavailable, a light echo-free template must be constructed by a complicated and usually subjective process. Here we present an application of the NN2 method of Barris et al. to extract the relative fluxes of light echoes across a range of epochs directly from a series of difference images. This method requires no privileged image and makes maximal use of the observational data. Statistical methods to estimate the zero-flux level and thus the absolute flux are also presented. The efficacy of the technique is demonstrated by an application to the light echoes around SN 1987A. The resulting images reveal new detail and faint light echo structures. This method can be adapted and applied to other extended variable light sources, such as stellar outflows and supernova remnants.

We assess the principal statistical and physical uncertainties associated with the determination of magnetic field strengths in clusters of galaxies from measurements of Faraday rotation (FR) and Compton-synchrotron emissions. In the former case a basic limitation is noted, that the relative uncertainty in the estimation of the mean-squared FR will generally be at least one third. Even greater uncertainty stems from the crucial dependence of the Faraday-deduced field on the coherence length scale characterizing its random orientation; we further elaborate this dependence, and argue that previous estimates of the field are likely to be too high by a factor of a few. Lack of detailed spatial information on the radio emission--and the recently deduced nonthermal X-ray emission in four clusters--has led to an underestimation of the mean value of the field in cluster cores. We conclude therefore that it is premature to draw definite quantitative conclusions from the previously-claimed seemingly-discrepant values of the field determined by these two methods.