We present the completed KMOS$^\mathrm{3D}$ survey $-$ an integral field spectroscopic survey of 739, $\log(M_{\star}/M_{\odot})>9$, galaxies at $0.6<z<2.7$ using the K-band Multi Object Spectrograph (KMOS) at the Very Large Telescope (VLT). KMOS$^\mathrm{3D}$ provides a population-wide census of kinematics, star formation, outflows, and nebular gas conditions both on and off the star-forming galaxy main sequence through the spatially resolved and integrated properties of H$\alpha$, [N II], and [S II] emission lines. We detect H$\alpha$ emission for 91% of galaxies on the main sequence of star-formation and 79% overall. The depth of the survey has allowed us to detect galaxies with star-formation rates below 1 M$_{\odot}$/ yr$^{-1}$, as well as to resolve 81% of detected galaxies with $\geq3$ resolution elements along the kinematic major axis. The detection fraction of H$\alpha$ is a strong function of both color and offset from the main sequence, with the detected and non-detected samples exhibiting different SED shapes. Comparison of H$\alpha$ and UV+IR star formation rates (SFRs) reveal that dust attenuation corrections may be underestimated by 0.5 dex at the highest masses ($\log(M_{\star}/M_{\odot})>10.5$). We confirm our first year results of a high rotation dominated fraction (monotonic velocity gradient and $v_\mathrm{rot}$/$\sigma_0 > \sqrt{3.36}$) of 77% for the full KMOS$^\mathrm{3D}$ H$\alpha$sample. The rotation-dominated fraction is a function of both stellar mass and redshift with the strongest evolution measured over the redshift range of the survey for galaxies with $\log(M_{\star}/M_{\odot})<10.5$. With this paper we include a final data release of all 739 observed objects.

We make publicly available a catalog of calibrated environmental measures for galaxies in the five 3D-HST/CANDELS deep fields. Leveraging the spectroscopic and grism redshifts from the 3D-HST survey, multi wavelength photometry from CANDELS, and wider field public data for edge corrections, we derive densities in fixed apertures to characterize the environment of galaxies brighter than $JH_{140} < 24$ mag in the redshift range $0.5<z<3.0$. By linking observed galaxies to a mock sample, selected to reproduce the 3D-HST sample selection and redshift accuracy, each 3D-HST galaxy is assigned a probability density function of the host halo mass, and a probability that is a central or a satellite galaxy. The same procedure is applied to a $z=0$ sample selected from SDSS. We compute the fraction of passive central and satellite galaxies as a function of stellar and halo mass, and redshift, and then derive the fraction of galaxies that were quenched by environment specific processes. Using the mock sample, we estimate that the timescale for satellite quenching is $t_{\rm quench} \sim 2-5$ Gyr; longer at lower stellar mass or lower redshift, but remarkably independent of halo mass. This indicates that, in the range of environments commonly found within the 3D-HST sample, satellites are quenched by exhaustion of their gas reservoir in absence of cosmological accretion. We find that the quenching times can be separated into a delay phase during which satellite galaxies behave similarly to centrals at fixed stellar mass, and a phase where the star formation rate drops rapidly ($\sim 0.4-0.6$ Gyr), as shown previously at $z=0$. We conclude that this scenario requires satellite galaxies to retain a large reservoir of multi-phase gas upon accretion, even at high redshift, and that this gas sustains star formation for the long quenching times observed.

We exploit deep integral-field spectroscopic observations with KMOS/VLT of 240 star-forming disks at 0.6 < z < 2.6 to dynamically constrain their mass budget. Our sample consists of massive ($\gtrsim 10^{9.8} M_\odot$) galaxies with sizes $R_e \gtrsim 2$ kpc. By contrasting the observed velocity and dispersion profiles to dynamical models, we find that on average the stellar content contributes $32^{+8}_{-7}\%$ of the total dynamical mass, with a significant spread among galaxies (68th percentile range f_star ~ 18 - 62%). Including molecular gas as inferred from CO- and dust-based scaling relations, the estimated baryonic mass adds up to $56^{+17}_{-12}\%$ of total for the typical galaxy in our sample, reaching ~ 90% at z > 2. We conclude that baryons make up most of the mass within the disk regions of high-redshift star-forming disk galaxies, with typical disks at z > 2 being strongly baryon-dominated within $R_e$. Substantial object-to-object variations in both stellar and baryonic mass fractions are observed among the galaxies in our sample, larger than what can be accounted for by the formal uncertainties in their respective measurements. In both cases, the mass fractions correlate most strongly with measures of surface density. High $\Sigma_{star}$ galaxies feature stellar mass fractions closer to unity, and systems with high inferred gas or baryonic surface densities leave less room for additional mass components other than stars and molecular gas. Our findings can be interpreted as more extended disks probing further (and more compact disks probing less far) into the dark matter halos that host them.

