Cold dark matter haloes consist of a relatively smooth dark matter component as well as a system of bound subhaloes. It is the prevailing practice to include all halo mass, including mass in subhaloes, in studies of halo density profiles. However, often in observational studies satellites are treated as having their own distinct dark matter density profiles in addition to the profile of the host. This difference makes comparisons between theoretical and observed results difficult. In this work we investigate density profiles of the smooth components of host haloes by excluding mass contained within subhaloes. We find that the density profiles of the smooth halo component (without subhaloes) differs substantially from the conventional halo density profile. Smooth profiles decline more rapidly at large radii and are not well characterised by the standard NFW profile. We also find that concentrations derived from smooth density profiles exhibit less scatter at fixed mass and a weaker mass dependence than standard concentrations. Both smooth and standard halo profiles can be described by a generalised Einasto profile, an Einasto profile with a modified central slope, with smaller residuals than either an NFW or Einasto profile. These results hold for both Milky Way-mass and cluster-mass haloes. This new characterisation of smooth halo profiles can be useful for many analyses, such as lensing and dark matter annihilation, in which the smooth and clumpy components of a halo should be accounted for separately.

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)

Using the CANDELS photometric catalogs for the HST/ACS and WFC3, we identified massive evolved galaxies at $3 < z < 4.5$, employing three different selection methods. We find the comoving number density of these objects to be $\sim 2 \times 10^{-5}$ and $8 \times 10^{-6}Mpc^{-3}$ after correction for completeness for two redshift bins centered at $z=3.4, 4.7$. We quantify a measure of how much confidence we should have for each candidate galaxy from different selections and what are the conservative error estimates propagated into our selection. Then we compare the evolution of the corresponding number densities and their stellar mass density with numerical simulations, semi-analytical models, and previous observational estimates, which shows slight tension at higher redshifts as the models tend to underestimate the number and mass densities. By estimating the average halo masses of the candidates ($M_h \approx 4.2, 1.9, 1.3 \times 10^{12} M_\odot$ for redshift bins centered at $z=3.4, 4.1, 4.7$), we find them to be consistent with halos that were efficient in turning baryons to stars and were relatively immune to the feedback effects and on the verge of transition into hot-mode accretion. This can suggest the relative cosmological starvation of the cold gas followed by an overconsumption phase in which the galaxy consumes the available cold gas rapidly as one of the possible drivers for the quenching of the massive evolved population at high redshift.

A photometric redshift sample of Luminous Red Galaxies (hereafter LRGs) obtained from The DECam Legacy Survey (DECaLS) is analysed to probe cosmic distances by exploiting the wedge approach of the two-point correlation function. Although the cosmological information is highly contaminated by the uncertainties existing in the photometric redshifts from the galaxy map, an angular diameter distance can be probed at the perpendicular configuration in which the measured correlation function is minimally contaminated. An ensemble of wedged correlation functions selected up to a given threshold based on having the least contamination was studied in the previous work (arXiv:1903.09651v2 [astro-ph.CO]) using simulations, and the extracted cosmological information was unbiased within this threshold. We apply the same methodology for analysing the LRG sample from DECaLS which will provide the optical imaging for targeting two-thirds of the DESI footprint and measure the angular diameter distances at $z=0.69$ and $z=0.87$ to be $D_{A}(0.697)=(1499 \pm 77\,\mathrm{Mpc})(r_{d}/r_{d,fid})$ and $D_{A}(0.874)=(1680 \pm 109\,\mathrm{Mpc})(r_{d}/r_{d,fid})$ with a fractional error of 5.14% and 6.48% respectively. We obtain a value of $H_{0}=67.59\pm5.52$ km/s/Mpc which supports the $H_0$ measured by all other BAO results and is consistent with $\Lambda$CDM model.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

The 1.5 Seyfert galaxy NGC 3227 has been observed by several X-ray missions. We carried out combined analysis of the data obtained by more recent major observations of this source - two observations performed by XMM-Newton in 2000 and 2006 and six observations performed by Suzaku in 2008. A unified model was constructed which is consistent with all eight of the observations by the two satellites with large intensity and spectral changes. The model consists of a hard power law with the spectral index of 1.35-1.65 which is interpreted as the Comptonized emission from the corona above an accretion disk. In the high flux states an additional soft excess component dominates, which is consistent with a model with either a steeper power law with an index of 3.5-3.6 or the warm Comptonization component. These emissions from the central engine are absorbed by a neutral partial covering material and warm absorbers. A reflection component and several emission lines are also present. We examined the relationship between the intrinsic luminosity and the absorbers' physical parameters such as the column density, which suggests that the source expanded significantly during the bright states where the soft excess is greatly enhanced.

The Milky Way provides an ideal laboratory to test our understanding of galaxy evolution, owing to our ability to observe our Galaxy over fine scales. However, connecting the Galaxy to the wider galaxy population remains difficult, due to the challenges posed by our internal perspective and to the different observational techniques employed. Here, we present a sample of galaxies identified as Milky Way Analogs (MWAs) on the basis of their stellar masses and bulge-to-total ratios, observed as part of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. We analyse the galaxies in terms of their stellar kinematics and populations as well as their ionised gas contents. We find our sample to contain generally young stellar populations in their outskirts. However, we find a wide range of stellar ages in their central regions, and we detect central AGN-like or composite-like activity in roughly half of the sample galaxies, with the other half consisting of galaxies with central star-forming emission or emission consistent with old stars. We measure gradients in gas metallicity and stellar metallicity that are generally flatter in physical units than those measured for the Milky Way; however, we find far better agreement with the Milky Way when scaling gradients by galaxies' disc scale lengths. From this, we argue much of the discrepancy in metallicity gradients to be due to the relative compactness of the Milky Way, with differences in observing perspective also likely to be a factor.

We present measurements of the redshift-dependent clustering of a DESI-like luminous red galaxy (LRG) sample selected from the Legacy Survey imaging dataset, and use the halo occupation distribution (HOD) framework to fit the clustering signal. The LRG sample contains 2.7 million objects over the redshift range of $0.4 < z < 0.9$ over 5655 sq. degrees. We have developed new photometric redshift (photo-$z$) estimates using the Legacy Survey DECam and WISE photometry, with $\sigma_{\mathrm{NMAD}} = 0.02$ precision for LRGs. We compute the projected correlation function using new methods that maximize signal-to-noise while incorporating redshift uncertainties. We present a novel algorithm for dividing irregular survey geometries into equal-area patches for jackknife resampling. For a 5-parameter HOD model fit using the MultiDark halo catalog, we find that there is little evolution in HOD parameters except at the highest-redshifts. The inferred large-scale structure bias is largely consistent with constant clustering amplitude over time. In an appendix, we explore limitations of MCMC fitting using stochastic likelihood estimates resulting from applying HOD methods to N-body catalogs, and present a new technique for finding best-fit parameters in this situation. Accompanying this paper we have released the PRLS (Photometric Redshifts for the Legacy Surveys) catalog of photo-$z$'s obtained by applying the methods used in this work to the full Legacy Survey Data Release 8 dataset. This catalog provides accurate photometric redshifts for objects with $z < 21$ over more than 16,000 square degrees of sky.

Many scientific investigations of photometric galaxy surveys require redshift estimates, whose uncertainty properties are best encapsulated by photometric redshift (photo-z) posterior probability density functions (PDFs). A plethora of photo-z PDF estimation methodologies abound, producing discrepant results with no consensus on a preferred approach. We present the results of a comprehensive experiment comparing twelve photo-z algorithms applied to mock data produced for the Large Synoptic Survey Telescope (LSST) Dark Energy Science Collaboration (DESC). By supplying perfect prior information, in the form of the complete template library and a representative training set as inputs to each code, we demonstrate the impact of the assumptions underlying each technique on the output photo-z PDFs. In the absence of a notion of true, unbiased photo-z PDFs, we evaluate and interpret multiple metrics of the ensemble properties of the derived photo-z PDFs as well as traditional reductions to photo-z point estimates. We report systematic biases and overall over/under-breadth of the photo-z PDFs of many popular codes, which may indicate avenues for improvement in the algorithms or implementations. Furthermore, we raise attention to the limitations of established metrics for assessing photo-z PDF accuracy; though we identify the conditional density estimate (CDE) loss as a promising metric of photo-z PDF performance in the case where true redshifts are available but true photo-z PDFs are not, we emphasize the need for science-specific performance metrics.

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $\kappa < 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ \chi < 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv.

We analyze the clustering of galaxies in the first public data release of the HSC Subaru Strategic Program. Despite the relatively small footprints of the observed fields, the data are an excellent proxy for the deep photometric datasets that will be acquired by LSST, and are therefore an ideal test bed for the analysis methods being implemented by the LSST DESC. We select a magnitude limited sample with $i<24.5$ and analyze it in four redshift bins covering $0.15\lesssim z \lesssim1.5$. We carry out a Fourier-space analysis of the two-point clustering of this sample, including all auto- and cross-correlations. We demonstrate the use of map-level deprojection methods to account for fluctuations in the galaxy number density caused by observational systematics. Through an HOD analysis, we place constraints on the characteristic halo masses of this sample, finding a good fit up to scales $k_{\rm max}=1\,{\rm Mpc}^{-1}$, including both auto- and cross-correlations. Our results show monotonically decreasing average halo masses, which can be interpreted in terms of the drop-out of red galaxies at high redshifts for a flux-limited sample. In terms of photometric redshift systematics, we show that additional care is needed in order to marginalize over uncertainties in the redshift distribution in galaxy clustering, and that these uncertainties can be constrained by including cross-correlations. We are able to make a $\sim3\sigma$ detection of lensing magnification in the HSC data. Our results are stable to variations in $\sigma_8$ and $\Omega_c$ and we find constraints that agree well with measurements from Planck and low-redshift probes. Finally, we use our pipeline to study the clustering of galaxies as a function of limiting flux, and provide a simple fitting function for the linear galaxy bias for magnitude limited samples as a function of limiting magnitude and redshift. [abridged]

We present a robust method, weighted von Mises kernel density estimation, along with boundary correction to reconstruct the underlying number density field of galaxies. We apply this method to galaxies brighter than $\rm HST/F160w\le 26$ AB mag at the redshift range of $0.4\leq z \leq 5$ in the five CANDELS fields (GOODS-N, GOODS-S, EGS, UDS, and COSMOS). We then use these measurements to explore the environmental dependence of the star formation activity of galaxies. We find strong evidence of environmental quenching for massive galaxies ($\rm M \gtrsim 10^{11} \rm {M}_\odot$) out to $z\sim 3.5$ such that an over-dense environment hosts $\gtrsim 20\%$ more massive quiescent galaxies compared to an under-dense region. We also find that environmental quenching efficiency grows with stellar mass and reaches $\sim 60\%$ for massive galaxies at $z\sim 0.5$. The environmental quenching is also more efficient in comparison to the stellar mass quenching for low mass galaxies ($\rm M \lesssim 10^{10} \rm {M}_\odot$) at low and intermediate redshifts ($z\lesssim 1.2$). Our findings concur thoroughly with the "over-consumption" quenching model where the termination of cool gas accretion (cosmological starvation) happens in an over-dense environment and the galaxy starts to consume its remaining gas reservoir in depletion time. The depletion time depends on the stellar mass and could explain the evolution of environmental quenching efficiency with the stellar mass.

We present results from a search for a radio transient associated with the LIGO/Virgo source S190814bv, a likely neutron star-black hole (NSBH) merger, with the Australian Square Kilometre Array Pathfinder. We imaged a $30\,{\rm deg}^2$ field at $\Delta T$=2, 9 and 33 days post-merger at a frequency of 944\,MHz, comparing them to reference images from the Rapid ASKAP Continuum Survey observed 110 days prior to the event. Each epoch of our observations covers $89\%$ of the LIGO/Virgo localisation region. We conducted an untargeted search for radio transients in this field, resulting in 21 candidates. For one of these, \object[AT2019osy]{AT2019osy}, we performed multi-wavelength follow-up and ultimately ruled out the association with S190814bv. All other candidates are likely unrelated variables, but we cannot conclusively rule them out. We discuss our results in the context of model predictions for radio emission from neutron star-black hole mergers and place constrains on the circum-merger density and inclination angle of the merger. This survey is simultaneously the first large-scale radio follow-up of an NSBH merger, and the most sensitive widefield radio transients search to-date.

