The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

We report the spectroscopic confirmation of a new protocluster in the COSMOS field at $z$ $\sim$ 2.2, COSMOS Cluster 2.2 (CC2.2), originally identified as an overdensity of narrowband selected H$\alpha$ emitting candidates. With only two masks of Keck/MOSFIRE near-IR spectroscopy in both $H$ ($\sim$ 1.47-1.81 $\mu$m) and $K$ ($\sim$ 1.92-2.40 $\mu$m) bands ($\sim$ 1.5 hour each), we confirm 35 unique protocluster members with at least two emission lines detected with S/N $>$ 3. Combined with 12 extra members from the zCOSMOS-deep spectroscopic survey (47 in total), we estimate a mean redshift and a line-of-sight velocity dispersion of $z_{mean}$=2.23224 $\pm$ 0.00101 and $\sigma_{los}$=645 $\pm$ 69 km s$^{-1}$ for this protocluster, respectively. Assuming virialization and spherical symmetry for the system, we estimate a total mass of $M_{vir}$ $\sim$ $(1-2) \times$10$^{14}$ $M_{\odot}$ for the structure. We evaluate a number density enhancement of $\delta_{g}$ $\sim$ 7 for this system and we argue that the structure is likely not fully virialized at $z$ $\sim$ 2.2. However, in a spherical collapse model, $\delta_{g}$ is expected to grow to a linear matter enhancement of $\sim$ 1.9 by $z$=0, exceeding the collapse threshold of 1.69, and leading to a fully collapsed and virialized Coma-type structure with a total mass of $M_{dyn}$($z$=0) $\sim$ 9.2$\times$10$^{14}$ $M_{\odot}$ by now. This observationally efficient confirmation suggests that large narrowband emission-line galaxy surveys, when combined with ancillary photometric data, can be used to effectively trace the large-scale structure and protoclusters at a time when they are mostly dominated by star-forming galaxies.

We explore synergies between the space-based Wide-Field Infrared Survey Telescope (WFIRST) and the ground-based Rubin Observatory Legacy Survey of Space and Time (LSST). In particular, we consider a scenario where the currently envisioned survey strategy for WFIRST's High Latitude Survey (HLS), i.e., 2000 square degrees in four narrow photometric bands is altered in favor of a strategy that combines rapid coverage of the LSST area (to full LSST depth) in one band. We find that a 5-month WFIRST survey in the W-band can cover the full LSST survey area providing high-resolution imaging for >95% of the LSST Year 10 gold galaxy sample. We explore a second, more ambitious scenario where WFIRST spends 1.5 years covering the LSST area. For this second scenario we quantify the constraining power on dark energy equation of state parameters from a joint weak lensing and galaxy clustering analysis, and compare it to an LSST-only survey and to the Reference WFIRST HLS survey. Our survey simulations are based on the WFIRST exposure time calculator and redshift distributions from the CANDELS catalog. Our statistical uncertainties account for higher-order correlations of the density field, and we include a wide range of systematic effects, such as uncertainties in shape and redshift measurements, and modeling uncertainties of astrophysical systematics, such as galaxy bias, intrinsic galaxy alignment, and baryonic physics. Assuming the 5-month WFIRST wide scenario, we find a significant increase in constraining power for the joint LSST+WFIRST wide survey compared to LSST Y10 (FoM(Wwide)= 2.4 FoM(LSST)) and compared to LSST+WFIRST HLS (FoM(Wwide)= 5.5 FoM(HLS)).

We simulate the scientific performance of the Wide-Field Infrared Survey Telescope (WFIRST) High Latitude Survey (HLS) on dark energy and modified gravity. The 1.6 year HLS Reference survey is currently envisioned to image 2000 deg$^2$ in multiple bands to a depth of $\sim$26.5 in Y, J, H and to cover the same area with slit-less spectroscopy beyond z=3. The combination of deep, multi-band photometry and deep spectroscopy will allow scientists to measure the growth and geometry of the Universe through a variety of cosmological probes (e.g., weak lensing, galaxy clusters, galaxy clustering, BAO, Type Ia supernova) and, equally, it will allow an exquisite control of observational and astrophysical systematic effects. In this paper we explore multi-probe strategies that can be implemented given WFIRST's instrument capabilities. We model cosmological probes individually and jointly and account for correlated systematics and statistical uncertainties due to the higher order moments of the density field. We explore different levels of observational systematics for the WFIRST survey (photo-z and shear calibration) and ultimately run a joint likelihood analysis in N-dim parameter space. We find that the WFIRST reference survey alone (no external data sets) can achieve a standard dark energy FoM of >300 when including all probes. This assumes no information from external data sets and realistic assumptions for systematics. Our study of the HLS reference survey should be seen as part of a future community driven effort to simulate and optimize the science return of WFIRST.