We analyze the angular momenta of massive star forming galaxies (SFGs) at the peak of the cosmic star formation epoch (z~0.8-2.6). Our sample of ~360 log(M*/Msun) ~ 9.3-11.8 SFGs is mainly based on the KMOS3D and SINS/zC-SINF surveys of H$\alpha$ kinematics, and collectively provides a representative subset of the massive star forming population. The inferred halo scale angular momentum distribution is broadly consistent with that theoretically predicted for their dark matter halos, in terms of mean spin parameter <$\lambda$> ~ 0.037 and its dispersion ($\sigma_{log(\lambda)}$~0.2). Spin parameters correlate with the disk radial scale, and with their stellar surface density, but do not depend significantly on halo mass, stellar mass, or redshift. Our data thus support the long-standing assumption that on average, even at high redshifts, the specific angular momentum of disk galaxies reflects that of their dark matter halos (j_d = j_DM). The lack of correlation between $\lambda$ x (j_d/j_DM) and the nuclear stellar density $\Sigma_{*}$(1kpc) favors a scenario where disk-internal angular momentum redistribution leads to "compaction" inside massive high-redshift disks. For our sample, the inferred average stellar-to-dark matter mass ratio is ~2%, consistent with abundance matching results. Including the molecular gas, the total baryonic disk-to-dark matter mass ratio is ~5% for halos near $10^{12}$ Msun, which corresponds to 31% of the cosmologically available baryons, implying that high-redshift disks are strongly baryon dominated.

We present measurements of the [NII]/Ha ratio as a probe of gas-phase oxygen abundance for a sample of 419 star-forming galaxies at z=0.6-2.7 from the KMOS3D near-IR multi-IFU survey. The mass-metallicity relation (MZR) is determined consistently with the same sample selection, metallicity tracer, and methodology over the wide redshift range probed by the survey. We find good agreement with long-slit surveys in the literature, except for the low-mass slope of the relation at z~2.3, where this sample is less biased than previous samples based on optical spectroscopic redshifts. In this regime we measure a steeper slope than some literature results. Excluding the AGN contribution from the MZR reduces sensitivity at the high mass end, but produces otherwise consistent results. There is no significant dependence of the [NII]/Ha ratio on SFR or environment at fixed redshift and stellar mass. The IFU data allow spatially resolved measurements of [NII]/Ha, from which we can infer abundance gradients for 180 galaxies, thus tripling the current sample in the literature. The observed gradients are on average flat, with only 15 gradients statistically offset from zero at >3sigma. We have modelled the effect of beam-smearing, assuming a smooth intrinsic radial gradient and known seeing, inclination and effective radius for each galaxy. Our seeing-limited observations can recover up to 70% of the intrinsic gradient for the largest, face-on disks, but only 30% for the smaller, more inclined galaxies. We do not find significant trends between observed or corrected gradients and any stellar population, dynamical or structural galaxy parameters, mostly in agreement with existing studies with much smaller sample sizes. In cosmological simulations, strong feedback is generally required to produce flat gradients at high redshift.

We investigate the stellar populations of 25 massive, galaxies ($\log[M_\ast/M_\odot] \geq 10.9$) at $1.5 < z < 2$ using data obtained with the K-band Multi-Object Spectrograph (KMOS) on the ESO VLT. Targets were selected to be quiescent based on their broadband colors and redshifts using data from the 3D-HST grism survey. The mean redshift of our sample is $\bar{z} = 1.75$, where KMOS YJ-band data probe age- and metallicity-sensitive absorption features in the rest-frame optical, including the $G$ band, Fe I, and high-order Balmer lines. Fitting simple stellar population models to a stack of our KMOS spectra, we derive a mean age of $1.03^{+0.13}_{-0.08}$ Gyr. We confirm previous results suggesting a correlation between color and age for quiescent galaxies, finding mean ages of $1.22^{+0.56}_{-0.19}$ Gyr and $0.85^{+0.08}_{-0.05}$ Gyr for the reddest and bluest galaxies in our sample. Combining our KMOS measurements with those obtained from previous studies at $0.2 < z < 2$ we find evidence for a $2-3$ Gyr spread in the formation epoch of massive galaxies. At $z < 1$ the measured stellar ages are consistent with passive evolution, while at $1 < z \lesssim2$ they appear to saturate at $\sim$1 Gyr, which likely reflects changing demographics of the (mean) progenitor population. By comparing to star-formation histories inferred for "normal" star-forming galaxies, we show that the timescales required to form massive galaxies at $z \gtrsim 1.5$ are consistent with the enhanced $\alpha$-element abundances found in massive local early-type galaxies.