Joint survey processing (JSP) is the pixel level combination of LSST, Euclid, and WFIRST datasets. By combining the high spatial resolution of the space-based datasets with deep, seeing-limited, ground-based images in the optical bands, systematics like source confusion and astrometric mismatch can be addressed to derive the highest precision optical/infrared photometric catalogs. This white paper highlights the scientific motivation, computational and algorithmic needs to build joint pixel level processing capabilities, which the individual projects by themselves will not be able to support. Through this white paper, we request that the Astro2020 decadal committee recognize the JSP effort as a multi-agency project with the natural outcome being a collaborative effort among groups which are normally supported by a single agency. JSP will allow the U.S. (and international) astronomical community to manipulate the flagship data sets and undertake innovative science investigations ranging from solar system object characterization, exoplanet detections, nearby galaxy rotation rates and dark matter properties, to epoch of reionization studies. It will also result in the ultimate constraints on cosmological parameters and the nature of dark energy, with far smaller uncertainties and a better handle on systematics than by any one survey alone.

The study of the space density of bright AGNs at $z>4$ has been subject to extensive effort given its importance for the estimate of the cosmological ionizing emissivity and growth of supermassive black holes. In this context we have recently derived high space densities of AGNs at $z\sim 4$ and $-25<M_{1450}<-23$ in the COSMOS field from a spectroscopically complete sample. In the present paper we attempt to extend the knowledge of the AGN space density at fainter magnitudes ($-22.5<M_{1450}<-18.5$) in the $4<z<6.1$ redshift interval by means of a multiwavelength sample of galaxies in the CANDELS GOODS-South, GOODS-North and EGS fields. We use an updated criterion to extract faint AGNs from a population of NIR (rest-frame UV) selected galaxies at photometric $z>4$ showing X-ray detection in deep Chandra images available for the three CANDELS fields. We have collected a photometric sample of 32 AGN candidates in the selected redshift interval, six of which having spectroscopic redshifts. Including our COSMOS sample as well as other bright QSO samples allows a first guess on the shape of the UV luminosity function at $z\sim 4.5$. The resulting emissivity and photoionization rate appear consistent with that derived from the photoionization level of the intergalactic medium at $z\sim 4.5$. An extrapolation to $z\sim 5.6$ suggests an important AGN contribution to the IGM ionization if there are no significant changes in the shape of the UV luminosity function.

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