The Euclid mission will observe well over a billion galaxies out to $z\sim6$ and beyond. This will offer an unrivalled opportunity to investigate several key questions for understanding galaxy formation and evolution. The first step for many of these studies will be the selection of a sample of quiescent and star-forming galaxies, as is often done in the literature by using well known colour techniques such as the `UVJ' diagram. However, given the limited number of filters available for the Euclid telescope, the recovery of such rest-frame colours will be challenging. We therefore investigate the use of observed Euclid colours, on their own and together with ground-based u-band observations, for selecting quiescent and star-forming galaxies. The most efficient colour combination, among the ones tested in this work, consists of the (u-VIS) and (VIS-J) colours. We find that this combination allows users to select a sample of quiescent galaxies complete to above $\sim70\%$ and with less than 15$\%$ contamination at redshifts in the range $0.75<z<1$. For galaxies at high-z or without the u-band complementary observations, the (VIS-Y) and (J-H) colours represent a valid alternative, with $>65\%$ completeness level and contamination below 20$\%$ at $1<z<2$ for finding quiescent galaxies. In comparison, the sample of quiescent galaxies selected with the traditional UVJ technique is only $\sim20\%$ complete at $z<3$, when recovering the rest-frame colours using mock Euclid observations. This shows that our new methodology is the most suitable one when only Euclid bands, along with u-band imaging, are available.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

Large photometric surveys provide a rich source of observations of quiescent galaxies, including a surprisingly large population at z>1. However, identifying large, but clean, samples of quiescent galaxies has proven difficult because of their near-degeneracy with interlopers such as dusty, star-forming galaxies. We describe a new technique for selecting quiescent galaxies based upon t-distributed stochastic neighbor embedding (t-SNE), an unsupervised machine learning algorithm for dimensionality reduction. This t-SNE selection provides an improvement both over UVJ, removing interlopers which otherwise would pass color selection, and over photometric template fitting, more strongly towards high redshift. Due to the similarity between the colors of high- and low-redshift quiescent galaxies, under our assumptions t-SNE outperforms template fitting in 63% of trials at redshifts where a large training sample already exists. It also may be able to select quiescent galaxies more efficiently at higher redshifts than the training sample.

Using the Horizon-AGN hydrodynamical simulation and self-organising maps (SOMs), we show how to compress the complex data structure of a cosmological simulation into a 2-d grid which is much easier to analyse. We first verify the tight correlation between the observed 0.3$\!-\!5\mu$m broad-band colours of Horizon-AGN galaxies and their high-resolution spectra. The correlation is found to extend to physical properties such as redshift, stellar mass, and star formation rate (SFR). This direct mapping from colour to physical parameter space is shown to work also after including photometric uncertainties that mimic the COSMOS survey. We then label the SOM grid with a simulated calibration sample and estimate redshift and SFR for COSMOS-like galaxies up to $z\sim3$. In comparison to state-of-the-art techniques based on synthetic templates, our method is comparable in performance but less biased at estimating redshifts, and significantly better at predicting SFRs. In particular our "data-driven" approach, in contrast to model libraries, intrinsically allows for the complexity of galaxy formation and can handle sample biases. We advocate that obtaining the calibration for this method should be one of the goals of next-generation galaxy surveys.

Wide-field imaging surveys such as the Dark Energy Survey (DES) rely on coarse measurements of spectral energy distributions in a few filters to estimate the redshift distribution of source galaxies. In this regime, sample variance, shot noise, and selection effects limit the attainable accuracy of redshift calibration and thus of cosmological constraints. We present a new method to combine wide-field, few-filter measurements with catalogs from deep fields with additional filters and sufficiently low photometric noise to break degeneracies in photometric redshifts. The multi-band deep field is used as an intermediary between wide-field observations and accurate redshifts, greatly reducing sample variance, shot noise, and selection effects. Our implementation of the method uses self-organizing maps to group galaxies into phenotypes based on their observed fluxes, and is tested using a mock DES catalog created from N-body simulations. It yields a typical uncertainty on the mean redshift in each of five tomographic bins for an idealized simulation of the DES Year 3 weak-lensing tomographic analysis of $\sigma_{\Delta z} = 0.007$, which is a 60% improvement compared to the Year 1 analysis. Although the implementation of the method is tailored to DES, its formalism can be applied to other large photometric surveys with a similar observing strategy.

High-multiplex and deep spectroscopic follow-up of upcoming panoramic deep-imaging surveys like LSST, Euclid, and WFIRST is a widely recognized and increasingly urgent necessity. No current or planned facility at a U.S. observatory meets the sensitivity, multiplex, and rapid-response time needed to exploit these future datasets. FOBOS, the Fiber-Optic Broadband Optical Spectrograph, is a near-term fiber-based facility that addresses these spectroscopic needs by optimizing depth over area and exploiting the aperture advantage of the existing 10m Keck II Telescope. The result is an instrument with a uniquely blue-sensitive wavelength range (0.31-1.0 um) at R~3500, high-multiplex (1800 fibers), and a factor 1.7 greater survey speed and order-of-magnitude greater sampling density than Subaru's Prime Focus Spectrograph (PFS). In the era of panoramic deep imaging, FOBOS will excel at building the deep, spectroscopic reference data sets needed to interpret vast imaging data. At the same time, its flexible focal plane, including a mode with 25 deployable integral-field units (IFUs) across a 20 arcmin diameter field, enables an expansive range of scientific investigations. Its key programmatic areas include (1) nested stellar-parameter training sets that enable studies of the Milky Way and M31 halo sub-structure, as well as local group dwarf galaxies, (2) a comprehensive picture of galaxy formation thanks to detailed mapping of the baryonic environment at z~2 and statistical linking of evolving populations to the present day, and (3) dramatic enhancements in cosmological constraints via precise photometric redshifts and determined redshift distributions. In combination with Keck I instrumentation, FOBOS also provides instant access to medium-resolution spectroscopy for transient sources with full coverage from the UV to the K-band.

Although photometric redshifts (photo-z's) are crucial ingredients for current and upcoming large-scale surveys, the high-quality spectroscopic redshifts currently available to train, validate, and test them are substantially non-representative in both magnitude and color. We investigate the nature and structure of this bias by tracking how objects from a heterogeneous training sample contribute to photo-z predictions as a function of magnitude and color, and illustrate that the underlying redshift distribution at fixed color can evolve strongly as a function of magnitude. We then test the robustness of the galaxy-galaxy lensing signal in 120 deg$^2$ of HSC-SSP DR1 data to spectroscopic completeness and photo-z biases, and find that their impacts are sub-dominant to current statistical uncertainties. Our methodology provides a framework to investigate how spectroscopic incompleteness can impact photo-z-based weak lensing predictions in future surveys such as LSST and WFIRST.

We show unsupervised machine learning techniques are a valuable tool for both visualizing and computationally accelerating the estimation of galaxy physical properties from photometric data. As a proof of concept, we use self organizing maps (SOMs) to visualize a spectral energy distribution (SED) model library in the observed photometry space. The resulting visual maps allow for a better understanding of how the observed data maps to physical properties and to better optimize the model libraries for a given set of observational data. Next, the SOMs are used to estimate the physical parameters of 14,000 z~1 galaxies in the COSMOS field and found to be in agreement with those measured with SED fitting. However, the SOM method is able to estimate the full probability distribution functions for each galaxy up to about a million times faster than direct model fitting. We conclude by discussing how this speed up and learning how the galaxy data manifold maps to physical parameter space and visualizing this mapping in lower dimensions helps overcome other challenges in galaxy formation and evolution.

In order for Wide-Field Infrared Survey Telescope (WFIRST) and other Stage IV dark energy experiments (e.g., Large Synoptic Survey Telescope; LSST, and Euclid) to infer cosmological parameters not limited by systematic errors, accurate redshift measurements are needed. This accuracy can be met by using spectroscopic subsamples to calibrate the photometric redshifts for the full sample. In this work we employ the Self Organizing Map (SOM) spectroscopic sampling technique, to find the minimal number of spectra required for the WFIRST weak lensing calibration. We use galaxies from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) to build the LSST+WFIRST lensing analog sample of ~36 k objects and train the LSST+WFIRST SOM. We find that 26% of the WFIRST lensing sample consists of sources fainter than the Euclid depth in the optical, 91% of which live in color cells already occupied by brighter galaxies. We demonstrate the similarity between faint and bright galaxies as well as the feasibility of redshift measurements at different brightness levels. Our results suggest that the spectroscopic sample acquired for calibration to the Euclid depth is sufficient for calibrating the majority of the WFIRST color-space. For the spectroscopic sample to fully represent the synthetic color-space of WFIRST, we recommend obtaining additional spectroscopy of ~0.2-1.2 k new sources in cells occupied by mostly faint galaxies. We argue that either the small area of the CANDELS fields and the small overall sample size or the large photometric errors might be the reason for no/less bright galaxies mapped to these cells. Acquiring the spectra of these sources will confirm the above findings and will enable the comprehensive calibration of the WFIRST color-redshift relation.

The Complete Calibration of the Color-Redshift Relation (C3R2) survey is a multi-institution, multi-instrument survey that aims to map the empirical relation of galaxy color to redshift to i~24.5 (AB), thereby providing a firm foundation for weak lensing cosmology with the Stage IV dark energy missions Euclid and WFIRST. Here we present 3171 new spectroscopic redshifts obtained in the 2016B and 2017A semesters with a combination of DEIMOS, LRIS, and MOSFIRE on the Keck telescopes. The observations come from all of the Keck partners: Caltech, NASA, the University of Hawaii, and the University of California. Combined with the 1283 redshifts published in DR1, the C3R2 survey has now obtained and published 4454 high quality galaxy redshifts. We discuss updates to the survey design and provide a catalog of photometric and spectroscopic data. Initial tests of the calibration method performance are given, indicating that the sample, once completed and combined with extensive data collected by other spectroscopic surveys, should allow us to meet the cosmology requirements for Euclid, and make significant headway toward solving the problem for WFIRST. We use the full spectroscopic sample to demonstrate that galaxy brightness is weakly correlated with redshift once a galaxy is localized in the Euclid or WFIRST color space, with potentially important implications for the spectroscopy needed to calibrate redshifts for faint WFIRST and LSST sources.

Two decades after its discovery, cosmic acceleration remains the most profound mystery in cosmology and arguably in all of physics. Either the Universe is dominated by a form of dark energy with exotic physical properties not predicted by standard model physics, or General Relativity is not an adequate description of gravity over cosmic distances. WFIRST emerged as a top priority of Astro2010 in part because of its ability to address the mystery of cosmic acceleration through both high precision measurements of the cosmic expansion history and the growth of cosmic structures with multiple and redundant probes. We illustrate in this white paper how mission design changes since Astro2010 have made WFIRST an even more powerful dark energy facility and have improved the ability of WFIRST to respond to changes in the experimental landscape. WFIRST is the space-based probe of DE the community needs in the mid-2020s.

The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.

In this paper we address the problem of finding optimal cosmic shear tomographic bins. We generalise the definition of a cosmic shear tomographic bin to be a set of commonly labelled voxels in photometric colour space; rather than bins defined directly in redshift. We explore this approach by using a self-organising map to define the multi-dimensional colour space, and a we define a 'label space' of connected regions on the self-organising map using overlapping elliptical disks. This allows us to then find optimal labelling schemes by searching the label space. We use a metric that is the signal-to-noise ratio of a dark energy equation of state measurement, and in this case we find that for up to five tomographic bins the optimal colour-space labelling is an approximation of an equally-spaced binning in redshift; that is in all cases the best configuration. We also show that such a redefinition is more robust to photometric redshift outliers than a standard tomographic bin selection.

Cosmic acceleration is the most surprising cosmological discovery in many decades. Testing and distinguishing among possible explanations requires cosmological measurements of extremely high precision probing the full history of cosmic expansion and structure growth and, ideally, compare and contrast matter and relativistic tracers of the gravity potential. This program is one of the defining objectives of the Wide-Field Infrared Survey Telescope (WFIRST), as set forth in the New Worlds, New Horizons report (NWNH) in 2010. The WFIRST mission has the ability to improve these measurements by 1-2 orders of magnitude compared to the current state of the art, while simultaneously extending their redshift grasp, greatly improving control of systematic effects, and taking a unified approach to multiple probes that provide complementary physical information and cross-checks of cosmological results. We describe in this annual report the activities of the Science Investigation Team (SIT) "Cosmology with the High Latitude Survey (HLS)" during the year 2017. This team was selected by NASA in December 2015 in order to address the stringent challenges of the WFIRST dark energy (DE) program through the Project's formulation phase. This SIT has elected to jointly address Galaxy Redshift Survey, Weak Lensing and Cluster Growth and thus fully embrace the fact that the imaging and spectroscopic elements of the HLS will be realized as an integrated observing program, and they jointly impose requirements on performance and operations. WFIRST is designed to be able to deliver a definitive result on the origin of cosmic acceleration. It is not optimized for Figure of Merit sensitivity but for control of systematic uncertainties and for having multiple techniques each with multiple cross-checks. Our SIT work focuses on understanding the potential systematics in the WFIRST DE measurements.

We present a catalog of 10718 objects in the COSMOS field observed through multi-slit spectroscopy with the Deep Imaging Multi-Object Spectrograph (DEIMOS) on the Keck II telescope in the wavelength range ~5500-9800A. The catalog contains 6617 objects with high-quality spectra (two or more spectral features), and 1798 objects with a single spectroscopic feature confirmed by the photometric redshift. For 2024 typically faint objects we could not obtain reliable redshifts. The objects have been selected from a variety of input catalogs based on multi-wavelength observations in the field, and thus have a diverse selection function, which enables the study of the diversity in the galaxy population. The magnitude distribution of our objects is peaked at I_AB~23 and K_AB~21, with a secondary peak at K_AB~24. We sample a broad redshift distribution in the range 0<z<6, with one peak at z~1, and another one around z~4. We have identified 13 redshift spikes at z>0.65 with chance probabilities <4xE-4$, some of which are clearly related to protocluster structures of sizes >10 Mpc. An object-to-object comparison with a multitude of other spectroscopic samples in the same field shows that our DEIMOS sample is among the best in terms of fraction of spectroscopic failures and relative redshift accuracy. We have determined the fraction of spectroscopic blends to about 0.8% in our sample. This is likely a lower limit and at any rate well below the most pessimistic expectations. Interestingly, we find evidence for strong lensing of Ly-alpha background emitters within the slits of 12 of our target galaxies, increasing their apparent density by about a factor of 4.

We present an empirical parameterization of the [NII]/H$\alpha$ flux ratio as a function of stellar mass and redshift valid at 0 < z < 2.7 and 8.5 < log(M) < 11.0. This description can easily be applied to (i) simulations for modeling [NII]$\lambda6584$ line emission, (ii) deblend [NII] and H$\alpha$ in current low-resolution grism and narrow-band observations to derive intrinsic H$\alpha$ fluxes, and (iii) to reliably forecast the number counts of H$\alpha$ emission-line galaxies for future surveys, such as those planned for Euclid and the Wide Field Infrared Survey Telescope (WFIRST). Our model combines the evolution of the locus on the Baldwin, Phillips & Terlevich (BPT) diagram measured in spectroscopic data out to z = 2.5 with the strong dependence of [NII]/H$\alpha$ on stellar mass and [OIII]/H$\beta$ observed in local galaxy samples. We find large variations in the [NII]/H$\alpha$ flux ratio at a fixed redshift due to its dependency on stellar mass; hence, the assumption of a constant [NII] flux contamination fraction can lead to a significant under- or overestimate of H$\alpha$ luminosities. Specifically, measurements of the intrinsic H$\alpha$ luminosity function derived from current low-resolution grism spectroscopy assuming a constant 29% contamination of [NII] can be overestimated by factors of ~8 at log(L) > 43.0 for galaxies at redshifts z = 1.5. This has implications for the prediction of H$\alpha$ emitters for Euclid and WFIRST. We also study the impact of blended H$\alpha$ and [NII] on the accuracy of measured spectroscopic redshifts.

We present a multi-wavelength catalog in the Subaru-XMM Deep Field (SXDF) as part of the Spitzer Large Area Survey with Hyper-Suprime-Cam (SPLASH). We include the newly acquired optical data from the Hyper-Suprime Cam Subaru Strategic Program, accompanied by IRAC coverage from the SPLASH survey. All available optical and near-infrared data is homogenized and resampled on a common astrometric reference frame. Source detection is done using a multi-wavelength detection image including the $u$-band to recover the bluest objects. We measure multi-wavelength photometry and compute photometric redshifts as well as stellar masses for $\sim$1.17 million objects over $\sim$4.2 deg$^2$ with $\sim$800,000 objects in the 2.4 deg$^2$ HSC-UltraDeep coverage. Using the available spectroscopic redshifts from various surveys over the range of $0<z<6$, we verify the performance of the photometric redshifts and we find a normalized median absolute deviation of 0.023 and outlier fraction of 3.2%. The SPLASH-SXDF catalog is a valuable, publicly available resource that is perfectly suited for studying galaxies in the early universe and tracing their evolution through cosmic time. The catalog is available for download from https://z.umn.edu/SXDF .

Knowledge of the number density of H$\alpha$ emitting galaxies is vital for assessing the scientific impact of the Euclid and WFIRST missions. In this work we present predictions from a galaxy formation model, Galacticus, for the cumulative number counts of H$\alpha$-emitting galaxies. We couple Galacticus to three different dust attenuation methods and examine the counts using each method. A $\chi^2$ minimisation approach is used to compare the model predictions to observed galaxy counts and calibrate the dust parameters. We find that weak dust attenuation is required for the Galacticus counts to be broadly consistent with the observations, though the optimum dust parameters return large values for $\chi^2$, suggesting that further calibration of Galacticus is necessary. The model predictions are also consistent with observed estimates for the optical depth and the H$\alpha$ luminosity function. Finally we present forecasts for the redshift distributions and number counts for two Euclid-like and one WFIRST-like survey. For a Euclid-like survey with redshift range $0.9\leqslant z\leqslant 1.8$ and H$\alpha+{\rm [NII]}$ blended flux limit of $2\times 10^{-16}{\rm erg}\,{\rm s}^{-1}\,{\rm cm}^{-2}$ we predict a number density between 3900--4800 galaxies per square degree. For a WFIRST-like survey with redshift range $1\leqslant z\leqslant 2$ and blended flux limit of $1\times 10^{-16}{\rm erg}\,{\rm s}^{-1}\,{\rm cm}^{-2}$ we predict a number density between 10400--15200 galaxies per square degree.

Recent studies have found a significant evolution and scatter in the IRX-$\beta$ relation at z > 4, suggesting different dust properties of these galaxies. The total far-infrared (FIR) luminosity is key for this analysis but poorly constrained in normal (main-sequence) star-forming z > 5 galaxies where often only one single FIR point is available. To better inform estimates of the FIR luminosity, we construct a sample of local galaxies and three low-redshift analogs of z > 5 systems. The trends in this sample suggest that normal high-redshift galaxies have a warmer infrared (IR) SED compared to average z < 4 galaxies that are used as prior in these studies. The blue-shifted peak and mid-IR excess emission could be explained by a combination of a larger fraction of the metal-poor inter-stellar medium (ISM) being optically thin to ultra-violet (UV) light and a stronger UV radiation field due to high star formation densities. Assuming a maximally warm IR SED suggests 0.6 dex increased total FIR luminosities, which removes some tension between dust attenuation models and observations of the IRX-$\beta$ relation at z > 5. Despite this, some galaxies still fall below the minimum IRX-$\beta$ relation derived with standard dust cloud models. We propose that radiation pressure in these highly star-forming galaxies causes a spatial offset between dust clouds and young star-forming regions within the lifetime of O/B stars. These offsets change the radiation balance and create viewing-angle effects that can change UV colors at fixed IRX. We provide a modified model that can explain the location of these galaxies on the IRX-$\beta$ diagram.

We measure the stellar mass function (SMF) of galaxies in the COSMOS field up to $z\sim6$. We select them in the near-IR bands of the COSMOS2015 catalogue, which includes ultra-deep photometry from UltraVISTA-DR2, SPLASH, and Subaru/Hyper-SuprimeCam. At $z>2.5$ we use new precise photometric redshifts with error $\sigma_z=0.03(1+z)$ and an outlier fraction of $12\%$, estimated by means of the unique spectroscopic sample of COSMOS. The increased exposure time in the DR2, along with our panchromatic detection strategy, allow us to improve the stellar mass completeness at high $z$ with respect to previous UltraVISTA catalogues. We also identify passive galaxies through a robust colour-colour selection, extending their SMF estimate up to $z=4$. Our work provides a comprehensive view of galaxy stellar mass assembly between $z=0.1$ and 6, for the first time using consistent estimates across the entire redshift range. We fit these measurements with a Schechter function, correcting for Eddington bias. We compare the SMF fit with the halo mass function predicted from $\Lambda$CDM simulations. We find that at $z>3$ both functions decline with a similar slope in the high-mass end. This feature could be explained assuming that the mechanisms that quench star formation in massive haloes become less effective at high redshift; however further work needs to be done to confirm this scenario. Concerning the SMF low-mass end, it shows a progressive steepening as moving towards higher redshifts, with $\alpha$ decreasing from $-1.47_{-0.02}^{+0.02}$ at $z\simeq0.1$ to $-2.11_{-0.13}^{+0.30}$ at $z\simeq5$. This slope depends on the characterisation of the observational uncertainties, which is crucial to properly remove the Eddington bias. We show that there is currently no consensus on the method to quantify such errors: different error models result in different best-fit Schechter parameters. [Abridged]

A key goal of the Stage IV dark energy experiments Euclid, LSST and WFIRST is to measure the growth of structure with cosmic time from weak lensing analysis over large regions of the sky. Weak lensing cosmology will be challenging: in addition to highly accurate galaxy shape measurements, statistically robust and accurate photometric redshift (photo-z) estimates for billions of faint galaxies will be needed in order to reconstruct the three-dimensional matter distribution. Here we present an overview of and initial results from the Complete Calibration of the Color-Redshift Relation (C3R2) survey, designed specifically to calibrate the empirical galaxy color-redshift relation to the Euclid depth. These redshifts will also be important for the calibrations of LSST and WFIRST. The C3R2 survey is obtaining multiplexed observations with Keck (DEIMOS, LRIS, and MOSFIRE), the Gran Telescopio Canarias (GTC; OSIRIS), and the Very Large Telescope (VLT; FORS2 and KMOS) of a targeted sample of galaxies most important for the redshift calibration. We focus spectroscopic efforts on under-sampled regions of galaxy color space identified in previous work in order to minimize the number of spectroscopic redshifts needed to map the color-redshift relation to the required accuracy. Here we present the C3R2 survey strategy and initial results, including the 1283 high confidence redshifts obtained in the 2016A semester and released as Data Release 1.

ALMA observations of the long wavelength dust continuum are used to estimate the interstellar medium (ISM) masses in a sample of 708 galaxies at z = 0.3 to 4.5 in the COSMOS field. The galaxy sample has known far-infrared luminosities and, hence, star formation rates (SFRs), and stellar masses (M$_{\rm *}$) from the optical-infrared spectrum fitting. The galaxies sample SFRs from the main sequence (MS) to 50 times above the MS. The derived ISM masses are used to determine the dependence of gas mass on redshift, M$_{\rm *}$, and specific SFR (sSFR) relative to the MS. The ISM masses increase approximately 0.63 power of the rate of increase in SFRs with redshift and the 0.32 power of the sSFR/sSFR$_MS$. The SF efficiencies also increase as the 0.36 power of the SFR redshift evolutionary and the 0.7 power of the elevation above the MS; thus the increased activities at early epochs are driven by both increased ISM masses and SF efficiency. Using the derived ISM mass function we estimate the accretion rates of gas required to maintain continuity of the MS evolution ($>100$ \msun yr$^{-1}$ at z $>$ 2.5). Simple power-law dependences are similarly derived for the gas accretion rates. We argue that the overall evolution of galaxies is driven by the rates of gas accretion. The cosmic evolution of total ISM mass is estimated and linked to the evolution of SF and AGN activity at early epochs.

We investigate the physical conditions of ionized gas in high-z star-forming galaxies using diagnostic diagrams based on the rest-frame optical emission lines. The sample consists of 701 galaxies with an Ha detection at $1.4\lesssim z\lesssim1.7$, from the FMOS-COSMOS survey, that represent the normal star-forming population over the stellar mass range $10^{9.6} \lesssim M_\ast/M_\odot \lesssim 10^{11.6}$ with those at $M_\ast>10^{11}~M_\odot$ being well sampled. We confirm an offset of the average location of star-forming galaxies in the BPT diagram ([OIII]/Hb vs. [NII]/Ha), primarily towards higher [OIII]/Hb, compared with local galaxies. Based on the [SII] ratio, we measure an electron density ($n_e=220^{+170}_{-130}~\mathrm{cm^{-3}}$), that is higher than that of local galaxies. Based on comparisons to theoretical models, we argue that changes in emission-line ratios, including the offset in the BPT diagram, are caused by a higher ionization parameter both at fixed stellar mass and at fixed metallicity with additional contributions from a higher gas density and possibly a hardening of the ionizing radiation field. Ionization due to AGNs is ruled out as assessed with Chandra. As a consequence, we revisit the mass-metallicity relation using [NII]/Ha and a new calibration including [NII]/[SII] as recently introduced by Dopita et al. Consistent with our previous results, the most massive galaxies ($M_\ast\gtrsim10^{11}~M_\odot$) are fully enriched, while those at lower masses have metallicities lower than local galaxies. Finally, we demonstrate that the stellar masses, metallicities and star formation rates of the FMOS sample are well fit with a physically-motivated model for the chemical evolution of star-forming galaxies.

We present, for the first time, the local [CII] 158 um emission line luminosity function measured using a sample of more than 500 galaxies from the Revised Bright Galaxy Sample (RBGS). [CII] luminosities are measured from the Herschel PACS observations of the Luminous Infrared Galaxies in the Great Observatories All-sky LIRG Survey (GOALS) and estimated for the rest of the sample based on the far-IR luminosity and color. The sample covers 91.3% of the sky and is complete at S_60 um > 5.24 Jy. We calculated the completeness as a function of [CII] line luminosity and distance, based on the far-IR color and flux densities. The [CII] luminosity function is constrained in the range ~10^(7-9) (Lo) from both the 1/V_max and a maximum likelihood methods. The shape of our derived [CII] emission line luminosity function agrees well with the IR luminosity function. For the CO(1-0) and [CII] luminosity functions to agree, we propose a varying ratio of [CII]/CO(1-0) as a function of CO luminosity, with larger ratios for fainter CO luminosities. Limited [CII] high redshift observations as well as estimates based on the IR and UV luminosity functions are suggestive of an evolution in the [CII] luminosity function similar to the evolution trend of the cosmic star formation rate density. Deep surveys using ALMA with full capability will be able to confirm this prediction.

Highly accurate models of the galaxy population over cosmological volumes are necessary in order to predict the performance of upcoming cosmological missions. We present a data-driven model of the galaxy population constrained by deep 0.1-8 $\rm \mu m$ imaging and spectroscopic data in the COSMOS survey, with the immediate goal of simulating the spectroscopic redshift performance of the proposed SPHEREx mission. SPHEREx will obtain over the full-sky $R\sim41$ spectrophotometry at moderate spatial resolution ($\sim6"$) over the wavelength range 0.75-4.18 $\rm \mu m$ and $R\sim135$ over the wavelength range 4.18-5 $\rm \mu m$. We show that our simulation accurately reproduces a range of known galaxy properties, encapsulating the full complexity of the galaxy population and enables realistic, full end-to-end simulations to predict mission performance. Finally, we discuss potential applications of the simulation framework to future cosmology missions and give a description of released data products.

The offset of high redshift star-forming galaxies in the [OIII]/H$\beta$ versus [NII]/H$\alpha$ (O3N2) diagram in comparison with the local star-forming galaxy sequence is now well established. The physical origin of the shift is the subject of some debate, and has important implications for metallicity measurements based on strong lines at all redshifts. To investigate the origin of the O3N2 offset, we use a sample of ~100,000 star-forming galaxies from SDSS DR12 to probe the empirical correlations between emission line diagnostics and measurable galaxy physical properties. In particular, we examine how surface density of star formation, ionization parameter, nitrogen-to-oxygen (N/O) ratio, and stellar mass drive position in two key diagnostic diagrams: O3N2 and [OIII]/H$\beta$ versus [SII]/H$\alpha$ (O3S2). We show that, at fixed [OIII]/H$\beta$, galaxies falling closer to the high-redshift locus in O3N2 have higher star formation density, stellar mass and N/O ratios. We also find a tight correspondence in the distributions of stellar mass and N/O in the diagnostic diagrams. This relation, spanning a range of galaxy evolutionary states, suggests that the N/O-$M_{*}$ relation is more fundamental than the N/O-metallicity relation. We argue that a tight N/O-$M_{*}$ relation is well-motivated physically, and that the observed correlation of N/O with O/H in the local universe is primarily a reflection of the existence of the mass-metallicity relation. Because the mass-metallicity relation evolves much more rapidly with redshift than N/O-$M_{*}$, the N/O ratios of high redshift galaxies are significantly elevated in comparison with local galaxies with the same gas-phase O/H. The O3N2 shift and elevated N/O ratios observed in high redshift galaxies therefore come about as a natural consequence of the N/O-$M_{*}$ relation combined with the evolution of the mass-metallicity relation.

The current hierarchical merging paradigm and $\Lambda$CDM predict that the $z \sim 4-8$ universe should be a time in which the most massive galaxies are transitioning from their initial halo assembly to the later baryonic evolution seen in star-forming galaxies and quasars. However, no evidence of this transition has been found in many high redshift galaxy surveys including CFHTLS, CANDELS and SPLASH, the first studies to probe the high-mass end at these redshifts. Indeed, if halo mass to stellar mass ratios estimated at lower-redshift continue to $z \sim 6-8$, CANDELS and SPLASH report several orders of magnitude more $M \sim 10^{12-13} M_\odot$ halos than are possible to have formed by those redshifts, implying these massive galaxies formed impossibly early. We consider various systematics in the stellar synthesis models used to estimate physical parameters and possible galaxy formation scenarios in an effort to reconcile observation with theory. Although known uncertainties can greatly reduce the disparity between recent observations and cold dark matter merger simulations, even taking the most conservative view of the observations, there remains considerable tension with current theory.

We measure the H{\alpha} and [OIII] emission line properties as well as specific star-formation rates (sSFR) of spectroscopically confirmed 3<z<6 galaxies in COSMOS from their observed colors vs. redshift evolution. Our model describes consistently the ensemble of galaxies including intrinsic properties (age, metallicity, star-formation history), dust-attenuation, and optical emission lines. We forward-model the measured H{\alpha} equivalent-widths (EW) to obtain the sSFR out to z~6 without stellar mass fitting. We find a strongly increasing rest-frame H{\alpha} EW that is flattening off above z~2.5 with average EWs of 300-600A at z~6. The sSFR is increasing proportional to (1+z)^2.4 at z<2.2 and (1+z)^1.5 at higher redshifts, indicative of a fast mass build-up in high-z galaxies within e-folding times of 100-200Myr at z~6. The redshift evolution at z>3 cannot be fully explained in a picture of cold accretion driven growth. We find a progressively increasing [OIII]{\lambda}5007/H{\beta} ratio out to z~6, consistent with the ratios in local galaxies selected by increasing H{\alpha} EW (i.e., sSFR). This demonstrates the potential of using "local high-z analogs" to investigate the spectroscopic properties and relations of galaxies in the re-ionization epoch.

We measure a relation between the depth of four prominent rest-UV absorption complexes and metallicity for local galaxies and verify it up to z~3. We then apply this relation to a sample of 224 galaxies at 3.5 < z < 6.0 (<z> = 4.8) in COSMOS, for which unique UV spectra from DEIMOS and accurate stellar masses from SPLASH are available. The average galaxy population at z~5 and log(M/Msun) > 9 is characterized by 0.3-0.4 dex (in units of 12+log(O/H)) lower metallicities than at z~2, but comparable to z~3.5. We find galaxies with weak/no Ly-alpha emission to have metallicities comparable to z~2 galaxies and therefore may represent an evolved sub-population of z~5 galaxies. We find a correlation between metallicity and dust in good agreement with local galaxies and an inverse trend between metallicity and star-formation rate (SFR) consistent with observations at z~2. The relation between stellar mass and metallicity (MZ relation) is similar to z~3.5, however, there are indications of it being slightly shallower, in particular for the young, Ly-alpha emitting galaxies. We show that, within a "bathtub" approach, a shallower MZ relation is expected in the case of a fast (exponential) build-up of stellar mass with an e-folding time of 100-200 Myr. Due to this fast evolution, the process of dust production and metal enrichment as a function of mass could be more stochastic in the first billion years of galaxy formation compared to later times.

Calibrating the photometric redshifts of >10^9 galaxies for upcoming weak lensing cosmology experiments is a major challenge for the astrophysics community. The path to obtaining the required spectroscopic redshifts for training and calibration is daunting, given the anticipated depths of the surveys and the difficulty in obtaining secure redshifts for some faint galaxy populations. Here we present an analysis of the problem based on the self-organizing map, a method of mapping the distribution of data in a high-dimensional space and projecting it onto a lower-dimensional representation. We apply this method to existing photometric data from the COSMOS survey selected to approximate the anticipated Euclid weak lensing sample, enabling us to robustly map the empirical distribution of galaxies in the multidimensional color space defined by the expected Euclid filters. Mapping this multicolor distribution lets us determine where - in galaxy color space - redshifts from current spectroscopic surveys exist and where they are systematically missing. Crucially, the method lets us determine whether a spectroscopic training sample is representative of the full photometric space occupied by the galaxies in a survey. We explore optimal sampling techniques and estimate the additional spectroscopy needed to map out the color-redshift relation, finding that sampling the galaxy distribution in color space in a systematic way can efficiently meet the calibration requirements. While the analysis presented here focuses on the Euclid survey, similar analysis can be applied to other surveys facing the same calibration challenge, such as DES, LSST, and WFIRST.

In order to investigate the origin of quasars, we estimate the bias factor for low-luminosity quasars at high redshift for the first time. In this study, we use the two-point angular cross-correlation function (CCF) for both low-luminosity quasars at $-24<M_{\rm 1450}<-22$ and Lyman-break galaxies (LBGs). Our sample consists of both 25 low-luminosity quasars (16 objects are spectroscopically confirmed low-luminosity quasars) in the redshift range $3.1<z<4.5$ and 835 color-selected LBGs with $z^{\prime}_{\rm LBG}<25.0$ at $z\sim4$ in the COSMOS field. We have made our analysis for the following two quasar samples; (1) the spectroscopic sample (the 16 quasars confirmed by spectroscopy), and (2) the total sample (the 25 quasars including 9 quasars with photometric redshifts). The bias factor for low-luminosity quasars at $z\sim4$ is derived by utilizing the quasar-LBG CCF and the LBG auto-correlation function. We then obtain the $86\%$ upper limits of the bias factors for low-luminosity quasars, that are 5.63 and 10.50 for the total and the spectroscopic samples, respectively. These bias factors correspond to the typical dark matter halo masses, log $(M_{\rm DM}/(h^{-1}M_{\odot}))=$$12.7$ and $13.5$, respectively. This result is not inconsistent with the predicted bias for quasars which is estimated by the major merger models.

Using a grid of $\sim 2$ million elements ($\Delta z = 0.005$) adapted from COSMOS photometric redshift (photo-z) searches, we investigate the general properties of template-based photo-z likelihood surfaces. We find these surfaces are filled with numerous local minima and large degeneracies that generally confound rapid but "greedy" optimization schemes, even with additional stochastic sampling methods. In order to robustly and efficiently explore these surfaces, we develop BAD-Z [Brisk Annealing-Driven Redshifts (Z)], which combines ensemble Markov Chain Monte Carlo (MCMC) sampling with simulated annealing to sample arbitrarily large, pre-generated grids in approximately constant time. Using a mock catalog of 384,662 objects, we show BAD-Z samples $\sim 40$ times more efficiently compared to a brute-force counterpart while maintaining similar levels of accuracy. Our results represent first steps toward designing template-fitting photo-z approaches limited mainly by memory constraints rather than computation time.

We present a spectroscopic survey of galaxies in the COSMOS field using the Fiber Multi-Object Spectrograph (FMOS), a near-infrared instrument on the Subaru Telescope. Our survey is specifically designed to detect the Halpha emission line that falls within the H-band (1.6-1.8 um) spectroscopic window from star-forming galaxies with 1.4 < z < 1.7 and M_stellar>~10^10 Msolar. With the high multiplex capability of FMOS, it is now feasible to construct samples of over one thousand galaxies having spectroscopic redshifts at epochs that were previously challenging. The high-resolution mode (R~2600) effectively separates Halpha and [NII]6585 thus enabling studies of the gas-phase metallicity and photoionization state of the interstellar medium. The primary aim of our program is to establish how star formation depends on stellar mass and environment, both recognized as drivers of galaxy evolution at lower redshifts. In addition to the main galaxy sample, our target selection places priority on those detected in the far-infrared by Herschel/PACS to assess the level of obscured star formation and investigate, in detail, outliers from the star formation rate - stellar mass relation. Galaxies with Halpha detections are followed up with FMOS observations at shorter wavelengths using the J-long (1.11-1.35 um) grating to detect Hbeta and [OIII]5008 that provides an assessment of extinction required to measure star formation rates not hampered by dust, and an indication of embedded Active Galactic Nuclei. With 460 redshifts measured from 1153 spectra, we assess the performance of the instrument with respect to achieving our goals, discuss inherent biases in the sample, and detail the emission-line properties. Our higher-level data products, including catalogs and spectra, are available to the community.

SPHEREx (Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer) ( http://spherex.caltech.edu ) is a proposed all-sky spectroscopic survey satellite designed to address all three science goals in NASA's Astrophysics Division: probe the origin and destiny of our Universe; explore whether planets around other stars could harbor life; and explore the origin and evolution of galaxies. SPHEREx will scan a series of Linear Variable Filters systematically across the entire sky. The SPHEREx data set will contain R=40 spectra fir 0.75$<\lambda<$4.1$\mu$m and R=150 spectra for 4.1$<\lambda<$4.8$\mu$m for every 6.2 arc second pixel over the entire-sky. In this paper, we detail the extra-galactic and cosmological studies SPHEREx will enable and present detailed systematic effect evaluations. We also outline the Ice and Galaxy Evolution Investigations.

This report describes the 2014 study by the Science Definition Team (SDT) of the Wide-Field Infrared Survey Telescope (WFIRST) mission. It is a space observatory that will address the most compelling scientific problems in dark energy, exoplanets and general astrophysics using a 2.4-m telescope with a wide-field infrared instrument and an optical coronagraph. The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope as its top priority for a new large space mission. As conceived by the decadal survey, WFIRST would carry out a dark energy science program, a microlensing program to determine the demographics of exoplanets, and a general observing program utilizing its ultra wide field. In October 2012, NASA chartered a Science Definition Team (SDT) to produce, in collaboration with the WFIRST Study Office at GSFC and the Program Office at JPL, a Design Reference Mission (DRM) for an implementation of WFIRST using one of the 2.4-m, Hubble-quality telescope assemblies recently made available to NASA. This DRM builds on the work of the earlier WFIRST SDT, reported by Green et al. (2012) and the previous WFIRST-2.4 DRM, reported by Spergel et. (2013). The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the 1.3-m and 1.1-m designs considered previously, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The addition of an on-axis coronagraphic instrument to the baseline design enables imaging and spectroscopic studies of planets around nearby stars.

We used the PdBI to map a sample of 15 SMGs in the COSMOS field at the wavelength of 1.3 mm. The target SMGs were originally discovered in the JCMT/AzTEC 1.1 mm continuum survey at S/N=4-4.5. This paper presents, for the first time, interferometric millimetre-wavelength observations of these sources. The angular resolution of our observations, 1.8", allowed us to accurately determine the positions of the target SMGs. Using a detection threshold of S/N>4.5 regardless of multiwavelength counterpart association, and 4<S/N<=4.5 if a multiwavelength counterpart within 1.5" is also present, the total number of detections in our survey is 22. Three of our detected SMGs (AzTEC21, 27, and 28; which corresponds to 20%) are marginally resolved at our angular resolution, and these sources are found to have elongated or clumpy morphologies and/or multiple components. Using optical to NIR photometric redshifts, available spectroscopic redshifts, and redshifts estimated from the radio-to-submm spectral index we infer a median redshift of $\tilde{z}=3.20\pm0.25$ for our sample. To study the overall multiplicity and redshift distribution of flux-limited samples of SMGs we combined these sources with the 15 brightest AzTEC SMGs detected at 1.1 mm, AzTEC1-15, and studied previously. This constitutes a complete, flux- and S/N-limited 1.1-mm selected sample. We find that the median redshift for the 15 brightest AzTEC SMGs ($\tilde{z}=3.05\pm0.44$) is consistent with that for AzTEC16-30. This conforms to recent observational findings that SMGs do not exhibit any significant trend between the redshift and (sub)mm flux density. For the combined AzTEC1-30 sample we derive a median redshift of $\tilde{z}=3.17\pm0.27$, consistent with previous results based on mm-selected samples. We further infer that within the combined AzTEC1-30 sample $\sim25\pm9\%$ of sources separate into multiple components.

We study the physical properties of a sample of 6 SMGs in the COSMOS field, spectroscopically confirmed to lie at z>4. We use new GMRT 325 MHz and 3 GHz JVLA data to probe the rest-frame 1.4 GHz emission at z=4, and to estimate the sizes of the star-forming (SF) regions of these sources, resp. Combining our size estimates with those available in the literature for AzTEC1 and AzTEC3 we infer a median radio-emitting size for our z>4 SMGs of (0.63"+/-0.12")x(0.35"+/-0.05") or 4.1x2.3 kpc^2 (major times minor axis; assuming z=4.5) or lower if we take the two marginally resolved SMGs as unresolved. This is consistent with the sizes of SF regions in lower-redshift SMGs, and local normal galaxies, yet higher than the sizes of SF regions of local ULIRGs. Our SMG sample consists of a fair mix of compact and more clumpy systems with multiple, perhaps merging, components. With an average formation time of ~280 Myr, derived through modeling of the UV-IR SEDs, the studied SMGs are young systems. The average stellar mass, dust temperature, and IR luminosity we derive are M*~1.4x10^11 M_sun, T_dust~43 K, and L_IR~1.3x10^13L_sun, resp. The average L_IR is up to an order of magnitude higher than for SMGs at lower redshifts. Our SMGs follow the correlation between dust temperature and IR luminosity as derived for Herschel-selected 0.1<z<2 galaxies. We study the IR-radio correlation for our sources and find a deviation from that derived for z<3 ULIRGs (<q_IR>=1.95+/-0.26 for our sample, compared to q~2.6 for IR luminous galaxies at z<2). In summary, we find that the physical properties derived for our z>4 SMGs put them at the high end of the L_IR-T_dust distribution of SMGs, and that our SMGs form a morphologically heterogeneous sample. Thus, further in-depth analyses of large, statistical samples of high-redshift SMGs are needed to fully understand their role in galaxy formation and evolution.

Using the first 50% of data collected for the Spitzer Large Area Survey with Hyper-Suprime-Cam (SPLASH) observations on the 1.8 deg$^2$ Cosmological Evolution Survey (COSMOS) we estimate the masses and star formation rates of 3398 $M_*>10^{10}M_\odot $ star-forming galaxies at $4 < z < 6$ with a substantial population up to $M_* \gtrsim 10^{11.5} M_\odot$. We find that the strong correlation between stellar mass and star formation rate seen at lower redshift (the "main sequence" of star-forming galaxies) extends to $z\sim6$. The observed relation and scatter is consistent with a continued increase in star formation rate at fixed mass in line with extrapolations from lower-redshift observations. It is difficult to explain this continued correlation, especially for the most massive systems, unless the most massive galaxies are forming stars near their Eddington-limited rate from their first collapse. Furthermore, we find no evidence for moderate quenching at higher masses, indicating quenching either has not occurred prior to $z \sim 6$ or else occurs rapidly, so that few galaxies are visible in transition between star-forming and quenched.

We investigate the relationships between stellar mass, gas-phase oxygen abundance (metallicity), star formation rate, and dust content of star-forming galaxies at z$\sim$1.6 using Subaru/FMOS spectroscopy in the COSMOS field. The mass-metallicity relation at $z\sim1.6$ is steeper than the relation observed in the local Universe. The steeper MZ relation at $z\sim1.6$ is mainly due to evolution in the stellar mass where the MZ relation begins to turnover and flatten. This turnover mass is 1.2 dex larger at $z\sim1.6$. The most massive galaxies at $z\sim1.6$ ($\sim 10^{11}M_\odot$) are enriched to the level observed in massive galaxies in the local Universe. The mass-metallicity relation we measure at $z\sim1.6$ supports the suggestion of an empirical upper metallicity limit that does not significantly evolve with redshift. We find an anti-correlation between metallicity and star formation rate for galaxies at a fixed stellar mass at $z\sim1.6$ which is similar to trends observed in the local Universe. We do not find a relation between stellar mass, metallicity and star formation rate that is independent of redshift; our data suggest that there is redshift evolution in this relation. We examine the relation between stellar mass, metallicity and dust extinction. We find that at a fixed stellar mass dustier galaxies tend to be more metal rich. From examination of the stellar masses, metallicities, SFRs and dust extinctions we conclude that stellar mass is most closely related to dust extinction.

We present the first results from a near-IR spectroscopic survey of the COSMOS field, using the Fiber Multi-Object Spectrograph on the Subaru telescope, designed to characterize the star-forming galaxy population at $1.4<z<1.7$. The high-resolution mode is implemented to detect H$\alpha$ in emission between $1.6{\rm -}1.8 \mathrm{\mu m}$ with $f_{\rm H\alpha}\gtrsim4\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$. Here, we specifically focus on 271 sBzK-selected galaxies that yield a H$\alpha$ detection thus providing a redshift and emission line luminosity to establish the relation between star formation rate and stellar mass. With further $J$-band spectroscopy for 89 of these, the level of dust extinction is assessed by measuring the Balmer decrement using co-added spectra. We find that the extinction ($0.6\lesssim A_\mathrm{H\alpha} \lesssim 2.5$) rises with stellar mass and is elevated at high masses compared to low-redshift galaxies. Using this subset of the spectroscopic sample, we further find that the differential extinction between stellar and nebular emission \hbox{$E_\mathrm{star}(B-V)/E_\mathrm{neb}(B-V)$} is 0.7--0.8, dissimilar to that typically seen at low redshift. After correcting for extinction, we derive an H$\alpha$-based main sequence with a slope ($0.81\pm0.04$) and normalization similar to previous studies at these redshifts.

We present results from near-infrared spectroscopy of 26 emission-line galaxies at z ~ 2 obtained with the FIRE spectrometer on the Magellan Baade telescope. The sample was selected from the WISP survey, which uses the near-infrared grism of the Hubble Space Telescope Wide Field Camera 3 to detect emission-line galaxies over 0.3 < z < 2.3. Our FIRE follow-up spectroscopy (R~5000) over 1.0-2.5 micron permits detailed measurements of physical properties of the z~2 emission-line galaxies. Dust-corrected star formation rates for the sample range from ~5-100 M_sun yr-1. We derive a median metallicity for the sample of ~0.45 Z_sun, and the estimated stellar masses range from ~10^8.5 - 10^9.5 M_sun. The average ionization parameters measured for the sample are typically much higher than what is found for local star-forming galaxies. We derive composite spectra from the FIRE sample, from which we infer typical nebular electron densities of ~100-400 cm^-3. Based on the location of the galaxies and composite spectra on BPT diagrams, we do not find evidence for significant AGN activity in the sample. Most of the galaxies as well as the composites are offset in the BPT diagram toward higher [O III]/H-beta at a given [N II]/H-alpha, in agreement with other observations of z > 1 star-forming galaxies, but composite spectra derived from the sample do not show an appreciable offset from the local star-forming sequence on the [O III]/H-beta versus [S II]/H-alpha diagram. We infer a high nitrogen-to-oxygen abundance ratio from the composite spectrum, which may contribute to the offset of the high-redshift galaxies from the local star-forming sequence in the [O III]/H-beta versus [N II]/H-alpha diagram. We speculate that the elevated nitrogen abundance could result from substantial numbers of Wolf-Rayet stars in starbursting galaxies at z~2. (Abridged)

We present near-infrared emission line counts and luminosity functions from the HST WFC3 Infrared Spectroscopic Parallels (WISP) program for 29 fields (0.037 deg^2) observed using both the G102 and G141 grisms. Altogether we identify 1048 emission line galaxies with observed equivalent widths greater than 40 Angstroms, 467 of which have multiple detected emission lines. The WISP survey is sensitive to fainter flux levels (3-5x10^-17 ergs/s/cm^2) than the future space near-infrared grism missions aimed at baryonic acoustic oscillation cosmology (1-4x10^-16 ergs/s/cm^2), allowing us to probe the fainter emission line galaxies that the shallower future surveys may miss. Cumulative number counts of 0.7<z<1.5 galaxies reach 10,000 deg^-2 above an H-alpha flux of 2x10^-16 ergs/s/cm^2. H-alpha-emitting galaxies with comparable [OIII] flux are roughly 5 times less common than galaxies with just H-alpha emission at those flux levels. Galaxies with low H-alpha/[OIII] ratios are very rare at the brighter fluxes that future near-infrared grism surveys will probe; our survey finds no galaxies with H-alpha/[OIII] < 0.95 that have H-alpha flux greater than 3x10^-16 ergs/s/cm^2. Our H-alpha luminosity function contains a comparable number density of faint line emitters to that found by the NICMOS near-infrared grism surveys, but significantly fewer (factors of 3-4 less) high luminosity emitters. We also find that our high redshift (z=0.9-1.5) counts are in agreement with the high redshift (z=1.47) narrow band H-alpha survey of HiZELS (Sobral et al. 2013), while our lower redshift luminosity function (z=0.3-0.9) falls slightly below their z=0.84 result. The evolution in both the H-alpha luminosity function from z=0.3--1.5 and the [OIII] luminosity function from z=0.7-2.3 is almost entirely in the L* parameter, which steadily increases with redshift over those ranges.

We present the first robust measurement of the high redshift mass-metallicity (MZ) relation at 10^{8}< M/M_{\sun} < 10^{10}, obtained by stacking spectra of 83 emission-line galaxies with secure redshifts between 1.3 < z < 2.3. For these redshifts, infrared grism spectroscopy with the Hubble Space Telescope Wide Field Camera 3 is sensitive to the R23 metallicity diagnostic: ([OII]3726,3729 + [OIII] 4959,5007)/H\beta. Using spectra stacked in four mass quartiles, we find a MZ relation that declines significantly with decreasing mass, extending from 12+log(O/H) = 8.8 at M=10^{9.8} M_{\sun} to 12+log(O/H)= 8.2 at M=10^{8.2} M_{\sun}. After correcting for systematic offsets between metallicity indicators, we compare our MZ relation to measurements from the stacked spectra of galaxies with M>10^{9.5} M_{\sun} and z~2.3. Within the statistical uncertainties, our MZ relation agrees with the z~2.3 result, particularly since our somewhat higher metallicities (by around 0.1 dex) are qualitatively consistent with the lower mean redshift z=1.76 of our sample. For the masses probed by our data, the MZ relation shows a steep slope which is suggestive of feedback from energy-driven winds, and a cosmological downsizing evolution where high mass galaxies reach the local MZ relation at earlier times. In addition, we show that our sample falls on an extrapolation of the star-forming main sequence (the SFR-M_{*} relation) at this redshift. This result indicates that grism emission-line selected samples do not have preferentially high SFRs. Finally, we report no evidence for evolution of the mass-metallicity-SFR plane; our stack-averaged measurements show excellent agreement with the local relation.

We combine Hubble Space Telescope (HST) G102 & G141 NIR grism spectroscopy with HST/WFC3-UVIS, HST/WFC3-IR and Spitzer/IRAC[3.6\mu m] photometry to assemble a sample of massive (log(M_star/M_sun) ~ 11) and quenched galaxies at z~1.5. Our sample of 41 galaxies is the largest with G102+G141 NIR spectroscopy for quenched sources at these redshifts. In contrast to the local Universe, z~1.5 quenched galaxies in the high-mass range have a wide range of stellar population properties. We find their SEDs are well fitted with exponentially decreasing SFHs, and short star-formation time-scales (\tau<100Myr). Quenched galaxies also show a wide distribution in ages, between 1-4Gyr. In the (u-r)_0-versus-mass space quenched galaxies have a large spread in rest-frame color at a given mass. Most quenched galaxies populate the z~1.5 red-sequence (RS), but an important fraction of them (32%) have substantially bluer colors. Although with a large spread, we find that the quenched galaxies ON the RS have older median ages (3.1Gyr) than the quenched galaxies OFF the RS (1.5Gyr). We also show that a rejuvenated SED cannot reproduce the observed stacked spectra of (the bluer) quenched galaxies OFF the RS. We derive the upper limit on the fraction of massive galaxies ON the RS at z~1.5 to be <43%. We speculate that the young quenched galaxies OFF the RS are in a transition phase between vigorous star formation at z>2 and the z~1.5 RS. According to their estimated ages, the time required for quenched galaxies OFF the RS to join their counterparts ON the z~1.5 RS is of the order of ~1Gyr.

We present an investigation into how well the properties of the accretion flow onto a supermassive black hole may be coupled to those of the overlying hot corona. To do so, we specifically measure the characteristic spectral index, Gamma, of a power-law energy distribution, over an energy range of 2 to 10 keV, for X-ray selected, broad-lined radio-quiet AGN up to z~2 in COSMOS and E-CDF-S. We test the previously reported dependence between Gamma and black hole mass, FWHM and Eddington ratio using a sample of AGN covering a broad range in these parameters based on both the Mg ii and H-alpha emission lines with the later afforded by recent near infrared spectroscopic observations using Subaru/FMOS. We calculate the Eddington ratios, lambda_Edd, for sources where a bolometric luminosity (L_Bol) has been presented in the literature, based on SED fitting, or, for sources where these data do not exist, we calculate L_Bol using a bolometric correction to the LX, derived from a relationship between the bolometric correction, and LX/L3000. From a sample of 69 X-ray bright sources (>250 counts), where Gamma can be measured with greatest precision, with an estimate of L_Bol, we find a statistically significant correlation between Gamma and lambda_Edd, which is highly significant with a chance probability of 6.59x10^-8. A statistically significant correlation between Gamma and the FWHM of the optical lines is confirmed, but at lower significance than with lambda_Edd indicating that lambda_Edd is the key parameter driving conditions in the corona. Linear regression analysis reveals that Gamma=(0.32+/-0.05)log10 lambda_Edd+(2.27+/-0.06) and Gamma=(-0.69+/-0.11)log10(FWHM/km/s)+(4.44+/-0.42). Our results on Gamma-lambda_Edd are in very good agreement with previous results. (ABRIDGED)

Large-scale structures (LSS) out to z $< 3.0$ are measured in the Cosmic Evolution Survey (COSMOS) using extremely accurate photometric redshifts (photoz). The Ks-band selected sample (from Ultra-Vista) is comprised of 155,954 galaxies. Two techniques -- adaptive smoothing and Voronoi tessellation -- are used to estimate the environmental densities within 127 redshift slices. Approximately 250 statistically significant overdense structures are identified out to z $= 3.0$ with shapes varying from elongated filamentary structures to more circularly symmetric concentrations. We also compare the densities derived for COSMOS with those based on semi-analytic predictions for a $\Lambda$CDM simulation and find excellent overall agreement between the mean densities as a function of redshift and the range of densities. The galaxy properties (stellar mass, spectral energy distributions (SEDs) and star formation rates (SFRs)) are strongly correlated with environmental density and redshift, particularly at z $< 1.0 - 1.2$. Classifying the spectral type of each galaxy using the rest-frame b-i color (from the photoz SED fitting), we find a strong correlation of early type galaxies (E-Sa) with high density environments, while the degree of environmental segregation varies systematically with redshift out to z $\sim 1.3$. In the highest density regions, 80% of the galaxies are early types at z=0.2 compared to only 20% at z = 1.5. The SFRs and the star formation timescales exhibit clear environmental correlations. At z $> 0.8$, the star formation rate density (SFRD) is uniformly distributed over all environmental density percentiles, while at lower redshifts the dominant contribution is shifted to galaxies in lower density environments.

We present an analysis of broad emission lines observed in moderate-luminosity active galactic nuclei (AGNs), typical of those found in X-ray surveys of deep fields, with the aim to test the validity of single-epoch virial black hole mass estimates. We have acquired near-infrared (NIR) spectra of AGNs up to z ~ 1.8 in the COSMOS and Extended Chandra Deep Field-South Survey, with the Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope. These low-resolution NIR spectra provide a significant detection of the broad Halpha line that has been shown to be a reliable probe of black hole mass at low redshift. Our sample has existing optical spectroscopy which provides a detection of MgII, a broad emission line typically used for black hole mass estimation at z > 1. We carry out a spectral-line fitting procedure using both Halpha and MgII to determine the virial velocity of gas in the broad line region, the monochromatic continuum luminosity at 3000 A, and the total Halpha line luminosity. With a sample of 43 AGNs spanning a range of two decades in luminosity (i.e., L ~ 10^44-46 ergs/s), we find a tight correlation between the continuum and line luminosity with a distribution characterized by <log(L_3000/L_Halpha)> = 1.52 and a dispersion sigma = 0.16. There is also a close one-to-one relationship between the FWHM of Halpha and of MgII up to 10000 km/s with a dispersion of 0.14 in the distribution of the logarithm of their ratios. Both of these then lead to there being very good agreement between Halpha- and MgII-based masses over a wide range in black hole mass (i.e., M_BH ~ 10^7-9 M_sun). We do find a small offset in MgII-based masses, relative to those based on Halpha, of +0.17 dex and a dispersion sigma = 0.32. In general, these results demonstrate that local scaling relations, using MgII or Halpha, are applicable for AGN at moderate luminosities and up to z ~ 2.

Spectroscopic observations of Halpha and Hbeta emission lines of 128 star-forming galaxies in the redshift range 0.75<z<1.5 are presented. These data were taken with slitless spectroscopy using the G102 and G141 grisms of the Wide-Field-Camera 3 (WFC3) on board the Hubble Space Telescope as part of the WFC3 Infrared Spectroscopic Parallel (WISP) survey. Interstellar dust extinction is measured from stacked spectra that cover the Balmer decrement (Halpha/Hbeta). We present dust extinction as a function of Halpha luminosity (down to 3 x 10^{41} erg/s), galaxy stellar mass (reaching 4 x 10^{8} Msun), and rest-frame Halpha equivalent width. The faintest galaxies are two times fainter in Halpha luminosity than galaxies previously studied at z~1.5. An evolution is observed where galaxies of the same Halpha luminosity have lower extinction at higher redshifts, whereas no evolution is found within our error bars with stellar mass. The lower Halpha luminosity galaxies in our sample are found to be consistent with no dust extinction. We find an anti-correlation of the [OIII]5007/Halpha flux ratio as a function of luminosity where galaxies with L_{Halpha}<5 x 10^{41} erg/s are brighter in [OIII]5007 than Halpha. This trend is evident even after extinction correction, suggesting that the increased [OIII]5007/Halpha ratio in low luminosity galaxies is likely due to lower metallicity and/or higher ionization parameters.

We report on the discovery of a Type 1 quasar, SDSS 0956+5128, with a surprising combination of extreme velocity offsets. SDSS 0956+5128 is a broad-lined quasar exhibiting emission lines at three substantially different redshifts: a systemic redshift of z ~ 0.714 based on narrow emission lines, a broad MgII emission line centered 1200 km/s bluer than the systemic velocity, at z ~ 0.707, and broad H\alpha and H\beta emission lines centered at z ~ 0.690. The Balmer line peaks are 4100 km/s bluer than the systemic redshift. There are no previously known objects with such an extreme difference between broad MgII and broad Balmer emission. The two most promising explanations are either an extreme disk emitter or a high-velocity black hole recoil. However, neither explanation appears able to explain all of the observed features of SDSS 0956+5128, so the object may provide a challenge to our general understanding of quasar physics.

We investigate spectroscopically measured Ly{\alpha} equivalent widths and escape fractions of 244 sources of which 95 are Lyman Break Galaxies (LBGs) and 106 Lyman Alpha Emitters (LAEs) at z~4.2, z~4.8, and z~5.6 selected from intermediate and narrow band observations. The sources were selected from the Cosmic Evolution Survey (COSMOS), and observed with the DEIMOS spectrograph. We find that the distribution of equivalent widths shows no evolution with redshift for both the LBG selected sources and the intermediate/narrowband LAEs. We also find that the Ly{\alpha} escape fraction of intermediate and narrow band LAEs is on average higher and has a larger variation than the escape fraction of LBG selected sources. The escape fraction does not show a dependence with redshift. Similar to what has been found for LAEs at low redshifts, the sources with the highest extinctions show the lowest escape fractions. The range of escape fractions increases with decreasing extinction. This is evidence that the dust extinction is the most important factor affecting the escape of Ly{\alpha} photons, but at low extinctions other factors such as HI covering fraction and gas kinematics can be just as effective at inhibiting the escape of Ly{\alpha} photons.

We investigate the high-redshift quasar luminosity function (QLF) down to an apparent magnitude of I(AB) = 25 in the Cosmic Evolution Survey (COSMOS). Careful analysis of the extensive COSMOS photometry and imaging data allows us to identify and remove stellar and low-redshift contaminants, enabling a selection that is nearly complete for type-1 quasars at the redshifts of interest. We find 155 likely quasars at z > 3.1, 39 of which have prior spectroscopic confirmation. We present our sample in detail and use these confirmed and likely quasars to compute the rest-frame UV QLF in the redshift bins 3.1 < z < 3.5 and 3.5 < z < 5. The space density of faint quasars decreases by roughly a factor of four from z \sim 3.2 to z \sim 4, with faint-end slopes of {\beta} \sim -1.7 at both redshifts. The decline in space density of faint optical quasars at z > 3 is similar to what has been found for more luminous optical and X-ray quasars. We compare the rest-frame UV luminosity functions found here with the X-ray luminosity function at z > 3, and find that they evolve similarly between z \sim 3.2 and z \sim 4; however, the different normalizations imply that roughly 75% of X-ray bright active galactic nuclei (AGN) at z \sim 3 - 4 are optically obscured. This fraction is higher than found at lower redshift and may imply that the obscured, type-2 fraction continues to increase with redshift at least to z \sim 4. Finally, the implications of the results derived here for the contribution of quasars to cosmic reionization are discussed.

We present the result of our low-luminosity quasar survey in the redshift range of 4.5 < z < 5.5 in the COSMOS field. Using the COSMOS photometric catalog, we selected 15 quasar candidates with 22 < i' < 24 at z~5, that are ~ 3 mag fainter than the SDSS quasars in the same redshift range. We obtained optical spectra for 14 of the 15 candidates using FOCAS on the Subaru Telescope and did not identify any low-luminosity type-1 quasars at z~5 while a low-luminosity type-2 quasar at z~5.07 was discovered. In order to constrain the faint end of the quasar luminosity function at z~5, we calculated the 1sigma confidence upper limits of the space density of type-1 quasars. As a result, the 1sigma confidence upper limits on the quasar space density are Phi< 1.33*10^{-7} Mpc^{-3} mag^{-1} for -24.52 < M_{1450} < -23.52 and Phi< 2.88*10^{-7} Mpc^{-3} mag^{-1} for -23.52 < M_{1450} < -22.52. The inferred 1sigma confidence upper limits of the space density are then used to provide constrains on the faint-end slope and the break absolute magnitude of the quasar luminosity function at z~5. We find that the quasar space density decreases gradually as a function of redshift at low luminosity (M_{1450} ~ -23), being similar to the trend found for quasars with high luminosity (M_{1450}<-26). This result is consistent with the so-called downsizing evolution of quasars seen at lower redshifts.

The Chandra COSMOS Survey (C-COSMOS) is a large, 1.8 Ms, Chandra program that has imaged the central 0.9 deg^2 of the COSMOS field down to limiting depths of 1.9 10^-16 erg cm^-2 s-1 in the 0.5-2 keV band, 7.3 10^-16 erg cm^-2 s^-1 in the 2-10 keV band, and 5.7 10^-16 erg cm^-2 s-1 in the 0.5-10 keV band. In this paper we report the i, K and 3.6micron identifications of the 1761 X-ray point sources. We use the likelihood ratio technique to derive the association of optical/infrared counterparts for 97% of the X-ray sources. For most of the remaining 3%, the presence of multiple counterparts or the faintness of the possible counterpart prevented a unique association. For only 10 X-ray sources we were not able to associate a counterpart, mostly due to the presence of a very bright field source close by. Only 2 sources are truly empty fields. Making use of the large number of X-ray sources, we update the "classic locus" of AGN and define a new locus containing 90% of the AGN in the survey with full band luminosity >10^42 erg/s. We present the linear fit between the total i band magnitude and the X-ray flux in the soft and hard band, drawn over 2 orders of magnitude in X-ray flux, obtained using the combined C-COSMOS and XMM-COSMOS samples. We focus on the X-ray to optical flux ratio (X/O) and we test its known correlation with redshift and luminosity, and a recently introduced anti-correlation with the concentration index (C). We find a strong anti-correlation (though the dispersion is of the order of 0.5 dex) between C and X/O, computed in the hard band, and that 90% of the obscured AGN in the sample with morphological information live in galaxies with regular morphology (bulgy and disky/spiral), suggesting that secular processes govern a significant fraction of the BH growth at X-ray luminosities of 10^43- 10^44.5 erg/s.

We present the discovery of three late type (>T4) brown dwarfs, including a probable Y dwarf, in the WFC3 Infrared Spectroscopic Parallels (WISP) Survey. We use the G141 grism spectra to determine the spectral types of the dwarfs and derive distance estimates based on a comparison with nearby T dwarfs with known parallaxes. These are the most distant spectroscopically confirmed T/Y dwarfs, with the farthest at an estimated distance of ~400 pc. We compare the number of cold dwarfs found in the WISP survey with simulations of the brown dwarf mass function. The number found is generally consistent with an initial stellar mass function dN/dM \propto M^{-\alpha} with \alpha = 0.0--0.5, although the identification of a Y dwarf is somewhat surprising and may be indicative of either a flatter absolute magnitude/spectral type relation than previously reported or an upturn in the number of very late type brown dwarfs in the observed volume.

We present the high redshift (3<z<5.3) 0.5-2 keV number counts and the 2-10 keV (rest frame) space density of X-ray selected AGNs detected in the Chandra COSMOS survey. The sample comprises 81 X-ray detected sources with available spectroscopic (31) and photometric (50) redshifts plus 20 sources with a formal z_phot<3 but with a broad photometric redshift probability distribution, such that z_phot+1sigma>3. 81 sources are selected in the 0.5-2 keV band, 14 are selected in the 2-10 keV and 6 in the 0.5-10 keV bands. We sample the high luminosity (log L_(2-10 keV)>44.15 erg/s) space density up to z~5 and a fainter luminosity range (43.5<log L(2-10 keV)<44.15 erg/s) than previous studies, up to z=3.5. We weighted the contribution to the number counts and the space density of the sources with photometric redshift by using their probability of being at z>3. We find that the space density of high-luminosity AGNs declines exponentially at all the redshifts, confirming the trend observed for optically selected quasars. At lower luminosity, the measured space density is not conclusive, and a larger sample of faint sources is needed. Comparisons with optical luminosity functions and BH formation models are presented together with prospects for future surveys.

With this paper, we release accurate photometric redshifts for 1692 counterparts to Chandra sources in the central square degree of the COSMOS field. The availability of a large training set of spectroscopic redshifts that extends to faint magnitudes enabled photometric redshifts comparable to the highest quality results presently available for normal galaxies. We demonstrate that morphologically extended, faint X-ray sources without optical variability are more accurately described by a library of normal galaxies (corrected for emission lines) than by AGN-dominated templates, even if these sources have AGN-like X-ray luminosities. Preselecting the library on the bases of the source properties allowed us to reach an accuracy sigma_(Delta z/(1+z_spec)) \sim0.015 with a fraction of outliers of 5.8% for the entire Chandra-COSMOS sample. In addition, we release revised photometric redshifts for the 1735 optical counterparts of the XMM-detected sources over the entire 2 sq. deg.of COSMOS. For 248 sources, our updated photometric redshift differs from the previous release by Delta z>0.2. These changes are predominantly due to the inclusion of newly available deep H-band photometry H_AB=24 mag. We illustrate once again the importance of a spectroscopic training sample and how an assumption about the nature of a source together with the number and the depth of the available bands influence the accuracy of the photometric redshifts determined for AGN. These considerations should be kept in mind when defining the observational strategies of upcoming large surveys targeting AGN, such as eROSITA at X-ray energies and ASKAP/EMU in the radio band.

We present ground-based spectroscopic verification of six Y dwarfs (see Cushing et al), eighty-nine T dwarfs, eight L dwarfs, and one M dwarf identified by the Wide-field Infrared Survey Explorer (WISE). Eighty of these are cold brown dwarfs with spectral types greater than or equal to T6, six of which have been announced earlier in Mainzer et al and Burgasser et al. We present color-color and color-type diagrams showing the locus of M, L, T, and Y dwarfs in WISE color space. Near-infrared classifications as late as early Y are presented and objects with peculiar spectra are discussed. After deriving an absolute WISE 4.6 um (W2) magnitude vs. spectral type relation, we estimate spectrophotometric distances to our discoveries. We also use available astrometric measurements to provide preliminary trigonometric parallaxes to four our discoveries, which have types of L9 pec (red), T8, T9, and Y0; all of these lie within 10 pc of the Sun. The Y0 dwarf, WISE 1541-2250, is the closest at 2.8 (+1.3,-0.6) pc; if this 2.8 pc value persists after continued monitoring, WISE 1541-2250 will become the seventh closest stellar system to the Sun. Another ten objects, with types between T6 and >Y0, have spectrophotometric distance estimates also placing them within 10 pc. The closest of these, the T6 dwarf WISE 1506+7027, is believed to fall at a distance of roughly 4.9 pc. WISE multi-epoch positions supplemented with positional info primarily from Spitzer/IRAC allow us to calculate proper motions and tangential velocities for roughly one half of the new discoveries. This work represents the first step by WISE to complete a full-sky, volume-limited census of late-T and Y dwarfs. Using early results from this census, we present preliminary, lower limits to the space density of these objects and discuss constraints on both the functional form of the mass function and the low-mass limit of star formation.

Based on broad/narrow-band photometry and Keck DEIMOS spectroscopy we report a redshift of z=4.64-0.08+0.06 for AzTEC/COSMOS 1, the brightest sub-mm galaxy in the AzTEC/COSMOS field. In addition to the COSMOS-survey X-ray to radio data, we report observations of the source with Herschel/PACS (100, 160 micron), CSO/SHARC II (350 micron), CARMA and PdBI (3 mm). We do not detect CO(5-4) line emission in the covered redshift ranges, 4.56-4.76 (PdBI/CARMA) and 4.94-5.02 (CARMA). If the line is within this bandwidth, this sets 3sigma upper limits on the gas mass to <~8x10^9 M_Sol and <~5x10^10 M_Sol, respectively (assuming similar conditions as observed in z~2 SMGs). This could be explained by a low CO-excitation in the source. Our analysis of the UV-IR spectral energy distribution of AzTEC 1 shows that it is an extremely young (<~50 Myr), massive (M*~10^11 M_Sol), but compact (<~2 kpc) galaxy forming stars at a rate of ~1300 M_Sol/yr. Our results imply that AzTEC 1 is forming stars in a 'gravitationally bound' regime in which gravity prohibits the formation of a superwind, leading to matter accumulation within the galaxy and further generations of star formation.

We present an interactive IDL program for viewing and analyzing astronomical spectra in the context of modern imaging surveys. SpecPro's interactive design lets the user simultaneously view spectroscopic, photometric, and imaging data, allowing for rapid object classification and redshift determination. The spectroscopic redshift can be determined with automated cross-correlation against a variety of spectral templates or by overlaying common emission and absorption features on the 1-D and 2-D spectra. Stamp images as well as the spectral energy distribution (SED) of a source can be displayed with the interface, with the positions of prominent photometric features indicated on the SED plot. Results can be saved to file from within the interface. In this paper we discuss key program features and provide an overview of the required data formats.