We present the KMOS^3D survey, a new integral field survey of over 600 galaxies at 0.7<z<2.7 using KMOS at the Very Large Telescope (VLT). The KMOS^3D survey utilizes synergies with multi-wavelength ground and space-based surveys to trace the evolution of spatially-resolved kinematics and star formation from a homogeneous sample over 5 Gyrs of cosmic history. Targets, drawn from a mass-selected parent sample from the 3D-HST survey, cover the star formation-stellar mass ($M_*$) and rest-frame $(U-V)-M_*$ planes uniformly. We describe the selection of targets, the observations, and the data reduction. In the first year of data we detect Halpha emission in 191 $M_*=3\times10^{9}-7\times10^{11}$ Msun galaxies at z=0.7-1.1 and z=1.9-2.7. In the current sample 83% of the resolved galaxies are rotation-dominated, determined from a continuous velocity gradient and $v_{rot}/\sigma>1$, implying that the star-forming 'main sequence' (MS) is primarily composed of rotating galaxies at both redshift regimes. When considering additional stricter criteria, the Halpha kinematic maps indicate at least ~70% of the resolved galaxies are disk-like systems. Our high-quality KMOS data confirm the elevated velocity dispersions reported in previous IFS studies at z>0.7. For rotation-dominated disks, the average intrinsic velocity dispersion decreases by a factor of two from 50 km/s at z~2.3 to 25 km/s at z~0.9 while the rotational velocities at the two redshifts are comparable. Combined with existing results spanning z~0-3, disk velocity dispersions follow an approximate (1+z) evolution that is consistent with the dependence of velocity dispersion on gas fractions predicted by marginally-stable disk theory.

In this paper we follow up on our previous detection of nuclear ionized outflows in the most massive (log(M*/Msun) >= 10.9) z~1-3 star-forming galaxies (Forster Schreiber et al.), by increasing the sample size by a factor of six (to 44 galaxies above log(M*/Msun) >= 10.9) from a combination of the SINS/zC-SINF, LUCI, GNIRS, and KMOS^3D spectroscopic surveys. We find a fairly sharp onset of the incidence of broad nuclear emission (FWHM in the Ha, [NII], and [SII] lines ~ 450-5300 km/s), with large [NII]/Ha ratios, above log(M*/Msun) ~ 10.9, with about two thirds of the galaxies in this mass range exhibiting this component. Broad nuclear components near and above the Schechter mass are similarly prevalent above and below the main sequence of star-forming galaxies, and at z~1 and ~2. The line ratios of the nuclear component are fit by excitation from active galactic nuclei (AGN), or by a combination of shocks and photoionization. The incidence of the most massive galaxies with broad nuclear components is at least as large as that of AGNs identified by X-ray, optical, infrared or radio indicators. The mass loading of the nuclear outflows is near unity. Our findings provide compelling evidence for powerful, high-duty cycle, AGN-driven outflows near the Schechter mass, and acting across the peak of cosmic galaxy formation.

We present the correlations between stellar mass, star formation rate (SFR) and [NII]/Ha flux ratio as indicator of gas-phase metallicity for a sample of 222 galaxies at 0.8 < z < 2.6 and log(M*/Msun)=9.0-11.5 from the LUCI, SINS/zC-SINF and KMOS3D surveys. This sample provides a unique analysis of the mass-metallicity relation (MZR) over an extended redshift range using consistent data analysis techniques and strong-line metallicity indicator. We find a constant slope at the low-mass end of the relation and can fully describe its redshift evolution through the evolution of the characteristic turnover mass where the relation begins to flatten at the asymptotic metallicity. At fixed mass and redshift, our data do not show a correlation between the [NII]/Ha ratio and SFR, which disagrees with the 0.2-0.3dex offset in [NII]/Ha predicted by the "fundamental relation" between stellar mass, SFR and metallicity discussed in recent literature. However, the overall evolution towards lower [NII]/Ha at earlier times does broadly agree with these predictions.

We report the detection of ubiquitous powerful nuclear outflows in massive (> 10^11 Msun) z~2 star-forming galaxies (SFGs), which are plausibly driven by an Active Galactic Nucleus (AGN). The sample consists of the eight most massive SFGs from our SINS/zC-SINF survey of galaxy kinematics with the imaging spectrometer SINFONI, six of which have sensitive high-resolution adaptive optics (AO) assisted observations. All of the objects are disks hosting a significant stellar bulge. The spectra in their central regions exhibit a broad component in Halpha and forbidden [NII] and [SII] line emission, with typical velocity FWHM ~ 1500 km/s, [NII]/Halpha ratio ~ 0.6, and intrinsic extent of 2 - 3 kpc. These properties are consistent with warm ionized gas outflows associated with Type 2 AGN, the presence of which is confirmed via independent diagnostics in half the galaxies. The data imply a median ionized gas mass outflow rate of ~ 60 Msun/yr and mass loading of ~ 3. At larger radii, a weaker broad component is detected but with lower FWHM ~ 485 km/s and [NII]/Halpha ~ 0.35, characteristic for star formation-driven outflows as found in the lower-mass SINS/zC-SINF galaxies. The high inferred mass outflow rates and frequent occurrence suggest the nuclear outflows efficiently expel gas out of the centers of the galaxies with high duty cycles, and may thus contribute to the process of star formation quenching in massive galaxies. Larger samples at high masses will be crucial to confirm the importance and energetics of the nuclear outflow phenomenon, and its connection to AGN activity and bulge growth.