We present catalogs of calibrated photometry and spectroscopic redshifts in the Extended Groth Strip, intended for studies of photometric redshifts (photo-z's). The data includes ugriz photometry from CFHTLS and Y-band photometry from the Subaru Suprime camera, as well as spectroscopic redshifts from the DEEP2, DEEP3 and 3D-HST surveys. These catalogs incorporate corrections to produce effectively matched-aperture photometry across all bands, based upon object size information available in the catalog and Moffat profile point spread function fits. We test this catalog with a simple machine learning-based photometric redshift algorithm based upon Random Forest regression, and find that the corrected aperture photometry leads to significant improvement in photo-z accuracy compared to the original SExtractor catalogs from CFHTLS and Subaru. The deep ugrizY photometry and spectroscopic redshifts are well-suited for empirical tests of photometric redshift algorithms for LSST. The resulting catalogs are publicly available. We include a basic summary of the strategy of the DEEP3 Galaxy Redshift Survey to accompany the recent public release of DEEP3 data.

We present a WFC3 F160W ($H$-band) selected catalog in the CANDELS/GOODS-N field containing photometry from the ultraviolet (UV) to the far-infrared (IR), photometric redshifts and stellar parameters derived from the analysis of the multi-wavelength data. The catalog contains 35,445 sources over the 171 arcmin$^{2}$ of the CANDELS F160W mosaic. The 5$\sigma$ detection limits (within an aperture of radius 0\farcs17) of the mosaic range between $H=27.8$, 28.2 and 28.7 in the wide, intermediate and deep regions, that span approximately 50\%, 15\% and 35\% of the total area. The multi-wavelength photometry includes broad-band data from UV (U band from KPNO and LBC), optical (HST/ACS F435W, F606W, F775W, F814W, and F850LP), near-to-mid IR (HST/WFC3 F105W, F125W, F140W and F160W, Subaru/MOIRCS Ks, CFHT/Megacam K, and \spitzer/IRAC 3.6, 4.5, 5.8, 8.0 $\mu$m) and far IR (\spitzer/MIPS 24$\mu$m, HERSCHEL/PACS 100 and 160$\mu$m, SPIRE 250, 350 and 500$\mu$m) observations. In addition, the catalog also includes optical medium-band data (R$\sim50$) in 25 consecutive bands, $\lambda=500$ to 950~nm, from the SHARDS survey and WFC3 IR spectroscopic observations with the G102 and G141 grisms (R$\sim210$ and 130). The use of higher spectral resolution data to estimate photometric redshifts provides very high, and nearly uniform, precision from $z=0-2.5$. The comparison to 1,485 good quality spectroscopic redshifts up to $z\sim3$ yields $\Delta z$/(1+$z_{\rm spec}$)$=$0.0032 and an outlier fraction of $\eta=$4.3\%. In addition to the multi-band photometry, we release added-value catalogs with emission line fluxes, stellar masses, dust attenuations, UV- and IR- based star formation rates and rest-frame colors.

This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `Outer Rim' run. It covers 440 deg^2 of sky area to a redshift of z=3 and is complete to a magnitude depth of 28 in the r-band. Each galaxy is characterized by a multitude of properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-driven empirical approaches with semi-analytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data.

MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at 2<z<5. A 6.5-m Magellan telescope will be coupled with DESI spectrographs to achieve multiplexing of 20,000. MegaMapper would be located at Las Campanas Observatory to fully access LSST imaging for target selection.

We present the status of the Dark Energy Spectroscopic Instrument (DESI) and its plans and opportunities for the coming decade. DESI construction and its initial five years of operations are an approved experiment of the US Department of Energy and is summarized here as context for the Astro2020 panel. Beyond 2025, DESI will require new funding to continue operations. We expect that DESI will remain one of the world's best facilities for wide-field spectroscopy throughout the decade. More about the DESI instrument and survey can be found at https://www.desi.lbl.gov.

The Large Synoptic Survey Telescope (LSST) can advance scientific frontiers beyond its groundbreaking 10-year survey. Here we explore opportunities for extended operations with proposal-based observing strategies, new filters, or transformed instrumentation. We recommend the development of a mid-decade community- and science-driven process to define next-generation LSST capabilities.

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.

The Maunakea Spectroscopic Explorer is a next-generation massively multiplexed spectroscopic facility currently under development in Hawaii. It is completely dedicated to large-scale spectroscopic surveys and will enable transformative science. In this white paper we summarize the science case and describe the current state of the project.

Hubble Space Telescope observations show that low-mass ($M_*=10^9-10^{10}M_{\odot}$) galaxies at high redshift ($z=1.0-2.5$) tend to be elongated (prolate) rather than disky (oblate) or spheroidal. This is explained in zoom-in cosmological hydrodynamical simulations by the fact that these galaxies are forming in cosmic web filaments where accretion happens preferentially along the direction of elongation. We ask whether the elongated morphology of these galaxies allows them to be used as effective tracers of cosmic web filaments at high redshift via their intrinsic alignments. Using mock lightcones and spectroscopically-confirmed galaxy pairs from the CANDELS survey, we test two types of alignments: (1) between the galaxy major axis and the direction to nearby galaxies of any mass, and (2) between the major axes of nearby pairs of low-mass, likely prolate, galaxies. The mock lightcones predict strong signals in 3D real space, 3D redshift space, and 2D projected redshift space for both types of alignments (assuming prolate galaxy orientations are the same as those of their host prolate halos), but we do not detect significant alignment signals in CANDELS observations. However, we show that spectroscopic redshifts have been obtained for only a small fraction of highly elongated galaxies, and accounting for spectroscopic incompleteness and redshift errors significantly degrades the 2D mock signal. This may partly explain the alignment discrepancy and highlights one of several avenues for future work.

On small scales there have been a number of claims of discrepancies between the standard Cold Dark Matter (CDM) model and observations. The 'missing satellites problem' infamously describes the overabundance of subhalos from CDM simulations compared to the number of satellites observed in the Milky Way. A variety of solutions to this discrepancy have been proposed; however, the impact of the specific properties of the Milky Way halo relative to the typical halo of its mass have yet to be explored. Motivated by recent studies that identified ways in which the Milky Way is atypical (e.g., Licquia et al. 2015), we investigate how the properties of dark matter halos with mass comparable to our Galaxy's --- including concentration, spin, shape, and scale factor of the last major merger --- correlate with the subhalo abundance. Using zoom-in simulations of Milky Way-like halos, we build two models of subhalo abundance as functions of host halo properties and conclude that the Milky Way should be expected to have 22%-44% fewer subhalos with low maximum rotation velocities ($V_{\rm max}^{\rm sat} \sim 10$kms$^{-1}$) at the 95% confidence level and up to 72% fewer than average subhalos with high rotation velocities ($V_{\rm max}^{\rm sat} \gtrsim 30$kms$^{-1}$, comparable to the Magellanic Clouds) than would be expected for a typical halo of the Milky Way's mass. Concentration is the most informative single parameter for predicting subhalo abundance. Our results imply that models tuned to explain the missing satellites problem assuming typical subhalo abundances for our Galaxy will be over-correcting.

This paper presents stellar mass functions and i-band luminosity functions for Sloan Digital Sky Survey (SDSS) galaxies at $i < 21$ using clustering redshifts, from which we also compute targeting completeness measurements for the Baryon Oscillation Spectroscopic Survey (BOSS). Clustering redshifts is a method of obtaining the redshift distribution of a sample of galaxies with only photometric information by measuring the angular crosscorrelation with a spectroscopic sample in different redshift bins. We construct a spectroscopic sample containing data from the BOSS + eBOSS surveys, allowing us to recover redshift distributions from photometric data out to $z\simeq 2.5$. We produce k-corrected i-band luminosity functions and stellar mass functions by applying clustering redshifts to SDSS DR8 galaxies in small bins of colour and magnitude. There is little evolution in the mass function between $0.2 < z < 0.8$, implying the most massive galaxies form most of their mass before $z = 0.8$. These mass functions are used to produce stellar mass completeness estimates for the Baryon Oscillation Spectroscopic Survey (BOSS), giving a stellar mass completeness of $80\%$ above $M_{\star} > 10^{11.4}$ between $0.2 < z < 0.7$, with completeness falling significantly at redshifts higher than 0.7, and at lower masses. Large photometric datasets will be available in the near future (DECaLS, DES, Euclid), so this, and similar techniques will become increasingly useful in order to fully utilise this data.

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

We propose an extension of the LSST survey to cover the northern sky to DEC < +30 (accessible at airmass <1.8). This survey will increase the LSST sky coverage by ~9,600 square degrees from 18,900 to 28,500 square degrees (a 50% increase) but use only 0.6-2.5% of the time depending on the synergies with other surveys. This increased area addresses a wide range of science cases that enhance all of the primary LSST science goals by significant amounts. The science enabled includes: increasing the area of the sky accessible for follow-up of multi-messenger transients including gravitational waves, mapping the milky way halo and halo dwarfs including discovery of RR Lyrae stars in the outer galactic halo, discovery of z>7 quasars in combination Euclid, enabling a second generation DESI and other spectroscopic surveys, and enhancing all areas of science by improving synergies with Euclid, WFIRST, and unique northern survey facilities. This white paper is the result of the Tri-Agency Working Group (TAG) appointed to develop synergies between missions and presents a unified plan for northern coverage. The range of time estimates reflects synergies with other surveys. If the modified DESC WFD survey, the ecliptic plane mini survey, and the north galactic spur mini survey are executed this plan would only need 0.6% of the LSST time, however if none of these are included the overall request is 2.5% of the 10 year survey life. In other words, the majority of these observations are already suggested as part of these other surveys and the intent of this white paper is to propose a unified baseline plan to carry out a broad range of objectives to facilitate a combination of multiple science objectives. A companion white paper gives Euclid specific science goals, and we support the white papers for southern extensions of the LSST survey.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

LSST will open new vistas for cosmology in the next decade, but it cannot reach its full potential without data from other telescopes. Cosmological constraints can be greatly enhanced using wide-field ($>20$ deg$^2$ total survey area), highly-multiplexed optical and near-infrared multi-object spectroscopy (MOS) on 4-15m telescopes. This could come in the form of suitably-designed large surveys and/or community access to add new targets to existing projects. First, photometric redshifts can be calibrated with high precision using cross-correlations of photometric samples against spectroscopic samples at $0 < z < 3$ that span thousands of sq. deg. Cross-correlations of faint LSST objects and lensing maps with these spectroscopic samples can also improve weak lensing cosmology by constraining intrinsic alignment systematics, and will also provide new tests of modified gravity theories. Large samples of LSST strong lens systems and supernovae can be studied most efficiently by piggybacking on spectroscopic surveys covering as much of the LSST extragalactic footprint as possible (up to $\sim20,000$ square degrees). Finally, redshifts can be measured efficiently for a high fraction of the supernovae in the LSST Deep Drilling Fields (DDFs) by targeting their hosts with wide-field spectrographs. Targeting distant galaxies, supernovae, and strong lens systems over wide areas in extended surveys with (e.g.) DESI or MSE in the northern portion of the LSST footprint or 4MOST in the south could realize many of these gains; DESI, 4MOST, Subaru/PFS, or MSE would all be well-suited for DDF surveys. The most efficient solution would be a new wide-field, highly-multiplexed spectroscopic instrument in the southern hemisphere with $>6$m aperture. In two companion white papers we present gains from deep, small-area MOS and from single-target imaging and spectroscopy.

This document summarizes the science cases related to cosmology studies with the MaunaKea Spectroscopic Explorer (MSE), a highly-multiplexed (4332 fibers), wide FOV (1.5 sq deg), large aperture (11.25 m in diameter), optical/NIR (360nm to 1300nm) facility. The MSE High-z Cosmology Survey is designed to probe a large volume of the Universe with a galaxy density sufficient to measure the extremely-large-scale density fluctuations required to explore primordial non-Gaussianity and therefore inflation. We expect a measurement of the local parameter $f_{NL}$ to a precision $\sigma(f_{NL}) = 1.8$. Combining the MSE High-z Cosmology Survey data with data from a next generation CMB stage 4 experiment and existing DESI data will provide the first $5\sigma$ confirmation of the neutrino mass hierarchy from astronomical observations. In addition, the Baryonic Acoustic Oscillations (BAO) observed within the sample will provide measurements of the distance-redshift relationship in six different redshift bins between $z=1.6$ and 4.0, each with an accuracy of $\sim0.6\%$. The simultaneous measurements of Redshift Space Distortions (RSD) constrain the amplitude of the fluctuations, at a level ranging from $1.9\%$ to $3.6\%$. The proposed survey covers 10,000 ${\rm deg}^2$, measuring redshifts for three classes of target objects: Emission Line Galaxies (ELGs) with $1.6<z<2.4$, Lyman Break Galaxies (LBGs) with $2.4<z<4.0$, and quasars $2.1<z<3.5$. The ELGs and LBGs will be used as direct tracers of the underlying density field, while the Lyman-$\alpha$ forests in the quasar spectra will be utilized to probe structure. Exposures of duration 1,800sec will guarantee a redshift determination efficiency of $90\%$ for ELGS and at least $50\%$ for LBGs. The survey will represent 100 nights per year for a 5-year MSE program. Finally, three ideas for additional projects of cosmological interest are proposed.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

Community access to deep (i ~ 25), highly-multiplexed optical and near-infrared multi-object spectroscopy (MOS) on 8-40m telescopes would greatly improve measurements of cosmological parameters from LSST. The largest gain would come from improvements to LSST photometric redshifts, which are employed directly or indirectly for every major LSST cosmological probe; deep spectroscopic datasets will enable reduced uncertainties in the redshifts of individual objects via optimized training. Such spectroscopy will also determine the relationship of galaxy SEDs to their environments, key observables for studies of galaxy evolution. The resulting data will also constrain the impact of blending on photo-z's. Focused spectroscopic campaigns can also improve weak lensing cosmology by constraining the intrinsic alignments between the orientations of galaxies. Galaxy cluster studies can be enhanced by measuring motions of galaxies in and around clusters and by testing photo-z performance in regions of high density. Photometric redshift and intrinsic alignment studies are best-suited to instruments on large-aperture telescopes with wider fields of view (e.g., Subaru/PFS, MSE, or GMT/MANIFEST) but cluster investigations can be pursued with smaller-field instruments (e.g., Gemini/GMOS, Keck/DEIMOS, or TMT/WFOS), so deep MOS work can be distributed amongst a variety of telescopes. However, community access to large amounts of nights for surveys will still be needed to accomplish this work. In two companion white papers we present gains from shallower, wide-area MOS and from single-target imaging and spectroscopy.

Single-object imaging and spectroscopy on telescopes with apertures ranging from ~4 m to 40 m have the potential to greatly enhance the cosmological constraints that can be obtained from LSST. Two major cosmological probes will benefit greatly from LSST follow-up: accurate spectrophotometry for nearby and distant Type Ia supernovae will expand the cosmological distance lever arm by unlocking the constraining power of high-z supernovae; and cosmology with time delays of strongly-lensed supernovae and quasars will require additional high-cadence imaging to supplement LSST, adaptive optics imaging or spectroscopy for accurate lens and source positions, and IFU or slit spectroscopy to measure detailed properties of lens systems. We highlight the scientific impact of these two science drivers, and discuss how additional resources will benefit them. For both science cases, LSST will deliver a large sample of objects over both the wide and deep fields in the LSST survey, but additional data to characterize both individual systems and overall systematics will be key to ensuring robust cosmological inference to high redshifts. Community access to large amounts of natural-seeing imaging on ~2-4 m telescopes, adaptive optics imaging and spectroscopy on 8-40 m telescopes, and high-throughput single-target spectroscopy on 4-40 m telescopes will be necessary for LSST time domain cosmology to reach its full potential. In two companion white papers we present the additional gains for LSST cosmology that will come from deep and from wide-field multi-object spectroscopy.

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

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

The nature of dark energy remains a profound mystery 20 years after the discovery of cosmic acceleration. A very high number density galaxy redshift survey over a wide area (HD GRS Wide) spanning the redshift range of 0.5<z<4 using the same tracer, carried out using massively parallel wide field multi-object slit spectroscopy from space, will provide definitive dark energy measurements with minimal observational systematics by design. The HD GRS Wide will illuminate the nature of dark energy, and lead to revolutionary advances in particle physics and cosmology. It will also trace the cosmic web of dark matter and provide key insight into large-scale structure in the Universe. The required observational facility can be established as part of the probe portfolio by NASA within the next decade.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

We use data from the All Wavelength Extended Groth Strip International Survey to construct stacked X-ray maps of optically bright active galaxies (AGN) and an associated control sample of galaxies at high redshift (z less than 0.6). From our analysis of the surface brightness profiles obtained from these X-ray maps, we find evidence of feedback from the active nuclei. We find that excluding galaxies and AGN, residing in group environments, from our samples enhances the significance of our detection. Our results support the tentative findings of Chatterjee et al. who use X-ray selected AGN for their analysis. We discuss the implications of these results in the context of quantifying AGN feedback and show that the current method can be used to extract X-ray source population in high redshift galaxies.

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

Cosmology is one of the four science pillars of LSST, which promises to be transformative for our understanding of dark energy and dark matter. The LSST Dark Energy Science Collaboration (DESC) has been tasked with deriving constraints on cosmological parameters from LSST data. Each of the cosmological probes for LSST is heavily impacted by the choice of observing strategy. This white paper is written by the LSST DESC Observing Strategy Task Force (OSTF), which represents the entire collaboration, and aims to make recommendations on observing strategy that will benefit all cosmological analyses with LSST. It is accompanied by the DESC DDF (Deep Drilling Fields) white paper (Scolnic et al.). We use a variety of metrics to understand the effects of the observing strategy on measurements of weak lensing, large-scale structure, clusters, photometric redshifts, supernovae, strong lensing and kilonovae. In order to reduce systematic uncertainties, we conclude that the current baseline observing strategy needs to be significantly modified to result in the best possible cosmological constraints. We provide some key recommendations: moving the WFD (Wide-Fast-Deep) footprint to avoid regions of high extinction, taking visit pairs in different filters, changing the 2x15s snaps to a single exposure to improve efficiency, focusing on strategies that reduce long gaps (>15 days) between observations, and prioritizing spatial uniformity at several intervals during the 10-year survey.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

The Wide-Field Infrared Survey Telescope (WFIRST) is expected to launch in the mid-2020s. With its wide-field near-infrared (NIR) camera, it will survey the sky to unprecedented detail. As part of normal operations and as the result of multiple expected dedicated surveys, WFIRST will produce several relatively wide-field (tens of square degrees) deep (limiting magnitude of 28 or fainter) fields. In particular, a planned supernova survey is expected to image 3 deep fields in the LSST footprint roughly every 5 days over 2 years. Stacking all data, this survey will produce, over all WFIRST supernova fields in the LSST footprint, ~12-25 deg^2 and ~5-15 deg^2 regions to depths of ~28 mag and ~29 mag, respectively. We suggest LSST undertake mini-surveys that will match the WFIRST cadence and simultaneously observe the supernova survey fields during the 2-year WFIRST supernova survey, achieving a stacked depth similar to that of the WFIRST data. We also suggest additional observations of these same regions throughout the LSST survey to get deep images earlier, have long-term monitoring in the fields, and produce deeper images overall. These fields will provide a legacy for cosmology, extragalactic, and transient/variable science.

Double-peaked narrow emission lines in active galactic nucleus (AGN) spectra can be produced by AGN outflows, rotation, or dual AGNs, which are AGN pairs in ongoing galaxy mergers. Consequently, double-peaked narrow AGN emission lines are useful tracers of the coevolution of galaxies and their supermassive black holes, as driven by AGN feedback and AGN fueling. We investigate this concept further with follow-up optical longslit observations of a sample of 95 Sloan Digital Sky Survey (SDSS) galaxies that have double-peaked narrow AGN emission lines. Based on a kinematic analysis of the longslit spectra, we confirm previous work that finds that the majority of double-peaked narrow AGN emission lines are associated with outflows. We also find that eight of the galaxies have companion galaxies with line-of-sight velocity separations < 500 km/s and physical separations <30 kpc. Since we find evidence of AGNs in both galaxies, all eight of these systems are compelling dual AGN candidates. Galaxies with double-peaked narrow AGN emission lines occur in such galaxy mergers at least twice as often as typical active galaxies. Finally, we conclude that at least 3% of SDSS galaxies with double-peaked narrow AGN emission lines are found in galaxy mergers where both galaxies are resolved in SDSS imaging.

To achieve a fuller understanding of galaxy evolution, SED fitting can be used to recover quantities beyond stellar masses (M$_*$) and star formation rates (SFRs). We use Star Formation Histories (SFHs) reconstructed via the Dense Basis method of Iyer \& Gawiser (2017) for a sample of $17,873$ galaxies at $0.5<z<6$ in the CANDELS GOODS-S field to study the nature and evolution of the SFR-M$_*$ correlation. The reconstructed SFHs represent trajectories in SFR-M$_*$ space, enabling us to study galaxies at epochs earlier than observed by propagating them backwards in time along these trajectories. We study the SFR-M$_*$ correlation at $z=1,2,3,4,5,6$ using both direct fits to galaxies observed at those epochs and SFR-M$_*$ trajectories of galaxies observed at lower redshifts. The SFR-M$_*$ correlations obtained using the two approaches are found to be consistent with each other through a KS test. Validation tests using SFHs from semi-analytic models and cosmological hydrodynamical simulations confirm the sensitivity of the method to changes in the slope, normalization and shape of the SFR-M$_*$ correlation. This technique allows us to further probe the low-mass regime of the correlation at high-z by $\sim 1$ dex and over an effective volume of $\sim 10\times$ larger than possible with just direct fits. We find that the SFR-M$_*$ correlation is consistent with being linear down to M$_*\sim 10^7 M_\odot$ at $z>4$. The evolution of the correlation is well described by $\log SFR= (0.80\pm 0.029 - 0.017\pm 0.010\times t_{univ})\log M_*$ $- (6.487\pm 0.282-0.039\pm 0.008\times t_{univ})$, where $t_{univ}$ is the age of the universe in Gyr.

The Large Synoptic Survey Telescope (LSST) Dark Energy Science Collaboration (DESC) will use five cosmological probes: galaxy clusters, large scale structure, supernovae, strong lensing, and weak lensing. This Science Requirements Document (SRD) quantifies the expected dark energy constraining power of these probes individually and together, with conservative assumptions about analysis methodology and follow-up observational resources based on our current understanding and the expected evolution within the field in the coming years. We then define requirements on analysis pipelines that will enable us to achieve our goal of carrying out a dark energy analysis consistent with the Dark Energy Task Force definition of a Stage IV dark energy experiment. This is achieved through a forecasting process that incorporates the flowdown to detailed requirements on multiple sources of systematic uncertainty. Future versions of this document will include evolution in our software capabilities and analysis plans along with updates to the LSST survey strategy.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

This is the first in a series of papers examining the demographics of star-forming galaxies at $0.2<z<2.5$ in CANDELS. We study 9,100 galaxies from GOODS-S and UDS having published values of redshifts, masses, star-formation rates (SFRs), and dust attenuation ($A_V$) derived from UV-optical SED fitting. In agreement with previous works, we find that the $UVJ$ colors of a galaxy are closely correlated with its specific star-formation rate (SSFR) and $A_V$. We define rotated $UVJ$ coordinate axes, termed $S_\mathrm{SED}$ and $C_\mathrm{SED}$, that are parallel and perpendicular to the star-forming sequence and derive a quantitative calibration that predicts SSFR from $C_\mathrm{SED}$ with an accuracy of ~0.2 dex. SFRs from UV-optical fitting and from UV+IR values based on Spitzer/MIPS 24 $\mu\mathrm{m}$ agree well overall, but systematic differences of order 0.2 dex exist at high and low redshifts. A novel plotting scheme conveys the evolution of multiple galaxy properties simultaneously, and dust growth, as well as star-formation decline and quenching, exhibit "mass-accelerated evolution" ("downsizing"). A population of transition galaxies below the star-forming main sequence is identified. These objects are located between star-forming and quiescent galaxies in $UVJ$ space and have lower $A_V$ and smaller radii than galaxies on the main sequence. Their properties are consistent with their being in transit between the two regions. The relative numbers of quenched, transition, and star-forming galaxies are given as a function of mass and redshift.

MaNGA provides the opportunity to make precise spatially resolved measurements of the IMF slope in galaxies owing to its unique combination of spatial resolution, wavelength coverage and sample size. We derive radial gradients in age, element abundances and IMF slope analysing optical and near-infrared absorption features from stacked spectra out to the half-light radius of 366 early-type galaxies with masses $9.9 - 10.8\;\log M/M_{\odot}$. We find flat gradients in age and [$\alpha$/Fe] ratio, as well as negative gradients in metallicity, consistent with the literature. We further derive significant negative gradients in the [Na/Fe] ratio with galaxy centres being well enhanced in Na abundance by up to 0.5 dex. Finally, we find a gradient in IMF slope with a bottom-heavy IMF in the centre (typical mass excess factor of 1.5) and a Milky Way-type IMF at the half-light radius. This pattern is mass-dependent with the lowest mass galaxies in our sample featuring only a shallow gradient around a Milky Way IMF. Our results imply the local IMF-$\sigma$ relation within galaxies to be even steeper than the global relation and hint towards the local metallicity being the dominating factor behind the IMF variations. We also employ different stellar population models in our analysis and show that a radial IMF gradient is found independently of the stellar population model used. A similar analysis of the Wing-Ford band provides inconsistent results and further evidence of the difficulty in measuring and modelling this particular feature.

The use of high-quality simulated sky catalogs is essential for the success of cosmological surveys. The catalogs have diverse applications, such as investigating signatures of fundamental physics in cosmological observables, understanding the effect of systematic uncertainties on measured signals and testing mitigation strategies for reducing these uncertainties, aiding analysis pipeline development and testing, and survey strategy optimization. The list of applications is growing with improvements in the quality of the catalogs and the details that they can provide. Given the importance of simulated catalogs, it is critical to provide rigorous validation protocols that enable both catalog providers and users to assess the quality of the catalogs in a straightforward and comprehensive way. For this purpose, we have developed the DESCQA framework for the Large Synoptic Survey Telescope Dark Energy Science Collaboration as well as for the broader community. The goal of DESCQA is to enable the inspection, validation, and comparison of an inhomogeneous set of synthetic catalogs via the provision of a common interface within an automated framework. In this paper, we present the design concept and first implementation of DESCQA. In order to establish and demonstrate its full functionality we use a set of interim catalogs and validation tests. We highlight several important aspects, both technical and scientific, that require thoughtful consideration when designing a validation framework, including validation metrics and how these metrics impose requirements on the synthetic sky catalogs.

The extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 sample includes 80,118 Luminous Red Galaxies. By combining these galaxies with the high-redshift tail of the BOSS galaxy sample, we form a sample of LRGs at an effective redshift $z=0.72$, covering an effective volume of 0.9~Gpc$^3$. We introduce new techniques to account for spurious fluctuations caused by targeting and by redshift failures which were validated on a set of mock catalogs. This analysis is sufficient to provide a $2.6$\% measurement of spherically averaged BAO, $D_V(z=0.72) = 2353^{+63}_{-61} (r_d/r_{d,\rm{fid}}) h^{-1}$Mpc, at 2.8$\sigma$ of significance. Together with the recent quasar-based BAO measurement at $z=1.5$, and forthcoming Emission Line Galaxy-based measurements, this measurement demonstrates that eBOSS is fulfilling its remit of extending the range of redshifts covered by such measurements, laying the ground work for forthcoming surveys such as the Dark Energy Spectroscopic Survey and Euclid.

The rate of major galaxy-galaxy merging is theoretically predicted to steadily increase with redshift during the peak epoch of massive galaxy development ($1{\leq}z{\leq}3$). We use close-pair statistics to objectively study the incidence of massive galaxies (stellar $M_{1}{\geq}2{\times}10^{10}M_{\odot}$) hosting major companions ($1{\leq}M_{1}/M_{2}{\leq}4$; i.e., $<$4:1) at six epochs spanning $0{<}z{<}3$. We select companions from a nearly complete, mass-limited ($\geq5{\times}10^{9}M_{\odot}$) sample of 23,696 galaxies in the five CANDELS fields and the SDSS. Using $5-50$ kpc projected separation and close redshift proximity criteria, we find that the major companion fraction $f_{\mathrm{mc}}(z)$ based on stellar mass-ratio (MR) selection increases from 6% ($z{\sim}0$) to 16% ($z{\sim}0.8$), then turns over at $z{\sim}1$ and decreases to 7% ($z{\sim}3$). Instead, if we use a major F160W flux ratio (FR) selection, we find that $f_{\mathrm{mc}}(z)$ increases steadily until $z=3$ owing to increasing contamination from minor (MR$>$4:1) companions at $z>1$. We show that these evolutionary trends are statistically robust to changes in companion proximity. We find disagreements between published results are resolved when selection criteria are closely matched. If we compute merger rates using constant fraction-to-rate conversion factors ($C_{\mathrm{merg,pair}}{=}0.6$ and $T_{\mathrm{obs,pair}}{=}0.65\mathrm{Gyr}$), we find that MR rates disagree with theoretical predictions at $z{>}1.5$. Instead, if we use an evolving $T_{\mathrm{obs,pair}}(z){\propto}(1+z)^{-2}$ from Snyder et al., our MR-based rates agree with theory at $0{<}z{<}3$. Our analysis underscores the need for detailed calibration of $C_{\mathrm{merg,pair}}$ and $T_{\mathrm{obs,pair}}$ as a function of redshift, mass and companion selection criteria to better constrain the empirical major merger history.

Euclid and the Large Synoptic Survey Telescope (LSST) are poised to dramatically change the astronomy landscape early in the next decade. The combination of high cadence, deep, wide-field optical photometry from LSST with high resolution, wide-field optical photometry and near-infrared photometry and spectroscopy from Euclid will be powerful for addressing a wide range of astrophysical questions. We explore Euclid/LSST synergy, ignoring the political issues associated with data access to focus on the scientific, technical, and financial benefits of coordination. We focus primarily on dark energy cosmology, but also discuss galaxy evolution, transient objects, solar system science, and galaxy cluster studies. We concentrate on synergies that require coordination in cadence or survey overlap, or would benefit from pixel-level co-processing that is beyond the scope of what is currently planned, rather than scientific programs that could be accomplished only at the catalog level without coordination in data processing or survey strategies. We provide two quantitative examples of scientific synergies: the decrease in photo-z errors (benefitting many science cases) when high resolution Euclid data are used for LSST photo-z determination, and the resulting increase in weak lensing signal-to-noise ratio from smaller photo-z errors. We briefly discuss other areas of coordination, including high performance computing resources and calibration data. Finally, we address concerns about the loss of independence and potential cross-checks between the two missions and potential consequences of not collaborating.

We present measurements of the Baryon Acoustic Oscillation (BAO) scale in redshift-space using the clustering of quasars. We consider a sample of 147,000 quasars from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) distributed over 2044 square degrees with redshifts $0.8 < z < 2.2$ and measure their spherically-averaged clustering in both configuration and Fourier space. Our observational dataset and the 1400 simulated realizations of the dataset allow us to detect a preference for BAO that is greater than 2.8$\sigma$. We determine the spherically averaged BAO distance to $z = 1.52$ to 3.8 per cent precision: $D_V(z=1.52)=3843\pm147 \left(r_{\rm d}/r_{\rm d, fid}\right)\ $Mpc. This is the first time the location of the BAO feature has been measured between redshifts 1 and 2. Our result is fully consistent with the prediction obtained by extrapolating the Planck flat $\Lambda$CDM best-fit cosmology. All of our results are consistent with basic large-scale structure (LSS) theory, confirming quasars to be a reliable tracer of LSS, and provide a starting point for numerous cosmological tests to be performed with eBOSS quasar samples. We combine our result with previous, independent, BAO distance measurements to construct an updated BAO distance-ladder. Using these BAO data alone and marginalizing over the length of the standard ruler, we find $\Omega_{\Lambda} > 0$ at 6.6$\sigma$ significance when testing a $\Lambda$CDM model with free curvature.

We present the first scientific results from the luminous red galaxy sample (LRG) of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We measure the small and intermediate scale clustering from a sample of more than 61,000 galaxies in the redshift range $0.6 < z < 0.9$. We interpret these measurements in the framework of the Halo Occupation Distribution. The bias of eBOSS LRGs is $2.30 \pm 0.03$, with a satellite fraction of $13\pm3$\% and a mean halo mass of $2.5\times10^{13}h^{-1}M_{\odot}$. These results are consistent with expectations, demonstrating that eBOSS galaxies will be reliable tracers of large scale structure at $z\sim 0.7$. The eBOSS galaxy bias implies a scatter of luminosity at fixed halo mass, $\sigma_{\log L}$, of 0.19 dex. Using the clustering of massive galaxies from BOSS-CMASS, BOSS-LOWZ, and SDSS, we find that $\sigma_{\log L}=0.19$ is consistent with observations over the full redshift range that these samples cover. The addition of eBOSS to previous surveys allows investigation of the evolution of massive galaxies over the past $\sim 7$ Gyr.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

Two narrow and dense rings (called C1R and C2R) were discovered around the Centaur object (10199) Chariklo during a stellar occultation observed on 2013 June 3. Following this discovery, we planned observations of several occultations by Chariklo's system in order to better characterize the physical properties of the ring and main body. Here, we use 12 successful occulations by Chariklo observed between 2014 and 2016. They provide ring profiles (physical width, opacity, edge structure) and constraints on the radii and pole position. Our new observations are currently consistent with the circular ring solution and pole position, to within the $\pm 3.3$ km formal uncertainty for the ring radii derived by Braga-Ribas et al. The six resolved C1R profiles reveal significant width variations from $\sim 5$ to 7.5 km. The width of the fainter ring C2R is less constrained, and may vary between 0.1 and 1 km. The inner and outer edges of C1R are consistent with infinitely sharp boundaries, with typical upper limits of one kilometer for the transition zone between the ring and empty space. No constraint on the sharpness of C2R's edges is available. A 1$\sigma$ upper limit of $\sim 20$ m is derived for the equivalent width of narrow (physical width <4 km) rings up to distances of 12,000 km, counted in the ring plane.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

The Large Synoptic Survey Telescope (LSST) will enable revolutionary studies of galaxies, dark matter, and black holes over cosmic time. The LSST Galaxies Science Collaboration has identified a host of preparatory research tasks required to leverage fully the LSST dataset for extragalactic science beyond the study of dark energy. This Galaxies Science Roadmap provides a brief introduction to critical extragalactic science to be conducted ahead of LSST operations, and a detailed list of preparatory science tasks including the motivation, activities, and deliverables associated with each. The Galaxies Science Roadmap will serve as a guiding document for researchers interested in conducting extragalactic science in anticipation of the forthcoming LSST era.

We explore observational and theoretical constraints on how galaxies might transition between the "star-forming main sequence" (SFMS) and varying "degrees of quiescence" out to $z=3$. Our analysis is focused on galaxies with stellar mass $M_*>10^{10}M_{\odot}$, and is enabled by GAMA and CANDELS observations, a semi-analytic model (SAM) of galaxy formation, and a cosmological hydrodynamical "zoom in" simulation with momentum-driven AGN feedback. In both the observations and the SAM, transition galaxies tend to have intermediate S\'ersic indices, half-light radii, and surface stellar mass densities compared to star-forming and quiescent galaxies out to $z=3$. We place an observational upper limit on the average population transition timescale as a function of redshift, finding that the average high-redshift galaxy is on a "fast track" for quenching whereas the average low-redshift galaxy is on a "slow track" for quenching. We qualitatively identify four physical origin scenarios for transition galaxies in the SAM: oscillations on the SFMS, slow quenching, fast quenching, and rejuvenation. Quenching timescales in both the SAM and the hydrodynamical simulation are not fast enough to reproduce the quiescent population that we observe at $z\sim3$. In the SAM, we do not find a clear-cut morphological dependence of quenching timescales, but we do predict that the mean stellar ages, cold gas fractions, SMBH masses, and halo masses of transition galaxies tend to be intermediate relative to those of star-forming and quiescent galaxies at $z<3$.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We present a 0.4-8$\mu$m multi-wavelength photometric catalog in the Extended Groth Strip (EGS) field. This catalog is built on the Hubble Space Telescope (HST) WFC3 and ACS data from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS), and it incorporates the existing HST data from the All-wavelength Extended Groth strip International Survey (AEGIS) and the 3D-HST program. The catalog is based on detections in the F160W band reaching a depth of F160W=26.62 AB (90% completeness, point-sources). It includes the photometry for 41457 objects over an area of $\approx 206$ arcmin$^2$ in the following bands: HST ACS F606W and F814W; HST WFC3 F125W, F140W and F160W; CFHT/Megacam $u^*$, $g'$, $r'$, $i'$ and $z'$; CFHT/WIRCAM $J$, $H$ and $K_\mathrm{S}$; Mayall/NEWFIRM $J1$, $J2$, $J3$, $H1$, $H2$, $K$; Spitzer IRAC $3.6\mu$m, $4.5\mu$m, $5.8\mu$m and $8.0\mu$m. We are also releasing value-added catalogs that provide robust photometric redshifts and stellar mass measurements. The catalogs are publicly available through the CANDELS repository.

We present a multi-wavelength photometric catalog in the COSMOS field as part of the observations by the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). The catalog is based on Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) and Advanced Camera for Surveys (ACS) observations of the COSMOS field (centered at RA: $10^h00^m28^s$, Dec:$+02^{\circ}12^{\prime}21^{\prime\prime}$). The final catalog has 38671 sources with photometric data in forty two bands from UV to the infrared ($\rm \sim 0.3-8\,\mu m$). This includes broad-band photometry from the HST, CFHT, Subaru, VISTA and Spitzer Space Telescope in the visible, near infrared and infrared bands along with intermediate and narrow-band photometry from Subaru and medium band data from Mayall NEWFIRM. Source detection was conducted in the WFC3 F160W band (at $\rm 1.6\,\mu m$) and photometry is generated using the Template FITting algorithm. We further present a catalog of the physical properties of sources as identified in the HST F160W band and measured from the multi-band photometry by fitting the observed spectral energy distributions of sources against templates.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

We present two wide-field catalogs of photometrically-selected emission line galaxies (ELGs) at z=0.8 covering about 2800 deg^2 over the south galactic cap. The catalogs were obtained using a Fisher discriminant technique described in a companion paper. The two catalogs differ by the imaging used to define the Fisher discriminant: the first catalog includes imaging from the Sloan Digital Sky Survey and the Wide-Field Infrared Survey Explorer, the second also includes information from the South Galactic Cap U-band Sky Survey (SCUSS). Containing respectively 560,045 and 615,601 objects, they represent the largest ELG catalogs available today and were designed for the ELG programme of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We study potential sources of systematic variation in the angular distribution of the selected ELGs due to fluctuations of the observational parameters. We model the influence of the observational parameters using a multivariate regression and implement a weighting scheme that allows effective removal of all of the systematic errors induced by the observational parameters. We show that fluctuations in the imaging zero-points of the photometric bands have minor impact on the angular distribution of objects in our catalogs. We compute the angular clustering of both catalogs and show that our weighting procedure effectively removes spurious clustering on large scales. We fit a model to the small scale angular clustering, showing that the selections have similar biases of 1.35/D_a(z) and 1.28/D_a(z). Both catalogs are publicly available.

We study the correlation of galaxy structural properties with their location relative to the SFR-M* correlation, also known as the star formation "main sequence" (SFMS), in the CANDELS and GAMA surveys and in a semi-analytic model (SAM) of galaxy formation. We first study the distribution of median Sersic index, effective radius, star formation rate (SFR) density and stellar mass density in the SFR-M* plane. We then define a redshift dependent main sequence and examine the medians of these quantities as a function of distance from this main sequence, both above (higher SFRs) and below (lower SFRs). Finally, we examine the distributions of distance from the main sequence in bins of these quantities. We find strong correlations between all of these galaxy structural properties and the distance from the SFMS, such that as we move from galaxies above the SFMS to those below it, we see a nearly monotonic trend towards higher median Sersic index, smaller radius, lower SFR density, and higher stellar density. In the semi-analytic model, bulge growth is driven by mergers and disk instabilities, and is accompanied by the growth of a supermassive black hole which can regulate or quench star formation via Active Galactic Nucleus (AGN) feedback. We find that our model qualitatively reproduces the trends described above, supporting a picture in which black holes and bulges co-evolve, and AGN feedback plays a critical role in moving galaxies off of the SFMS.

It is crucial to understand how the Milky Way, the galaxy we can study in the most intimate detail, fits in amongst other galaxies. Key examples include the Tully-Fisher relation (TFR) --- i.e., the tight correlation between luminosity ($L$) and rotational velocity ($V_\textrm{rot}$) --- and the 3-dimensional luminosity-velocity-radius ($LVR$) scaling relation. Several past studies have characterized the MW as a 1--1.5$\sigma$ outlier to the TFR. This study reexamines such comparisons using new estimates of MW properties that are robust to many of the systematic uncertainties that have been a problem in the past and are based on assumptions consistent with those used for other spiral galaxies. Comparing to scaling relations derived from modern extragalactic data, we find that our Galaxy's properties are in excellent agreement with TFRs defined using any SDSS-filter absolute magnitude, stellar mass, or baryonic mass as the $L$ proxy. We next utilize disk scale length ($R_\textrm{d}$) measurements to extend this investigation to the $LVR$ relation. Here we find that our Galaxy lies farther from the relation than $\sim$90\% of other spiral galaxies, yielding $\sim$9.5$\sigma$ evidence that it is unusually compact for its $L$ and $V_\textrm{rot}$ (based on MW errors alone), a result that holds for all of the $L$ proxies considered. The expected $R_\textrm{d}$ for the MW from the $LVR$ relation is $\sim$5 kpc, nearly twice as large as the observed value, with error estimates placing the two in tension at the $\sim$1.4$\sigma$ level. The compact scale length of the Galactic disk could be related to other ways in which the MW has been found to be anomalous.

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

We present quantified visual morphologies of approximately 48,000 galaxies observed in three Hubble Space Telescope legacy fields by the Cosmic And Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and classified by participants in the Galaxy Zoo project. 90% of galaxies have z < 3 and are observed in rest-frame optical wavelengths by CANDELS. Each galaxy received an average of 40 independent classifications, which we combine into detailed morphological information on galaxy features such as clumpiness, bar instabilities, spiral structure, and merger and tidal signatures. We apply a consensus-based classifier weighting method that preserves classifier independence while effectively down-weighting significantly outlying classifications. After analysing the effect of varying image depth on reported classifications, we also provide depth-corrected classifications which both preserve the information in the deepest observations and also enable the use of classifications at comparable depths across the full survey. Comparing the Galaxy Zoo classifications to previous classifications of the same galaxies shows very good agreement; for some applications the high number of independent classifications provided by Galaxy Zoo provides an advantage in selecting galaxies with a particular morphological profile, while in others the combination of Galaxy Zoo with other classifications is a more promising approach than using any one method alone. We combine the Galaxy Zoo classifications of "smooth" galaxies with parametric morphologies to select a sample of featureless disks at 1 < z < 3, which may represent a dynamically warmer progenitor population to the settled disk galaxies seen at later epochs.

The Large Synoptic Survey Telescope (LSST) will be a discovery machine for the astronomy and physics communities, revealing astrophysical phenomena from the Solar System to the outer reaches of the observable Universe. While many discoveries will be made using LSST data alone, taking full scientific advantage of LSST will require ground-based optical-infrared (OIR) supporting capabilities, e.g., observing time on telescopes, instrumentation, computing resources, and other infrastructure. This community-based study identifies, from a science-driven perspective, capabilities that are needed to maximize LSST science. Expanding on the initial steps taken in the 2015 OIR System Report, the study takes a detailed, quantitative look at the capabilities needed to accomplish six representative LSST-enabled science programs that connect closely with scientific priorities from the 2010 decadal surveys. The study prioritizes the resources needed to accomplish the science programs and highlights ways that existing, planned, and future resources could be positioned to accomplish the science goals.

The exponential scale length ($L_d$) of the Milky Way's (MW's) disk is a critical parameter for describing the global physical size of our Galaxy, important both for interpreting other Galactic measurements and helping us to understand how our Galaxy fits into extragalactic contexts. Unfortunately, current estimates span a wide range of values and often are statistically incompatible with one another. Here, we perform a Bayesian meta-analysis to determine an improved, aggregate estimate for $L_d$, utilizing a mixture-model approach to account for the possibility that any one measurement has not properly accounted for all statistical or systematic errors. Within this machinery we explore a variety of ways of modeling the nature of problematic measurements, and then employ a Bayesian model averaging technique to derive net posterior distributions that incorporate any model-selection uncertainty. Our meta-analysis combines 29 different (15 visible and 14 infrared) photometric measurements of $L_d$ available in the literature; these involve a broad assortment of observational datasets, MW models and assumptions, and methodologies, all tabulated herein. Analyzing the visible and infrared measurements separately yields estimates for $L_d$ of $2.71^{+0.22}_{-0.20}$ kpc and $2.51^{+0.15}_{-0.13}$ kpc, respectively, whereas considering them all combined yields $2.64\pm0.13$ kpc. The ratio between the visible and infrared scale lengths determined here is very similar to that measured in external spiral galaxies. We use these results to update the model of the Galactic disk from our previous work, constraining its stellar mass to be $4.8^{+1.5}_{-1.1}\times10^{10}$ $\textrm{M}_\odot$, and the MW's total stellar mass to be $5.7^{+1.5}_{-1.1}\times10^{10}$ $\textrm{M}_\odot$.

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) is an optical fiber-bundle integral-field unit (IFU) spectroscopic survey that is one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV). With a spectral coverage of 3622 - 10,354 Angstroms and an average footprint of ~ 500 arcsec^2 per IFU the scientific data products derived from MaNGA will permit exploration of the internal structure of a statistically large sample of 10,000 low redshift galaxies in unprecedented detail. Comprising 174 individually pluggable science and calibration IFUs with a near-constant data stream, MaNGA is expected to obtain ~ 100 million raw-frame spectra and ~ 10 million reduced galaxy spectra over the six-year lifetime of the survey. In this contribution, we describe the MaNGA Data Reduction Pipeline (DRP) algorithms and centralized metadata framework that produces sky-subtracted, spectrophotometrically calibrated spectra and rectified 3-D data cubes that combine individual dithered observations. For the 1390 galaxy data cubes released in Summer 2016 as part of SDSS-IV Data Release 13 (DR13), we demonstrate that the MaNGA data have nearly Poisson-limited sky subtraction shortward of ~ 8500 Angstroms and reach a typical 10-sigma limiting continuum surface brightness mu = 23.5 AB/arcsec^2 in a five arcsec diameter aperture in the g band. The wavelength calibration of the MaNGA data is accurate to 5 km/s rms, with a median spatial resolution of 2.54 arcsec FWHM (1.8 kpc at the median redshift of 0.037) and a median spectral resolution of sigma = 72 km/s.

Emission line galaxies are one of the main tracers of the large-scale structure to be targeted by the next-generation dark energy surveys. To provide a better understanding of the properties and statistics of these galaxies, we have collected spectroscopic data from the VVDS and DEEP2 deep surveys and estimated the galaxy luminosity functions (LFs) of three distinct emission lines, [OII] ($0.5 < z < 1.3$), H{\beta} ($0.3 < z < 0.8$) and [OIII] ($0.3 < z < 0.8$). Our measurements are based on 35,639 emission line galaxies and cover a volume of $\sim10^7$Mpc$^3$. We present the first measurement of the H{\beta} LF at these redshifts. We have also compiled LFs from the literature that were based on independent data or covered different redshift ranges, and we fit the entire set over the whole redshift range with analytic Schechter and Saunders models, assuming a natural redshift dependence of the parameters. We find that the characteristic luminosity ($L_*$) and density ($\phi_*$) of all LFs increase with redshift. Using the Schechter model over the redshift ranges considered, we find that, for [OII] emitters, the characteristic luminosity $L_*(z=0.5)=3.2\times10^{41}$ erg/s increases by a factor of $2.7 \pm 0.2$ from z=0.5 to 1.3; for H{\beta} emitters $L_*(z=0.3)=1.3\times10^{41}$ erg/s increases by a factor of $2.0 \pm 0.2$ from z=0.3 to 0.8; and for [OIII] emitters $L_*(z=0.3)=7.3\times10^{41}$ erg/s increases by a factor of $3.5 \pm 0.4$ from z=0.3 to 0.8.

The Sloan Digital Sky Survey IV extended Baryonic Oscillation Spectroscopic Survey (SDSS-IV/eBOSS) will observe 195,000 emission-line galaxies (ELGs) to measure the Baryonic Acoustic Oscillation standard ruler (BAO) at redshift 0.9. To test different ELG selection algorithms, 9,000 spectra were observed with the SDSS spectrograph as a pilot survey based on data from several imaging surveys. First, using visual inspection and redshift quality flags, we show that the automated spectroscopic redshifts assigned by the pipeline meet the quality requirements for a reliable BAO measurement. We also show the correlations between sky emission, signal-to-noise ratio in the emission lines, and redshift error. Then we provide a detailed description of each target selection algorithm we tested and compare them with the requirements of the eBOSS experiment. As a result, we provide reliable redshift distributions for the different target selection schemes we tested. Finally, we determine an target selection algorithms that is best suited to be applied on DECam photometry because they fulfill the eBOSS survey efficiency requirements.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

Dust attenuation affects nearly all observational aspects of galaxy evolution, yet very little is known about the form of the dust-attenuation law in the distant Universe. Here, we model the spectral energy distributions (SEDs) of galaxies at z = 1.5--3 from CANDELS with rest-frame UV to near-IR imaging under different assumptions about the dust law, and compare the amount of inferred attenuated light with the observed infrared (IR) luminosities. Some individual galaxies show strong Bayesian evidence in preference of one dust law over another, and this preference agrees with their observed location on the plane of infrared excess (IRX, $L_{\text{TIR}}/L_{\text{UV}}$) and UV slope ($\beta$). We generalize the shape of the dust law with an empirical model, $A_{\lambda,\delta}=E(B-V)\ k_\lambda\ (\lambda/\lambda_V)^\delta$ where $k_\lambda$ is the dust law of Calzetti et al. (2000), and show that there exists a correlation between the color excess ${E(B-V)}$ and tilt $\delta$ with ${\delta=(0.62\pm0.05)\log(E(B-V))}$+ ${(0.26~\pm~0.02)}$. Galaxies with high color excess have a shallower, starburst-like law, and those with low color excess have a steeper, SMC-like law. Surprisingly, the galaxies in our sample show no correlation between the shape of the dust law and stellar mass, star-formation rate, or $\beta$. The change in the dust law with color excess is consistent with a model where attenuation is caused by by scattering, a mixed star-dust geometry, and/or trends with stellar population age, metallicity, and dust grain size. This rest-frame UV-to-near-IR method shows potential to constrain the dust law at even higher ($z>3$) redshifts.

We present the stellar mass ($M_{*}$)--gas-phase metallicity relation (MZR) and its scatter at intermediate redshifts ($0.5\leq z\leq0.7$) for 1381 field galaxies collected from deep spectroscopic surveys. The star formation rate (SFR) and color at a given $M_{*}$ of this magnitude-limited ($R\lesssim24$ AB) sample are representative of normal star-forming galaxies. For masses below $10^9 M_\odot$, our sample of 237 galaxies is $\sim$10 times larger than those in previous studies beyond the local universe. This huge gain in sample size enables superior constraints on the MZR and its scatter in the low-mass regime. We find a power-law MZR at $10^{8} M_\odot < M_{*} < 10^{11} M_\odot$: ${12+log(O/H) = (5.83\pm0.19) + (0.30\pm0.02)log(M_{*}/M_\odot)}$. Our MZR shows good agreement with others measured at similar redshifts in the literature in the intermediate and massive regimes, but is shallower than the extrapolation of the MZRs of others to masses below $10^{9} M_\odot$. The SFR dependence of the MZR in our sample is weaker than that found for local galaxies (known as the Fundamental Metallicity Relation). Compared to a variety of theoretical models, the slope of our MZR for low-mass galaxies agrees well with predictions incorporating supernova energy-driven winds. Being robust against currently uncertain metallicity calibrations, the scatter of the MZR serves as a powerful diagnostic of the stochastic history of gas accretion, gas recycling, and star formation of low-mass galaxies. Our major result is that the scatter of our MZR increases as $M_{*}$ decreases. Our result implies that either the scatter of the baryonic accretion rate or the scatter of the $M_{*}$--$M_{halo}$ relation increases as $M_{*}$ decreases. Moreover, our measures of scatter at $z=0.7$ appears consistent with that found for local galaxies.

We present results from a multi-chord Pluto stellar occultation observed on 29 June 2015 from New Zealand and Australia. This occurred only two weeks before the NASA New Horizons flyby of the Pluto system and serves as a useful comparison between ground-based and space results. We find that Pluto's atmosphere is still expanding, with a significant pressure increase of 5$\pm$2\% since 2013 and a factor of almost three since 1988. This trend rules out, as of today, an atmospheric collapse associated with Pluto's recession from the Sun. A central flash, a rare occurrence, was observed from several sites in New Zealand. The flash shape and amplitude are compatible with a spherical and transparent atmospheric layer of roughly 3~km in thickness whose base lies at about 4~km above Pluto's surface, and where an average thermal gradient of about 5 K~km$^{-1}$ prevails. We discuss the possibility that small departures between the observed and modeled flash are caused by local topographic features (mountains) along Pluto's limb that block the stellar light. Finally, using two possible temperature profiles, and extrapolating our pressure profile from our deepest accessible level down to the surface, we obtain a possible range of 11.9-13.7~$\mu$bar for the surface pressure.

We present a science forecast for the eBOSS survey, part of the SDSS-IV project, which is a spectroscopic survey using multiple tracers of large-scale structure, including luminous red galaxies (LRGs), emission line galaxies (ELGs) and quasars (both as a direct probe of structure and through the Ly-$\alpha$ forest). Focusing on discrete tracers, we forecast the expected accuracy of the baryonic acoustic oscillation (BAO), the redshift-space distortion (RSD) measurements, the $f_{\rm NL}$ parameter quantifying the primordial non-Gaussianity, the dark energy and modified gravity parameters. We also use the line-of-sight clustering in the Ly-$\alpha$ forest to constrain the total neutrino mass. We find that eBOSS LRGs ($0.6<z<1.0$) (combined with the BOSS LRGs at $z>0.6$), ELGs ($0.6<z<1.2$) and Clustering Quasars (CQs) ($0.6<z<2.2$) can achieve a precision of 1%, 2.2% and 1.6% precisions, respectively, for spherically averaged BAO distance measurements. Using the same samples, the constraint on $f\sigma_8$ is expected to be 2.5%, 3.3% and 2.8% respectively. For primordial non-Gaussianity, eBOSS alone can reach an accuracy of $\sigma(f_{\rm NL})\sim10-15$, depending on the external measurement of the galaxy bias and our ability to model large-scale systematic errors. eBOSS can at most improve the dark energy Figure of Merit (FoM) by a factor of $3$ for the Chevallier-Polarski-Linder (CPL) parametrisation, and can well constrain three eigenmodes for the general equation-of-state parameter (Abridged).

We present a new selection technique of producing spectroscopic target catalogues for massive spectroscopic surveys for cosmology. This work was conducted in the context of the extended Baryon Oscillation Spectroscopic Survey (eBOSS), which will use ~200 000 emission line galaxies (ELGs) at 0.6<zspec<1.0 to obtain a precise baryon acoustic oscillation measurement. Our proposed selection technique is based on optical and near-infrared broad-band filter photometry. We used a training sample to define a quantity, the Fisher discriminant (linear combination of colours), which correlates best with the desired properties of the target: redshift and [OII] flux. The proposed selections are simply done by applying a cut on magnitudes and this Fisher discriminant. We used public data and dedicated SDSS spectroscopy to quantify the redshift distribution and [OII] flux of our ELG target selections. We demonstrate that two of our selections fulfil the initial eBOSS/ELG redshift requirements: for a target density of 180 deg^2, ~70% of the selected objects have 0.6<zspec<1.0 and only ~1% of those galaxies in the range 0.6<zspec<1.0 are expected to have a catastrophic zspec estimate. Additionally, the stacked spectra and stacked deep images for those two selections show characteristic features of star-forming galaxies. The proposed approach using the Fisher discriminant could, however, be used to efficiently select other galaxy populations, based on multi-band photometry, providing that spectroscopic information is available. This technique could thus be useful for other future massive spectroscopic surveys such as PFS, DESI, and 4MOST.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

As part of the Sloan Digital Sky Survey IV the extended Baryon Oscillation Spectroscopic Survey (eBOSS) will improve measurements of the cosmological distance scale by applying the Baryon Acoustic Oscillation (BAO) method to quasar samples. eBOSS will adopt two approaches to target quasars over 7500 sq. deg. First, a "CORE" quasar sample will combine optical selection in ugriz using a likelihood-based routine called XDQSOz, with a mid-IR-optical color-cut. eBOSS CORE selection (to g < 22 OR r < 22) should return ~ 70 quasars per sq. deg. at redshifts 0.9 < z < 2.2 and ~7 z > 2.1 quasars per sq. deg. Second, a selection based on variability in multi-epoch imaging from the Palomar Transient Factory should recover an additional ~3-4 z > 2.1 quasars per sq. deg. to g < 22.5. A linear model of how imaging systematics affect target density recovers the angular distribution of eBOSS CORE quasars over 96.7% (76.7%) of the SDSS North (South) Galactic Cap area. The eBOSS CORE quasar sample should thus be sufficiently dense and homogeneous over 0.9 < z < 2.2 to yield the first few-percent-level BAO constraint near z~1.5. eBOSS quasars at z > 2.1 will be used to improve BAO measurements in the Lyman-alpha Forest. Beyond its key cosmological goals, eBOSS should be the next-generation quasar survey, comprising > 500,000 new quasars and > 500,000 uniformly selected spectroscopically confirmed 0.9 < z < 2.2 quasars. At the conclusion of eBOSS, the SDSS will have provided unique spectra of over 800,000 quasars.

Using data from four deep fields (COSMOS, AEGIS, ECDFS, and CDFN), we study the correlation between the position of galaxies in the star formation rate (SFR) versus stellar mass plane and local environment at $z<1.1$. To accurately estimate the galaxy SFR, we use the deepest available Spitzer/MIPS 24 and Herschel/PACS datasets. We distinguish group environments ( $M_{halo}\sim$10$^{12.5-14.2}$$M_{\odot}$) based on the available deep X-ray data and lower halo mass environments based on the local galaxy density. We confirm that the Main Sequence (MS) of star forming galaxies is not a linear relation and there is a flattening towards higher stellar masses ( $M_*>10^{10.4-10.6}$ $M_{\odot}$), across all environments. At high redshift ( $0.5<z<1.1$ ), the MS varies little with environment. At low redshift ( $0.15<z<0.5$ ), group galaxies tend to deviate from the mean MS towards the region of quiescence with respect to isolated galaxies and less-dense environments. We find that the flattening of the MS toward low SFR is due to an increased fraction of bulge dominated galaxies at high masses. Instead, the deviation of group galaxies from the MS at low redshift is caused by a large fraction of red disk dominated galaxies which are not present in the lower density environments. Our results suggest that above a mass threshold ( $\sim10^{10.4}-10^{10.6}$$M_{\odot}$ ) stellar mass, morphology and environment act together in driving the evolution of the SF activity towards lower level. The presence of a dominating bulge and the associated quenching processes are already in place beyond $z\sim$1. The environmental effects appear, instead, at lower redshifts and have a long time-scale.

We present improved estimates of several global properties of the Milky Way, including its current star formation rate (SFR), the stellar mass contained in its disk and bulge+bar components, as well as its total stellar mass. We do so by combining previous measurements from the literature using a hierarchical Bayesian (HB) statistical method that allows us to account for the possibility that any value may be incorrect or have underestimated errors. We show that this method is robust to a wide variety of assumptions about the nature of problems in individual measurements or error estimates. Ultimately, our analysis yields a SFR for the Galaxy of $\dot{\mathrm{M}}_\star=1.65\pm0.19$ $\textrm{M}_\odot \textrm{yr}^{-1}$, assuming a Kroupa initial mass function (IMF). By combining HB methods with Monte Carlo simulations that incorporate the latest estimates of the Galactocentric radius of the Sun, $R_0$, the exponential scale length of the disk, $L_d$, and the local surface density of stellar mass, $\Sigma_\star(R_0)$, we show that the mass of the Galactic bulge+bar is $\textrm{M}_\star^B=0.91\pm0.07\times10^{10}$ $\textrm{M}_\odot$, the disk mass is $\textrm{M}_\star^D=5.17\pm1.11\times10^{10}$ $\textrm{M}_\odot$, and their combination yields a total stellar mass of $\textrm{M}_\star=6.08\pm1.14\times10^{10}$ $\textrm{M}_\odot$ (assuming a Kroupa IMF and an exponential disk profile). This analysis is based upon a new compilation of literature bulge mass estimates, normalized to common assumptions about the stellar initial mass function and Galactic disk properties, presented herein. We additionally find a bulge-to-total mass ratio for the Milky Way of $B/T=0.150^{+0.028}_{-0.019}$ and a specific star formation rate of $\dot{\mathrm{M}}_\star/\textrm{M}_\star=2.71\pm0.59\times10^{-11} \textrm{yr}^{-1}$.

We demonstrate a new statistical method of determining the global photometric properties of the Milky Way (MW) to an unprecedented degree of accuracy, allowing our Galaxy to be compared directly to objects measured in extragalactic surveys. Capitalizing on the high-quality imaging and spectroscopy dataset from the Sloan Digital Sky Survey (SDSS), we exploit the inherent dependence of galaxies' luminosities and colors on their total stellar mass, $\mathrm{M}_\star$, and star formation rate (SFR), $\mathrm{\dot{M}}_\star$, by selecting a sample of $Milky$ $Way$ $analog$ $galaxies$ designed to reproduce the best Galactic $\mathrm{M}_\star$ and $\mathrm{\dot{M}}_\star$ measurements, including all measurement uncertainties. Making the Copernican assumption that the MW is not extraordinary amongst galaxies of similar stellar mass and SFR, we then analyze the photometric properties of this matched sample, constraining the characteristics of our Galaxy without suffering interference from interstellar dust. We explore a variety of potential systematic errors that could affect this method, and find that they are subdominant to random uncertainties. We present both SDSS $ugriz$ absolute magnitudes and colors in both rest-frame $z$=0 and $z$=0.1 passbands for the MW, which are in agreement with previous estimates but can have up to $\sim$3$\times$ lower errors. We find the MW to have absolute magnitude $^0\!M_r-5\log h=-21.00_{-0.37}^{+0.38}$ and integrated color $^0(g-r)=0.682_{-0.056}^{+0.066}$, indicating that it may belong to the green-valley region in color-magnitude space and ranking it amongst the brightest and reddest of spiral galaxies. We also present new estimates of global stellar mass-to-light ratios for our Galaxy. This work will help relate our in-depth understanding of the Galaxy to studies of more distant objects.

We describe the algorithm used to select the Luminous Red Galaxy (LRG) sample for the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digital Sky Survey IV (SDSS-IV) using photometric data from both the SDSS and the Wide-Field Infrared Survey Explorer (WISE). LRG targets are required to meet a set of color selection criteria and have z-band and i-band MODEL magnitudes z < 19.95 and 19.9 < i < 21.8, respectively. Our algorithm selects roughly 50 LRG targets per square degree, the great majority of which lie in the redshift range 0.6 < z < 1.0 (median redshift 0.71). We demonstrate that our methods are highly effective at eliminating stellar contamination and lower-redshift galaxies. We perform a number of tests using spectroscopic data from SDSS-III/BOSS to determine the redshift reliability of our target selection and its ability to meet the science requirements of eBOSS. The SDSS spectra are of high enough signal-to-noise ratio that at least 89% of the target sample yields secure redshift measurements. We also present tests of the uniformity and homogeneity of the sample, demonstrating that it should be clean enough for studies of the large-scale structure of the universe at higher redshifts than SDSS-III/BOSS LRGs reached.

We present the rest-frame near-ultraviolet (NUV) spectroscopy of star-forming galaxies (SFGs) at 0.6<z<1.2 from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) in SDSS-IV. One of the eBOSS programs is to obtain 2 arcsec (about 15 kpc) fiber spectra of about 200,000 emission-line galaxies (ELGs) at redshift z>0.6. We use the data from the pilot observations of this program, including 8620 spectra of SFGs at 0.6<z<1.2. The median composite spectra of these SFGs at 2200 Ang < \lambda < 4000 Ang feature asymmetric, preferentially blueshifted non-resonant emission, Fe II*, and blueshifted resonant absorption, e.g., Fe II and Mg II, indicating ubiquitous outflows driven by star formation at these redshifts. For the absorption lines, we find a variety of velocity profiles with different degrees of blueshift. Comparing our new observations with the literature, we do not observe the non-resonant emission in the small-aperture (<40 pc) spectra of local star-forming regions with the Hubble Space Telescope, and find the observed line ratios in the SFG spectra to be different from those in the spectra of local star-forming regions, as well as those of quasar absorption-line systems in the same redshift range. We introduce an outflow model that can simultaneously explain the multiple observed properties and suggest that the variety of absorption velocity profiles and the line ratio differences are caused by scattered fluorescent emission filling in on top of the absorption in the large-aperture eBOSS spectra. We develop an observation-driven, model-independent method to correct the emission-infill to reveal the true absorption profiles. Our results show that eBOSS and future dark-energy surveys (e.g., DESI and PFS) will provide rich datasets of NUV spectroscopy for astrophysical applications.

We present a robust measurement and analysis of the rest-frame ultraviolet (UV) luminosity function at z=4-8. We use deep Hubble Space Telescope imaging over the CANDELS/GOODS fields, the Hubble Ultra Deep Field and the Year 1 Hubble Frontier Field deep parallel observations. These surveys provides an effective volume of 0.6-1.2 x 10^6 Mpc^3 over this epoch, allowing us to perform a robust search for faint (M_UV=-18) and bright (M_UV < -21) galaxies. We select candidate galaxies using a well-tested photometric redshift technique with careful screening of contaminants, finding a sample of 7446 galaxies at 3.5<z<8.5, with >1000 galaxies at z~6-8. We measure the luminosity function using a Markov Chain Monte Carlo analysis to measure robust uncertainties. At the faint end our results agree with previous studies, yet we find a higher abundance of UV-bright galaxies at z>6, with M* ~ -21 at z>5, different than that inferred based on previous trends at lower redshift. At z=8, a single power-law provides an equally good fit to the UV luminosity function, while at z=6 and 7, an exponential cutoff at the bright-end is moderately preferred. We compare to semi-analytical models, and find that the lack of evolution in M* is consistent with models where the impact of dust attenuation on the bright-end of the luminosity function decreases at higher redshift. We measure the evolution of the cosmic star-formation rate density, correcting for dust attenuation, and find that it declines as (1+z)^(-4.3 +/- 0.5) at z>4, consistent with observations at z>9. Our observations are consistent with a reionization history that starts at z>10, completes at z>6, and reaches a midpoint (x_HII = 0.5) at 6.7<z<9.4. Finally, our observations predict that the abundance of bright z=9 galaxies is likely higher than previous constraints, though consistent with recent estimates of bright z~10 galaxies. [abridged]

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

Using data from the All Wavelength Extended Groth Strip International Survey (AEGIS) we statistically detect the extended X-ray emission in the interstellar medium (ISM)/intra-cluster medium (ICM) in both active and normal galaxies at 0.3 <= z <= 1.3. For both active galactic nuclei (AGN) host galaxy and normal galaxy samples that are matched in restframe color, luminosity, and redshift distribution, we tentatively detect excess X-ray emission at scales of 1--10 arcsec at a few sigma significance in the surface brightness profiles. The exact significance of this detection is sensitive to the true characterization of Chandra's point spread function. The observed excess in the surface brightness profiles is suggestive of lower extended emission in AGN hosts compared to normal galaxies. This is qualitatively similar to theoretical predictions of the X-ray surface brightness profile from AGN feedback models, where feedback from AGN is likely to evacuate the gas from the center of the galaxy/cluster. We propose that AGN that are intrinsically under-luminous in X-rays, but have equivalent bolometric luminosities to our sources will be the ideal sample to study more robustly the effect of AGN feedback on diffuse ISM/ICM gas.

In this paper we perform a comprehensive study of the main sources of random and systematic errors in stellar mass measurement for galaxies using their Spectral Energy Distributions (SEDs). We use mock galaxy catalogs with simulated multi-waveband photometry (from U-band to mid-infrared) and known redshift, stellar mass, age and extinction for individual galaxies. Given different parameters affecting stellar mass measurement (photometric S/N ratios, SED fitting errors, systematic effects, the inherent degeneracies and correlated errors), we formulated different simulated galaxy catalogs to quantify these effects individually. We studied the sensitivity of stellar mass estimates to the codes/methods used, population synthesis models, star formation histories, nebular emission line contributions, photometric uncertainties, extinction and age. For each simulated galaxy, the difference between the input stellar masses and those estimated using different simulation catalogs, $\Delta\log(M)$, was calculated and used to identify the most fundamental parameters affecting stellar masses. We measured different components of the error budget, with the results listed as follows: (1). no significant bias was found among different codes/methods, with all having comparable scatter; (2). A source of error is found to be due to photometric uncertainties and low resolution in age and extinction grids; (3). The median of stellar masses among different methods provides a stable measure of the mass associated with any given galaxy; (4). The deviations in stellar mass strongly correlate with those in age, with a weaker correlation with extinction; (5). the scatter in the stellar masses due to free parameters are quantified, with the sensitivity of the stellar mass to both the population synthesis codes and inclusion of nebular emission lines studied.

Using data from the DEEP2 galaxy redshift survey and the All Wavelength Extended Groth Strip International Survey we obtain stacked X-ray maps of galaxies at 0.7 < z < 1.0 as a function of stellar mass. We compute the total X-ray counts of these galaxies and show that in the soft band (0.5--2,kev) there exists a significant correlation between galaxy X-ray counts and stellar mass at these redshifts. The best-fit relation between X-ray counts and stellar mass can be characterized by a power law with a slope of 0.58 +/- 0.1. We do not find any correlation between stellar mass and X-ray luminosities in the hard (2--7,kev) and ultra-hard (4--7,kev) bands. The derived hardness ratios of our galaxies suggest that the X-ray emission is degenerate between two spectral models, namely point-like power-law emission and extended plasma emission in the interstellar medium. This is similar to what has been observed in low redshift galaxies. Using a simple spectral model where half of the emission comes from power-law sources and the other half from the extended hot halo we derive the X-ray luminosities of our galaxies. The soft X-ray luminosities of our galaxies lie in the range 10^39-8x10^40, ergs/s. Dividing our galaxy sample by the criteria U-B > 1, we find no evidence that our results for X-ray scaling relations depend on optical color.

We present the discovery of two galaxy overdensities in the Hubble Space Telescope UDF: a proto-cluster, HUDFJ0332.4-2746.6 at $z = 1.84 \pm 0.01$, and a group, HUDFJ0332.5-2747.3 at $z =1.90 \pm 0.01$. Assuming viralization, the velocity dispersion of HUDFJ0332.4-2746.6 implies a mass of $M_{200}= (2.2 \pm 1.8) \times 10^{14} M_{\odot}$, consistent with the lack of extended X-ray emission. Neither overdensity shows evidence of a red sequence. About $50\%$ of their members show interactions and/or disturbed morphologies, which are signatures of merger remnants or disk instability. Most of their ETGs have blue colors and show recent star-formation. These observations reveal for the first time large fractions of spectroscopically confirmed star-forming blue ETGs in proto-clusters at $z\approx 2$. These star-forming ETGs are most likely among the progenitors of the quiescent population in clusters at more recent epochs. Their mass-size relation is consistent with that of passive ETGs in clusters at $z\sim0.7-1.5$. If these galaxies are the progenitors of cluster ETGs at these lower redshifts, their size would evolve according to a similar mass-size relation. It is noteworthy that quiescent ETGs in clusters at $z=1.8-2$ also do not show any significant size evolution over this redshift range, contrary to field ETGs. The ETG fraction is $\lesssim 50\%$, compared to the typical quiescent ETG fraction of $\approx 80\%$ in cluster cores at $z< 1$. The fraction, masses, and colors of the newly discovered ETGs imply that other cluster ETGs will be formed/accreted at a later time.

We describe a new method of combining optical and infrared photometry to select Luminous Red Galaxies (LRGs) at redshifts $z > 0.6$. We explore this technique using a combination of optical photometry from CFHTLS and HST, infrared photometry from the WISE satellite, and spectroscopic or photometric redshifts from the DEEP2 Galaxy Redshift Survey or COSMOS. We present a variety of methods for testing the success of our selection, and present methods for optimization given a set of rest-frame color and redshift requirements. We have tested this selection in two different regions of the sky, the COSMOS and Extended Groth Strip (EGS) fields, to reduce the effect of cosmic/sample variance. We have used these methods to assemble large samples of LRGs for two different ancillary programs as a part of the SDSS-III/ BOSS spectroscopic survey. This technique is now being used to select $\sim$600,000 LRG targets for SDSS-IV/eBOSS, which began observations in Fall 2014, and will be adapted for the proposed DESI survey. We have found these methods can select high-redshift LRGs efficiently with minimal stellar contamination; this is extremely difficult to achieve with selections that rely on optical photometry alone.

We present the results of deep \chandra\ imaging of the central region of the Extended Groth Strip, the AEGIS-X Deep (AEGIS-XD) survey. When combined with previous \chandra\ observations of a wider area of the strip, AEGIS-X Wide (AEGIS-XW; Laird et~al. 2009), these provide data to a nominal exposure depth of 800ks in the three central ACIS-I fields, a region of approximately $0.29$~deg$^{2}$. This is currently the third deepest X-ray survey in existence, a factor $\sim 2-3$ shallower than the Chandra Deep Fields (CDFs) but over an area $\sim 3$ times greater than each CDF. We present a catalogue of 937 point sources detected in the deep \chandra\ observations. We present identifications of our X-ray sources from deep ground-based, Spitzer, GALEX and HST imaging. Using a likelihood ratio analysis, we associate multi band counterparts for 929/937 of our X-ray sources, with an estimated 95~\% reliability, making the identification completeness approximately 94~\% in a statistical sense. Reliable spectroscopic redshifts for 353 of our X-ray sources are provided predominantly from Keck (DEEP2/3) and MMT Hectospec, so the current spectroscopic completeness is $\sim 38$~per cent. For the remainder of the X-ray sources, we compute photometric redshifts based on multi-band photometry in up to 35 bands from the UV to mid-IR. Particular attention is given to the fact that the vast majority the X-ray sources are AGN and require hybrid templates. Our photometric redshifts have mean accuracy of $\sigma=0.04$ and an outlier fraction of approximately 5\%, reaching $\sigma=0.03$ with less than 4\% outliers in the area covered by CANDELS . The X-ray, multi-wavelength photometry and redshift catalogues are made publicly available.

In order to derive the AGN contribution to the cosmological ionizing emissivity we have selected faint AGN candidates at $z>4$ in the CANDELS GOODS-South field which is one of the deepest fields with extensive multiwavelength coverage from Chandra, HST, Spitzer and various groundbased telescopes. We have adopted a relatively novel criterion. As a first step high redshift galaxies are selected in the NIR $H$ band down to very faint levels ($H\leq27$) using reliable photometric redshifts. This corresponds at $z>4$ to a selection criterion based on the galaxy rest-frame UV flux. AGN candidates are then picked up from this parent sample if they show X-ray fluxes above a threshold of $F_X\sim 1.5\times 10^{-17}$ cgs (0.5-2 keV). We have found 22 AGN candidates at $z>4$ and we have derived the first estimate of the UV luminosity function in the redshift interval $4<z<6.5$ and absolute magnitude interval $-22.5\lesssim M_{1450} \lesssim -18.5$ typical of local Seyfert galaxies. The faint end of the derived luminosity function is about two/four magnitudes fainter at $z\sim 4-6$ than that derived from previous UV surveys. We have then estimated ionizing emissivities and hydrogen photoionization rates in the same redshift interval under reasonable assumptions and after discussion of possible caveats, the most important being the large uncertainties involved in the estimate of photometric redshift for sources with featureless, almost power-law SEDs and/or low average escape fraction of ionizing photons from the AGN host galaxies. We argue that, under reasonable evaluations of possible biases, the probed AGN population can produce at $z=4-6.5$ photoionization rates consistent with that required to keep highly ionized the intergalactic medium observed in the Lyman-$\alpha$ forest of high redshift QSO spectra, providing an important contribution to the cosmic reionization.

The focus of this report is on the opportunities enabled by the combination of LSST, Euclid and WFIRST, the optical surveys that will be an essential part of the next decade's astronomy. The sum of these surveys has the potential to be significantly greater than the contributions of the individual parts. As is detailed in this report, the combination of these surveys should give us multi-wavelength high-resolution images of galaxies and broadband data covering much of the stellar energy spectrum. These stellar and galactic data have the potential of yielding new insights into topics ranging from the formation history of the Milky Way to the mass of the neutrino. However, enabling the astronomy community to fully exploit this multi-instrument data set is a challenging technical task: for much of the science, we will need to combine the photometry across multiple wavelengths with varying spectral and spatial resolution. We identify some of the key science enabled by the combined surveys and the key technical challenges in achieving the synergies.

We present the public release of the stellar mass catalogs for the GOODS-S and UDS fields obtained using some of the deepest near-IR images available, achieved as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) project. We combine the effort from ten different teams, who computed the stellar masses using the same photometry and the same redshifts. Each team adopted their preferred fitting code, assumptions, priors, and parameter grid. The combination of results using the same underlying stellar isochrones reduces the systematics associated with the fitting code and other choices. Thanks to the availability of different estimates, we can test the effect of some specific parameters and assumptions on the stellar mass estimate. The choice of the stellar isochrone library turns out to have the largest effect on the galaxy stellar mass estimates, resulting in the largest distributions around the median value (with a semi interquartile range larger than 0.1 dex). On the other hand, for most galaxies, the stellar mass estimates are relatively insensitive to the different parameterizations of the star formation history. The inclusion of nebular emission in the model spectra does not have a significant impact for the majority of galaxies (less than a factor of 2 for ~80% of the sample). Nevertheless, the stellar mass for the subsample of young galaxies (age < 100 Myr), especially in particular redshift ranges (e.g., 2.2 < z < 2.4, 3.2 < z < 3.6, and 5.5 < z < 6.5), can be seriously overestimated (by up to a factor of 10 for < 20 Myr sources) if nebular contribution is ignored.

We provide an overview of the science benefits of combining information from the Square Kilometre Array (SKA) and the Large Synoptic Survey Telescope (LSST). We first summarise the capabilities and timeline of the LSST and overview its science goals. We then discuss the science questions in common between the two projects, and how they can be best addressed by combining the data from both telescopes. We describe how weak gravitational lensing and galaxy clustering studies with LSST and SKA can provide improved constraints on the causes of the cosmological acceleration. We summarise the benefits to galaxy evolution studies of combining deep optical multi-band imaging with radio observations. Finally, we discuss the excellent match between one of the most unique features of the LSST, its temporal cadence in the optical waveband, and the time resolution of the SKA.

Most massive passive galaxies are compact at high redshifts, but similarly compact massive galaxies are rare in the local universe. The most common interpretation of this phenomenon is that massive galaxies have grown in size by a factor of about five since redshift z=2. An alternative explanation is that recently quenched massive galaxies are larger (a "progenitor bias"). In this paper we explore the importance of progenitor bias by looking for systematic differences in the stellar populations of compact early-type galaxies in the DEEP2 survey as a function of size. Our analysis is based on applying the statistical technique of bootstrap resampling to constrain differences in the median ages of our samples and to begin to characterize the distribution of stellar populations in our co-added spectra. The light-weighted ages of compact early-type galaxies at redshifts 0.5 < z < 1.4 are compared to those of a control sample of larger galaxies at similar redshifts. We find that massive compact early-type galaxies selected on the basis of red color and high bulge-to-total ratio are younger than similarly selected larger galaxies, suggesting that size growth in these objects is not driven mainly by progenitor bias, and that individual galaxies grow as their stellar populations age. However, compact early-type galaxies selected on the basis of image smoothness and high bulge-to-total ratio are older than a control sample of larger galaxies. Progenitor bias will play a significant role in defining the apparent size changes of early-type galaxies if they are selected on the basis of the smoothness of their light distributions.

This white paper summarizes the conclusions of the Snowmass White Paper "Spectroscopic Needs for Imaging Dark Energy Experiments" (arXiv:1309.5384) which are relevant to the calibration of LSST photometric redshifts; i.e., the accurate characterization of biases and uncertainties in photo-z's. Any significant miscalibration will lead to systematic errors in photo-z's, impacting nearly all extragalactic science with LSST. As existing deep redshift samples have failed to yield highly-secure redshifts for a systematic 20%-60% of their targets, it is a strong possibility that future deep spectroscopic samples will not solve the calibration problem on their own. The best options in this scenario are provided by cross-correlation methods that utilize clustering with objects from spectroscopic surveys (which need not be fully representative) to trace the redshift distribution of the full sample. For spectroscopy, the eBOSS survey would enable a basic calibration of LSST photometric redshifts, while the expected LSST-DESI overlap would be more than sufficient for an accurate calibration at z>0.2. A DESI survey of nearby galaxies conducted in bright time would enable accurate calibration down to z~0. The expanded areal coverage provided by the transfer of the DESI instrument (or duplication of it at the Blanco Telescope) would enable the best possible calibration from cross-correlations, in addition to other science gains.

This white paper summarizes those conclusions of the Snowmass White Paper "Spectroscopic Needs for Imaging Dark Energy Experiments" (arXiv:1309.5384) which are relevant to the training of LSST photometric redshifts; i.e., the use of spectroscopic redshifts to improve algorithms and reduce photo-z errors. The larger and more complete the available training set is, the smaller the RMS error in photo-z estimates should be, increasing LSST's constraining power. Among the better US-based options for this work are the proposed MANIFEST fiber feed for the Giant Magellan Telescope or (with lower survey speed) the WFOS spectrograph on the Thirty Meter Telescope (TMT). Due to its larger field of view and higher multiplexing, the PFS spectrograph on Subaru would be able to obtain a baseline training sample faster than TMT; comparable performance could be achieved with a highly-multiplexed spectrograph on Gemini with at least a 20 arcmin diameter field of view.