As part of the SINS/zC-SINF surveys of high-z galaxy kinematics, we derive the radial distributions of H-alpha surface brightness, stellar mass surface density, and dynamical mass at ~2 kpc resolution in 19 z~2 star-forming disks with deep SINFONI AO spectroscopy at the ESO VLT. From these data we infer the radial distribution of the Toomre Q-parameter for these main-sequence star forming galaxies (SFGs), covering almost two decades of stellar mass (10^9.6 to 10^11.5 solar masses). In more than half of our SFGs, the H-alpha distributions cannot be fit by a centrally peaked distribution, such as an exponential, but are better described by a ring, or the combination of a ring and an exponential. At the same time the kinematic data indicate the presence of a mass distribution more centrally concentrated than a single exponential distribution for 5 of the 19 galaxies. The resulting Q-distributions are centrally peaked for all, and significantly exceed unity there for three quarters of the SFGs. The occurrence of H-alpha rings and of large nuclear Q-values is strongly correlated, and is more common for the more massive SFGs. While our sample is small and there remain substantial uncertainties and caveats, our observations are consistent with a scenario in which cloud fragmentation and global star formation are secularly suppressed in gas rich high-z disks from the inside out, as the central stellar mass density of the disks grows.

We combine IRAM Plateau de Bure Interferometer and Herschel PACS and SPIRE measurements to study the dust and gas contents of high-redshift star forming galaxies. We present new observations for a sample of 17 lensed galaxies at z=1.4-3.1, which allow us to directly probe the cold ISM of normal star-forming galaxies with stellar masses of ~10^10Msun, a regime otherwise not (yet) accessible by individual detections in Herschel and molecular gas studies. The lensed galaxies are combined with reference samples of sub-millimeter and normal z~1-2 star-forming galaxies with similar far-infrared photometry to study the gas and dust properties of galaxies in the SFR-M*-redshift parameter space. The mean gas depletion timescale of main sequence galaxies at z>2 is measured to be only ~450Myr, a factor of ~1.5 (~5) shorter than at z=1 (z=0), in agreement with a (1+z)^-1 scaling. The mean gas mass fraction at z=2.8 is 40+/-15% (44% after incompleteness correction), suggesting a flattening or even a reversal of the trend of increasing gas fractions with redshift recently observed up to z~2. The depletion timescale and gas fractions of the z>2 normal star-forming galaxies can be explained under the "equilibrium model" for galaxy evolution, in which the gas reservoir of galaxies is the primary driver of the redshift evolution of specific star formation rates. Due to their high star formation efficiencies and low metallicities, the z>2 lensed galaxies have warm dust despite being located on the star formation main sequence. At fixed metallicity, they also have a gas-to-dust ratio 1.7 times larger than observed locally when using the same standard techniques, suggesting that applying the local calibration of the relation between gas-to-dust ratio and metallicity to infer the molecular gas mass of high redshift galaxies may lead to systematic differences with CO-based estimates.

We investigate the correlation between the mass of a central supermassive black hole and the total gravitational mass of the host galaxy (M_tot). The results are based on 43 galaxy-scale strong gravitational lenses from the Sloan Lens ACS (SLACS) Survey whose black hole masses were estimated through two scaling relations: the relation between black hole mass and Sersic index (M_bh - n) and the relation between black hole mass and stellar velocity dispersion (M_bh - sigma). We use the enclosed mass within R_200, the radius within which the density profile of the early type galaxy exceeds the critical density of the Universe by a factor of 200, determined by gravitational lens models fitted to HST imaging data, as a tracer of the total gravitational mass. The best fit correlation, where M_bh is determined from M_bh - sigma relation, is log(M_bh) = (8.18 +/- 0.11) + (1.55 +/- 0.31) (log(M_tot) - 13.0) over 2 orders of magnitude in M_bh. From a variety of tests, we find that we cannot reliably infer a connection between M_bh and M_tot from the M_bh - n relation. The M_bh - M_tot relation provides some of the first, direct observational evidence to test the prediction that supermassive black hole properties are determined by the halo properties of the host galaxy.