The Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) program monitors 849 active galactic nuclei (AGN) both spectroscopically and photometrically. The photometric observations used in this work span over four years and provide an excellent baseline for variability studies of these objects. We present the photometric light curves from 2014 to 2017 obtained by the Steward Observatory's Bok telescope and the CFHT telescope with MegaCam. We provide details on the data acquisition and processing of the data from each telescope, the difference imaging photometry used to produce the light curves, and the calculation of a variability index to quantify each AGN's variability. We find that the Welch-Stetson J-index provides a useful characterization of AGN variability and can be used to select AGNs for further study.

It is well known that reverberation mapping of active galactic nuclei (AGN) reveals a relationship between AGN luminosity and the size of the broad-line region, and that use of this relationship, combined with the Doppler width of the broad emission line, enables an estimate of the mass of the black hole at the center of the active nucleus based on a single spectrum. This has been discussed in numerous papers over the last two decades. An unresolved key issue is the choice of parameter used to characterize the line width; generally, most researchers use FWHM in favor of line dispersion (the square root of the second moment of the line profile) because the former is easier to measure, less sensitive to blending with other features, and usually can be measured with greater precision. However, use of FWHM introduces a bias, stretching the mass scale such that high masses are overestimated and low masses are underestimated. Here we investigate estimation of black hole masses in AGN based on individual or "single epoch" observations, with a particular emphasis in comparing mass estimates based on line dispersion and FWHM. We confirm the recent findings that, in addition to luminosity and line width, a third parameter is required to obtain accurate masses and that parameter seems to be Eddington ratio. We present simplified empirical formulae for estimating black hole masses from the Hbeta (4861 A) and C IV (1549 A) emission lines.

Quasar broad emission lines are largely powered by photoionization from the accretion continuum. Increased central luminosity will enhance line emissivity in more distant clouds, leading to increased average distance of the broad-line-emitting clouds and decreased averaged line width, known as the broad-line region (BLR) "breathing". However, different lines breathe differently, and some high-ionization lines, such as C IV, can even show "anti-breathing" where the line broadens when luminosity increases. Using multi-year photometric and spectroscopic monitoring data from the Sloan Digital Sky Survey Reverberation Mapping project, we quantify the breathing effect ($\Delta$log W=$\alpha\Delta$log L) of broad H$\alpha$, H$\beta$, Mg II, C IV, and C III] for statistical quasar samples over $z\approx 0.1-2.5$. We found that H$\beta$ displays the most consistent normal breathing expected from the virial relation ($\alpha\sim-0.25$), Mg II and H$\alpha$ on average show no breathing ($\alpha\sim 0$), and C IV (and similarly C III] and Si IV mostly shows anti-breathing ($\alpha>0$). The anti-breathing of C IV can be well understood by the presence of a non-varying core component in addition to a reverberating broad-base component, consistent with earlier findings. The deviation from canonical breathing introduces extra scatter (a luminosity-dependent bias) in single-epoch virial BH mass estimates due to intrinsic quasar variability, which underlies the long argued caveats of C IV single-epoch masses. Using the line dispersion instead of FWHM leads to less, albeit still substantial, deviations from canonical breathing in most cases. Our results strengthen the need for reverberation mapping to provide reliable quasar BH masses, and quantify the level of variability-induced bias in single-epoch BH masses based on various lines.

We present the results of a pilot Hubble Space Telescope (HST) imaging study of the host galaxies of ten quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Probing more than an order of magnitude in BH and stellar masses, our sample is the first statistical sample to study the BH-host correlations beyond z>0.3 with reliable BH masses from reverberation mapping rather than from single-epoch spectroscopy. We perform image decomposition in two HST bands (UVIS-F606W and IR-F110W) to measure host colors and estimate stellar masses using empirical relations between broad-band colors and the mass-to-light ratio. The stellar masses of our targets are mostly dominated by a bulge component. The BH masses and stellar masses of our sample broadly follow the same correlations found for local RM AGN and quiescent bulge-dominant galaxies, with no strong evidence of evolution in the M_BH-M_*,bulge relation to z~0.6. We further compare the host light fraction from HST imaging decomposition to that estimated from spectral decomposition. We found a good correlation between the host fractions derived with both methods. However, the host fraction derived from spectral decomposition is systematically smaller than that from imaging decomposition by ~30%, indicating different systematics in both approaches. This study paves the way for upcoming more ambitious host galaxy studies of quasars with direct RM-based BH masses at high redshift.

The origin of the correlations between mass, morphology, quenched fraction, and formation history in galaxies is difficult to define, primarily due to the uncertainties in galaxy star-formation histories. Star-formation histories are better constrained for higher redshift galaxies, observed closer to their formation and quenching epochs. Here we use "non-parametric" star-formation histories and a nested sampling method to derive constraints on the formation and quenching timescales of quiescent galaxies at $0.7<z<2.5$. We model deep HST grism spectroscopy and photometry from the CLEAR (CANDELS Lyman$-\alpha$ Emission at Reionization) survey. The galaxy formation redshifts, $z_{50}$ (defined as the point where they had formed 50\% of their stellar mass) range from $z_{50}\sim 2$ (shortly prior to the observed epoch) up to $z_{50} \simeq 5-8$. \editone{We find that early formation redshifts are correlated with high stellar-mass surface densities, $\log \Sigma_1 / (M_\odot\ \mathrm{kpc}^{-2}) >$10.25, where $\Sigma_1$ is the stellar mass within 1~pkpc (proper kpc). Quiescent galaxies with the highest stellar-mass surface density, $\log\Sigma_1 / (M_\odot\ \mathrm{kpc}^{-2}) > 10.25$, } show a \textit{minimum} formation redshift: all such objects in our sample have $z_{50} > 2.9$. Quiescent galaxies with lower surface density, $\log \Sigma_1 / (M_\odot\ \mathrm{kpc}^{-2}) = 9.5 - 10.25$, show a range of formation epochs ($z_{50} \simeq 1.5 - 8$), implying these galaxies experienced a range of formation and assembly histories. We argue that the surface density threshold $\log\Sigma_1/(M_\odot\ \mathrm{kpc}^{-2})>10.25$ uniquely identifies galaxies that formed in the first few Gyr after the Big Bang, and we discuss the implications this has for galaxy formation models.

We present reverberation mapping results for the MgII 2800 A broad emission line in a sample of 193 quasars at 0.35<z<1.7 with photometric and spectroscopic monitoring observations from the Sloan Digital Sky Survey Reverberation Mapping project during 2014 - 2017. We find significant time lags between the MgII and continuum lightcurves for 57 quasars and define a "gold sample" of 24 quasars with the most reliable lag measurements. We estimate false-positive rates for each lag that range from 1-24%, with an average false-positive rate of 11% for the full sample and 8% for the gold sample. There are an additional ~40 quasars with marginal MgII lag detections which may yield reliable lags after additional years of monitoring. The MgII lags follow a radius -- luminosity relation with a best-fit slope that is consistent with alpha=0.5 but with an intrinsic scatter of 0.36dex that is significantly larger than found for the Hb radius -- luminosity relation. For targets with SDSS-RM lag measurements of other emission lines, we find that our MgII lags are similar to the Hb lags and ~2-3 times larger than the CIV lags. This work significantly increases the number of MgII broad-line lags and provides additional reverberation-mapped black hole masses, filling the redshift gap at the peak of supermassive black hole growth between the Hb and CIV emission lines in optical spectroscopy.

Results from a few decades of reverberation mapping (RM) studies have revealed a correlation between the radius of the broad-line emitting region (BLR) and the continuum luminosity of active galactic nuclei. This "radius-luminosity" relation enables survey-scale black-hole mass estimates across cosmic time, using relatively inexpensive single-epoch spectroscopy, rather than intensive RM time monitoring. However, recent results from newer reverberation mapping campaigns challenge this widely used paradigm, reporting quasar BLR sizes that differ significantly from the previously established radius-luminosity relation. Using simulations of the radius--luminosity relation with the observational parameters of SDSS-RM, we find that this difference is not likely due to observational biases. Instead, it appears that previous RM samples were biased to a subset of quasar properties, and the broader parameter space occupied by the SDSS-RM quasar sample has a genuinely wider range of BLR sizes. We examine the correlation between the deviations from the radius-luminosity relation and several quasar parameters; the most significant correlations indicate that the deviations depend on UV/optical SED and the relative amount of ionizing radiation. Our results indicate that single-epoch black-hole mass estimates that do not account for the diversity of quasars in the radius-luminosity relation could be overestimated by an average of ~0.3 dex.

Active galactic nuclei (AGNs) have long been observed to "twinkle" (i.e., their brightness varies with time) on timescales from days to years in the UV/optical bands. Such AGN UV/optical variability is essential for probing the physics of supermassive black holes (SMBHs), the accretion disk, and the broad-line region. Here we show that the temperature fluctuations of an AGN accretion disk, which is magnetically coupled with the corona, can account for observed high-quality AGN optical light curves. We calculate the temperature fluctuations by considering the gas physics of the accreted matter near the SMBH. We find that the resulting simulated AGN UV/optical light curves share the same statistical properties as the observed ones as long as the dimensionless viscosity parameter $\alpha$, which is widely believed to be controlled by magnetohydrodynamic (MHD) turbulence in the accretion disk, is about $0.01$---$0.2$. Moreover, our model can simultaneously explain the larger-than-expected accretion disk sizes and the dependence of UV/optical variability upon wavelength for NGC 5548. Our model also has the potential to explain some other observational facts of AGN UV/optical variability, including the timescale-dependent bluer-when-brighter color variability and the dependence of UV/optical variability on AGN luminosity and black hole mass. Our results also demonstrate a promising way to infer the black-hole mass, the accretion rate, and the radiative efficiency, thereby facilitating understanding of the gas physics and MHD turbulence near the SMBH and its cosmic mass growth history by fitting the AGN UV/optical light curves in the era of time-domain astronomy.

We investigate the performance of different methodologies that measure the time lag between broad-line and continuum variations in reverberation mapping data using simulated light curves that probe a range of cadence, time baseline, and signal-to-noise ratio in the flux measurements. We compare three widely-adopted lag measuring methods: the Interpolated Cross-Correlation Function (ICCF), the z-transformed Discrete Correlation Function (ZDCF) and the MCMC code JAVELIN, for mock data with qualities typical of multi-object spectroscopic reverberation mapping (MOS-RM) surveys that simultaneously monitor hundreds of quasars. We quantify the overall lag detection efficiency, the rate of false detections, and the quality of lag measurements for each of these methods and under different survey designs (e.g., observing cadence and depth) using mock quasar light curves. Overall JAVELIN and ICCF outperform ZDCF in essentially all tests performed. Compared with ICCF, JAVELIN produces higher quality lag measurements, is capable of measuring more lags with timescales shorter than the observing cadence, is less susceptible to seasonal gaps and S/N degradation in the light curves, and produces more accurate lag uncertainties. We measure the Hbeta broad-line region size-luminosity (R-L) relation with each method using the simulated light curves to assess the impact of selection effects of the design of MOS-RM surveys. The slope of the R-L relation measured by JAVELIN is the least biased among the three methods, and is consistent across different survey designs. These results demonstrate a clear preference for JAVELIN over the other two non-parametric methods for MOS-RM programs, particularly in the regime of limited light curve quality as expected from most MOS-RM programs.

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.

We investigate the effects of extended multi-year light curves (9-year photometry and 5-year spectroscopy) on the detection of time lags between the continuum variability and broad-line response of quasars at z>~1.5, and compare with the results using 4-year photometry+spectroscopy presented in a companion paper. We demonstrate the benefits of the extended light curves in three cases: (1) lags that are too long to be detected by the shorter-duration data but can be detected with the extended data; (2) lags that are recovered by the extended light curves but missed in the shorter-duration data due to insufficient light curve quality; and (3) lags for different broad line species in the same object. These examples demonstrate the importance of long-term monitoring for reverberation mapping to detect lags for luminous quasars at high-redshift, and the expected performance of the final dataset from the Sloan Digital Sky Survey Reverberation Mapping project that will have 11-year photometric and 7-year spectroscopic baselines.

We present accretion-disk structure measurements from continuum lags in the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Lags are measured using the \texttt{JAVELIN} software from the first-year SDSS-RM $g$ and $i$ photometry, resulting in well-defined lags for 95 quasars, 33 of which have lag SNR $>$ 2$\sigma$. We also estimate lags using the \texttt{ICCF} software and find consistent results, though with larger uncertainties. Accretion-disk structure is fit using a Markov Chain Monte Carlo approach, parameterizing the measured continuum lags as a function of disk size normalization, wavelength, black hole mass, and luminosity. In contrast with previous observations, our best-fit disk sizes and color profiles are consistent (within 1.5~$\sigma$) with the \citet{SS73} analytic solution. We also find that more massive quasars have larger accretion disks, similarly consistent with the analytic accretion-disk model. The data are inconclusive on a correlation between disk size and continuum luminosity, with results that are consistent with both no correlation and with the \citet{SS73} expectation. The continuum lag fits have a large excess dispersion, indicating that our measured lag errors are underestimated and/or our best-fit model may be missing the effects of orientation, spin, and/or radiative efficiency. We demonstrate that fitting disk parameters using only the highest-SNR lag measurements biases best-fit disk sizes to be larger than the disk sizes recovered using a Bayesian approach on the full sample of well-defined lags.

The width of the broad emission lines in quasars is commonly characterized either by the full-width-at-half-maximum (FWHM) or the square root of the second moment of the line profile ($\sigma_{\rm line}$), and used as an indicator of the virial velocity of the broad-line region (BLR) in the estimation of black hole (BH) mass. We measure FWHM and $\sigma_{\rm line}$ for H$\alpha$, H$\beta$ and Mg II broad lines in both the mean and root-mean-square (rms) spectra of a large sample of quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. We introduce a new quantitative recipe to measure $\sigma_{\rm line}$ that is reproducible, less susceptible to noise and blending in the wings, and scales with the intrinsic width of the line. We compare the four definitions of line width (FWHM and $\sigma_{\rm line}$ in mean and rms spectra, respectively) for each of the three broad lines and among different lines. There are strong correlations among different width definitions for each line, providing justification for using the line width measured in single-epoch spectroscopy as a virial velocity indicator. There are also strong correlations among different lines, suggesting alternative lines to H$\beta$ can be used to estimate virial BH masses. We further investigate the correlations between virial BH masses using different line width definitions and the stellar velocity dispersion of the host galaxies, and the dependence of line shape (characterized by the ratio FWHM/$\sigma_{\rm line}$) on physical properties of the quasar. Our results provide further evidence that FWHM is more sensitive to the orientation of a flattened BLR geometry than $\sigma_{\rm line}$, but the overall comparison between the virial BH mass and host stellar velocity dispersion does not provide conclusive evidence that one particular width definition is significantly better than the others.

We analyze extensive spectroscopic and photometric data of the hypervariable quasar SDSS J131424+530527 (RMID 017) at z=0.456, an optical "changing look" quasar from the Sloan Digital Sky Survey Reverberation Mapping project that increased in optical luminosity by a factor of 10 between 2014 and 2017. The observed broad emission lines all respond in luminosity and width to the changing optical continuum, as expected for photoionization in a stratified, virialized broad emission line region. The luminosity changes therefore result from intrinsic changes in accretion power rather than variable obscuration. The variability is continuous and apparently stochastic, disfavoring an origin as a discrete event such as a tidal disruption flare or microlensing event. It is coordinated on day timescales with blue leading red, consistent with reprocessing powering the entire optical SED. We show that this process cannot work in a standard thin disk geometry on energetic grounds, and would instead require a large covering factor reprocessor. Disk instability models could potentially also explain the data, provided that the instability sets in near the inner radius of a geometrically thick accretion disk.

We present reverberation-mapping lags and black-hole mass measurements using the CIV 1549 broad emission line from a sample of 349 quasars monitored as a part of the Sloan Digital Sky Survey Reverberation Mapping Project. Our data span four years of spectroscopic and photometric monitoring for a total baseline of 1300 days. We report significant time delays between the continuum and the CIV 1549 emission line in 52 quasars, with an estimated false-positive detection rate of 10%. Our analysis of marginal lag measurements indicates that there are on the order of 100 additional lags that should be recoverable by adding more years of data from the program. We use our measurements to calculate black-hole masses and fit an updated CIV radius-luminosity relationship. Our results significantly increase the sample of quasars with CIV RM results, with the quasars spanning two orders of magnitude in luminosity toward the high-luminosity end of the CIV radius-luminosity relation. In addition, these quasars are located at among the highest redshifts (z~1.4-2.8) of quasars with black hole masses measured with reverberation mapping. This work constitutes the first large sample of CIV reverberation-mapping measurements in more than a dozen quasars, demonstrating the utility of multi-object reverberation mapping campaigns.

We present a detailed characterization of the 849 broad-line quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Our quasar sample covers a redshift range of 0.1<z<4.5 and is flux-limited to i_PSF<21.7 without any other cuts on quasar properties. The main sample characterization includes: 1) spectral measurements of the continuum and broad emission lines for individual objects from the coadded first-season spectroscopy in 2014; 2) identification of broad and narrow absorption lines in the spectra; 3) optical variability properties for continuum and broad lines from multi-epoch spectroscopy. We provide improved systemic redshift estimates for all quasars, and demonstrate the effects of signal-to-noise ratio on the spectral measurements. We compile measured properties for all 849 quasars along with supplemental multi-wavelength data for subsets of our sample from other surveys. The SDSS-RM sample probes a diverse range in quasar properties, and shows well detected continuum and broad-line variability for many objects from first-season monitoring data. The compiled properties serve as the benchmark for follow-up work based on SDSS-RM data. The spectral fitting tools are made public along with this work.

This white paper specifies the footprints, cadence requirements, and total-depth requirements needed to allow the most-successful AGN studies in the four currently selected LSST Deep-Drilling Fields (DDFs): ELAIS-S1, XMM-LSS, CDF-S, and COSMOS. The information provided on cadence and total-depth requirements will also likely be applicable to enabling effective AGN science in any additional DDFs that are chosen.

We systematically investigate short-timescale ($<$10-day rest-frame) CIV broad absorption-line (BAL) variability to constrain quasar-wind properties and provide insights into BAL-variability mechanisms in quasars. We employ data taken by the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, as the rapid cadence of these observations provides a novel opportunity to probe BAL variability on shorter rest-frame timescales than have previously been explored. In a sample of 27 quasars with a median of 58 spectral epochs per quasar, we have identified 15 quasars ($55^{+18}_{-14}$%), 19 of 37 CIV BAL troughs ($51^{+15}_{-12}$%), and 54 of 1460 epoch pairs ($3.7 \pm 0.5$%) that exhibit significant CIV BAL equivalent-width variability on timescales of less than 10 days in the quasar rest frame. These frequencies indicate that such variability is common among quasars and BALs, though somewhat rare among epoch pairs. Thus, models describing BALs and their behavior must account for variability on timescales down to less than a day in the quasar rest frame. We also examine a variety of spectral characteristics and find that in some cases, BAL variability is best described by ionization-state changes, while other cases are more consistent with changes in covering fraction or column density. We adopt a simple model to constrain the density and radial distance of two outflows appearing to vary by ionization-state changes, yielding outflow density lower limits consistent with previous work.

We present reverberation mapping results from the first year of combined spectroscopic and photometric observations of the Sloan Digital Sky Survey Reverberation Mapping Project. We successfully recover reverberation time delays between the $g+i$-band emission and the broad H$\beta$ emission line for a total of 44 quasars, and for the broad H$\alpha$ emission line in 18 quasars. Time delays are computed using the JAVELIN and CREAM software and the traditional interpolated cross-correlation function (ICCF): Using well defined criteria, we report measurements of 32 H$\beta$ and 13 H$\alpha$ lags with JAVELIN, 42 H$\beta$ and 17 H$\alpha$ lags with CREAM, and 16 H$\beta$ and 8 H$\alpha$ lags with the ICCF. Lag values are generally consistent among the three methods, though we typically measure smaller uncertainties with JAVELIN and CREAM than with the ICCF, given the more physically motivated light curve interpolation and more robust statistical modeling of the former two methods. The median redshift of our H$\beta$-detected sample of quasars is 0.53, significantly higher than that of the previous reverberation-mapping sample. We find that in most objects, the time delay of the H$\alpha$ emission is consistent with or slightly longer than that of H$\beta$. We measure black hole masses using our measured time delays and line widths for these quasars. These black hole mass measurements are mostly consistent with expectations based on the local M-sigma relationship, and are also consistent with single-epoch black hole mass measurements. This work increases the current sample size of reverberation-mapped active galaxies by about two-thirds and represents the first large sample of reverberation mapping observations beyond the local universe (z < 0.3).

Recent multi-band variability studies have revealed that active galactic nucleus (AGN) accretion disc sizes are generally larger than the predictions of the classical thin disc by a factor of $2\sim 3$. This hints at some missing key ingredient in the classical thin disc theory: here, we propose an accretion disc wind. For a given bolometric luminosity, in the outer part of an accretion disc, the effective temperature in the wind case is higher than that in the no-wind one; meanwhile, the radial temperature profile of the wind case is shallower than the no-wind one. In presence of winds, for a given band, blackbody emission from large radii can contribute more to the observed luminosity than the no-wind case. Therefore, the disc sizes of the wind case can be larger than those of the no-wind case. We demonstrate that a model with the accretion rate scaling as $\dot{M}_0 (R/R_{\mathrm{S}})^{\beta}$ (i.e., the accretion rate declines with decreasing radius due to winds) can match both the inter-band time lags and the spectral energy distribution of NGC 5548. Our model can also explain the inter-band time lags of other sources. Therefore, our model can help decipher current and future continuum reverberation mapping observations.

We use deep \textit{Hubble Space Telescope} spectroscopy to constrain the metallicities and (\editone{light-weighted}) ages of massive ($\log M_\ast/M_\odot\gtrsim10$) galaxies selected to have quiescent stellar populations at $1.0<z<1.8$. The data include 12--orbit depth coverage with the WFC3/G102 grism covering $\sim$ $8,000<\lambda<11,500$~\AA\, at a spectral resolution of $R\sim 210$ taken as part of the CANDELS Lyman-$\alpha$ Emission at Reionization (CLEAR) survey. At $1.0<z<1.8$, the spectra cover important stellar population features in the rest-frame optical. We simulate a suite of stellar population models at the grism resolution, fit these to the data for each galaxy, and derive posterior likelihood distributions for metallicity and age. We stack the posteriors for subgroups of galaxies in different redshift ranges that include different combinations of stellar absorption features. Our results give \editone{light-weighted ages of $t_{z \sim 1.1}= 3.2\pm 0.7$~Gyr, $t_{z \sim 1.2}= 2.2\pm 0.6$~Gyr, $t_{z\sim1.3}= 3.1\pm 0.6$~Gyr, and $t_{z\sim1.6}= 2.0 \pm 0.6$~Gyr, \editone{for galaxies at $z\sim 1.1$, 1.2, 1.3, and 1.6. This} implies that most of the massive quiescent galaxies at $1<z<1.8$ had formed $>68$\% of their stellar mass by a redshift of $z>2$}. The posteriors give metallicities of \editone{$Z_{z\sim1.1}=1.16 \pm 0.29$~$Z_\odot$, $Z_{z\sim1.2}=1.05 \pm 0.34$~$Z_\odot$, $Z_{z\sim1.3}=1.00 \pm 0.31$~$Z_\odot$, and $Z_{z\sim1.6}=0.95 \pm 0.39$~$Z_\odot$}. This is evidence that massive galaxies had enriched rapidly to approximately Solar metallicities as early as $z\sim3$.

It is well known that environment affects galaxy evolution, which is broadly related to supermassive black hole (SMBH) growth. We investigate whether SMBH evolution also depends on host-galaxy local (sub-Mpc) and global ($\approx 1-10$ Mpc) environment. We construct the surface-density field (local environment) and cosmic web (global environment) in the COSMOS field at $z=0.3-3.0$. The environments in COSMOS range from the field to clusters ($M_\mathrm{halo} \lesssim 10^{14}\ M_\odot$), covering the environments where ${\approx 99\%}$ of galaxies in the Universe reside. We measure sample-averaged SMBH accretion rate ($\overline{\mathrm{BHAR}}$) from X-ray observations, and study its dependence on overdensity and cosmic-web environment at different redshifts while controlling for galaxy stellar mass ($M_\star$). Our results show that $\overline{\mathrm{BHAR}}$ does not significantly depend on overdensity or cosmic-web environment once $M_\star$ is controlled, indicating that environment-related physical mechanisms (e.g. tidal interaction and ram-pressure stripping) might not significantly affect SMBH growth. We find that $\overline{\mathrm{BHAR}}$ is strongly related to host-galaxy $M_\star$, regardless of environment.

We present the rest-frame UV and optical photometry and morphology of low-redshift broad-line quasar host galaxies from the Sloan Digital Sky Survey Reverberation Mapping project. Our sample consists of 103 quasars at $z<0.8$, spanning a luminosity range of $-25\le M_g\le -17$ mag. We stack the multi-epoch images in the $g$ and $i$ bands taken by the Canada-France-Hawaii Telescope. The combined $g$-band ($i$-band) images reach a $5\sigma$ depth of 26.2 (25.2) mag, with a typical PSF size of $0.7"$ ($0.6"$). Each quasar is decomposed into a PSF and a S\'ersic profile, representing the central AGN and the host galaxy components, respectively. The systematic errors of the measured host galaxy flux in the two bands are 0.23 and 0.18 mag. The relative errors of the measured galaxy half-light radii ($R_e$) are about 13%. We estimate the rest-frame $u$- and $g$-band flux of the host galaxies, and find that the AGN-to-galaxy flux ratios in the $g$ band are between 0.9 to 4.4 (68.3% confidence). These galaxies have high stellar masses $M_\ast = 10^{10}\sim10^{11}\, M_\odot$. They have similar color with star-forming galaxies at similar redshifts, in consistent with AGN positive feedback in these quasars. We find that the $M_*-M_\text{BH}$ relation in our sample is shallower than the local $M_\text{Bulge}-M_\text{BH}$ relation. The S\'ersic indices and the $M_*-R_e$ relation indicate that the majority of the host galaxies are disk-like.

We present an X-ray point-source catalog from the XMM-Large Scale Structure survey region (XMM-LSS), one of the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS) fields. We target the XMM-LSS region with $1.3$ Ms of new XMM-Newton AO-15 observations, transforming the archival X-ray coverage in this region into a 5.3 deg$^2$ contiguous field with uniform X-ray coverage totaling $2.7$ Ms of flare-filtered exposure, with a $46$ ks median PN exposure time. We provide an X-ray catalog of 5242 sources detected in the soft (0.5-2 keV), hard (2-10 keV), and/or full (0.5-10 keV) bands with a 1% expected spurious fraction determined from simulations. A total of 2381 new X-ray sources are detected compared to previous source catalogs in the same area. Our survey has flux limits of $1.7\times10^{-15}$, $1.3\times10^{-14}$, and $6.5\times10^{-15}$ erg cm$^{-2}$ s$^{-1}$ over 90% of its area in the soft, hard, and full bands, respectively, which is comparable to those of the XMM-COSMOS survey. We identify multiwavelength counterpart candidates for 99.9% of the X-ray sources, of which 93% are considered as reliable based on their matching likelihood ratios. The reliabilities of these high-likelihood-ratio counterparts are further confirmed to be $\approx 97$% reliable based on deep Chandra coverage over $\approx 5$% of the XMM-LSS region. Results of multiwavelength identifications are also included in the source catalog, along with basic optical-to-infrared photometry and spectroscopic redshifts from publicly available surveys. We compute photometric redshifts for X-ray sources in 4.5 deg$^2$ of our field where forced-aperture multi-band photometry is available; $>70$% of the X-ray sources in this subfield have either spectroscopic or high-quality photometric redshifts.

Broad absorption lines (BALs) in the spectra of quasi-stellar objects (QSOs) originate from outflowing winds along our line of sight; winds are thought to originate from the inner regions of the QSO accretion disk, close to the central supermassive black hole (SMBH). Winds likely play a role in galaxy evolution and aid the accretion mechanism onto the SMBH. BAL equivalent widths can change on typical timescales from months to years; such variability is generally attributed to changes in the covering factor and/or in the ionization level of the gas. We investigate BAL variability, focusing on BAL disappearance. We analyze multi-epoch spectra of more than 1500 QSOs -the largest sample ever used for such a study- observed by different programs from the Sloan Digital Sky Survey-I/II/III (SDSS), and search for disappearing C IV BALs. The spectra rest-frame time baseline ranges from 0.28 to 4.9 yr; the source redshifts range from 1.68 to 4.27. We detect 73 disappearing BALs in the spectra of 67 sources. This corresponds to 3.9% of disappearing BALs, and 5.1% of our BAL QSOs exhibit at least one disappearing BAL. We estimate the average lifetime of a BAL along our line of sight (~ 80-100 yr), which appears consistent with the accretion disk orbital time at distances where winds are thought to originate. We inspect properties of the disappearing BALs and compare them to the properties of our main sample. We also investigate the existence of a correlation in the variability of multiple troughs in the same spectrum, and find it persistent at large velocity offsets between BAL pairs, suggesting that a mechanism extending on a global scale is necessary to explain the phenomenon. We select a more reliable sample of disappearing BALs following Filiz Ak et al. (2012), where a subset of our sample was analyzed, and compare the findings from the two works, obtaining generally consistent results.

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.

We present the discovery of KELT-21b, a hot Jupiter transiting the $V=10.5$ A8V star HD 332124. The planet has an orbital period of $P=3.6127647\pm0.0000033$ days and a radius of $1.586_{-0.040}^{+0.039}$ $R_J$. We set an upper limit on the planetary mass of $M_P<3.91$ $M_J$ at $3\sigma$ confidence. We confirmed the planetary nature of the transiting companion using this mass limit and Doppler tomographic observations to verify that the companion transits HD 332124. These data also demonstrate that the planetary orbit is well-aligned with the stellar spin, with a sky-projected spin-orbit misalignment of $\lambda=-5.6_{-1.9}^{+1.7 \circ}$. The star has $T_{\mathrm{eff}}=7598_{-84}^{+81}$ K, $M_*=1.458_{-0.028}^{+0.029}$ $M_{\odot}$, $R_*=1.638\pm0.034$ $R_{\odot}$, and $v\sin I_*=146$ km s$^{-1}$, the highest projected rotation velocity of any star known to host a transiting hot Jupiter. The star also appears to be somewhat metal-poor and $\alpha$-enhanced, with [Fe/H]$=-0.405_{-0.033}^{+0.032}$ and [$\alpha$/Fe]$=0.145 \pm 0.053$; these abundances are unusual, but not extraordinary, for a young star with thin-disk kinematics like KELT-21. High-resolution imaging observations revealed the presence of a pair of stellar companions to KELT-21, located at a separation of 1.2" and with a combined contrast of $\Delta K_S=6.39 \pm 0.06$ with respect to the primary. Although these companions are most likely physically associated with KELT-21, we cannot confirm this with our current data. If associated, the candidate companions KELT-21 B and C would each have masses of $\sim0.12$ $M_{\odot}$, a projected mutual separation of $\sim20$ AU, and a projected separation of $\sim500$ AU from KELT-21. KELT-21b may be one of only a handful of known transiting planets in hierarchical triple stellar systems.

We use the multi-epoch spectra of 362 quasars from the Sloan Digital Sky Survey Reverberation Mapping project to investigate the dependence of the blueshift of CIV relative to MgII on quasar properties. We confirm that high-blueshift sources tend to have low CIV equivalent widths (EWs), and that the low-EW sources span a range of blueshift. Other high-ionization lines, such as HeII, also show similar blueshift properties. The ratio of the line width (measured as both the full-width at half maximum and the velocity dispersion) of CIV to that of MgII increases with blueshift. Quasar variability might enhance the connection between the CIV blueshift and quasar properties (e.g., EW). The variability of the MgII line center (i.e., the wavelength that bisects the cumulative line flux) increases with blueshift. In contrast, the CIV line center shows weaker variability at the extreme blueshifts. Quasars with the high-blueshift CIV lines tend to have less variable continuum emission, when controlling for EW, luminosity, and redshift. Our results support the scenario that high-blueshift sources tend to have large Eddington ratios.

While major mergers have long been proposed as a driver of both AGN activity and the M-sigma relation, studies of moderate to high redshift Seyfert-luminosity AGN hosts have found little evidence for enhanced rates of interactions. However, both theory and observation suggest that while these AGN may be fueled by stochastic accretion and secular processes, high-luminosity, high-redshift, and heavily obscured AGN are the AGN most likely to be merger-driven. To better sample this population of AGN, we turn to infrared selection in the CANDELS/COSMOS field. Compared to their lower-luminosity and less obscured X-ray-only counterparts, IR-only AGN (luminous, heavily obscured AGN) are more likely to be classified as either irregular (50$^{+12}_{-12}$% vs. 9$^{+5}_{-2}$%) or asymmetric (69$^{+9}_{-13}$% vs. 17$^{+6}_{-4}$%) and are less likely to have a spheroidal component (31$^{+13}_{-9}$% vs. 77$^{+4}_{-6}$%). Furthermore, IR-only AGN are also significantly more likely than X-ray-only AGN (75$^{+8}_{-13}$% vs. 31$^{+6}_{-6}$%) to be classified either as interacting or merging in a way that significantly disturbs the host galaxy or disturbed though not clearly interacting or merging, which potentially represents the late stages of a major merger. This suggests that while major mergers may not contribute significantly to the fueling of Seyfert luminosity AGN, interactions appear to play a more dominant role in the triggering and fueling of high-luminosity heavily obscured AGN.

We present composite broad-line region (BLR) reverberation-mapping lag measurements for \halpha, \hbeta, \HeII\,$\lambda4686$ and \MgII\ for a sample of 144, $z\lesssim 1$ quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Using only the 32-epoch spectroscopic light curves in the first 6-month season of SDSS-RM observations, we compile correlation-function measurements for individual objects and then coadd them to allow the measurement of the average lags for our sample at mean redshifts of $0.4$ (for \halpha) and $\sim 0.65$ (for the other lines). At similar quasar luminosities and redshifts, the sample-averaged lag decreases in the order of \MgII, \halpha, \hbeta\ and \HeII. This decrease in lags is accompanied by an increase in the mean line width of the four lines, and is roughly consistent with the virialized motion for BLR gas in photoionization equilibrium. These are among the first RM measurements of stratified BLR structure at $z>0.3$. Dividing our sample by luminosity, \halpha\ shows clear evidence of increasing lags with luminosity, consistent with the expectation from the measured BLR size-luminosity relation based on \hbeta. The other three lines do not show a clear luminosity trend in their average lags due to the limited dynamic range of luminosity probed and the poor average correlation signals in the divided samples, a situation that will be improved with the incorporation of additional photometric and spectroscopic data from SDSS-RM. We discuss the utility and caveats of composite-lag measurements for large statistical quasar samples with reverberation-mapping data.

Previous studies suggest that the growth of supermassive black holes (SMBHs) may be fundamentally related to host-galaxy stellar mass ($M_\star$). To investigate this SMBH growth-$M_\star$ relation in detail, we calculate long-term SMBH accretion rate as a function of $M_\star$ and redshift [$\overline{\rm BHAR}(M_\star, z)$] over ranges of $\log(M_\star/M_\odot)=\text{9.5--12}$ and $z=\text{0.4--4}$. Our $\overline{\rm BHAR}(M_\star, z)$ is constrained by high-quality survey data (GOODS-South, GOODS-North, and COSMOS), and by the stellar mass function and the X-ray luminosity function. At a given $M_\star$, $\overline{\rm BHAR}$ is higher at high redshift. This redshift dependence is stronger in more massive systems (for $\log(M_\star/M_\odot)\approx 11.5$, $\overline{\rm BHAR}$ is three decades higher at $z=4$ than at $z=0.5$), possibly due to AGN feedback. Our results indicate that the ratio between $\overline{\rm BHAR}$ and average star formation rate ($\overline{\rm SFR}$) rises toward high $M_\star$ at a given redshift. This $\overline{\rm BHAR}/\overline{\rm SFR}$ dependence on $M_\star$ does not support the scenario that SMBH and galaxy growth are in lockstep. We calculate SMBH mass history [$M_{\rm BH}(z)$] based on our $\overline{\rm BHAR}(M_\star, z)$ and the $M_\star(z)$ from the literature, and find that the $M_{\rm BH}$-$M_\star$ relation has weak redshift evolution since $z\approx 2$. The $M_{\rm BH}/M_\star$ ratio is higher toward massive galaxies: it rises from $\approx 1/5000$ at $\log M_\star\lesssim 10.5$ to $\approx 1/500$ at $\log M_\star \gtrsim 11.2$. Our predicted $M_{\rm BH}/M_\star$ ratio at high $M_\star$ is similar to that observed in local giant ellipticals, suggesting that SMBH growth from mergers is unlikely to dominate over growth from accretion.

We examine the fraction of massive ($M_{*}>10^{10} M_{\odot}$), compact star-forming galaxies (cSFGs) that host an active galactic nucleus (AGN) at $z\sim2$. These cSFGs are likely the direct progenitors of the compact quiescent galaxies observed at this epoch, which are the first population of passive galaxies to appear in large numbers in the early Universe. We identify cSFGs that host an AGN using a combination of Hubble WFC3 imaging and Chandra X-ray observations in four fields: the Chandra Deep Fields, the Extended Groth Strip, and the UKIDSS Ultra Deep Survey field. We find that $39.2^{+3.9}_{-3.6}$\% (65/166) of cSFGs at $1.4<z<3.0$ host an X-ray detected AGN. This fraction is 3.2 times higher than the incidence of AGN in extended star-forming galaxies with similar masses at these redshifts. This difference is significant at the $6.2\sigma$ level. Our results are consistent with models in which cSFGs are formed through a dissipative contraction that triggers a compact starburst and concurrent growth of the central black hole. We also discuss our findings in the context of cosmological galaxy evolution simulations that require feedback energy to rapidly quench cSFGs. We show that the AGN fraction peaks precisely where energy injection is needed to reproduce the decline in the number density of cSFGs with redshift. Our results suggest that the first abundant population of massive, quenched galaxies emerged directly following a phase of elevated supermassive black hole growth and further hints at a possible connection between AGN and the rapid quenching of star formation in these galaxies.

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.

We investigate the dependence of black-hole accretion rate (BHAR) on host-galaxy star formation rate (SFR) and stellar mass ($M_*$) in the CANDELS/GOODS-South field in the redshift range of $0.5\leq z < 2.0$. Our sample consists of $\approx 18000$ galaxies, allowing us to probe galaxies with $0.1 \lesssim \mathrm{SFR} \lesssim 100\ M_\odot\ \mathrm{yr}^{-1}$ and/or $10^8 \lesssim M_* \lesssim 10^{11}\ M_{\odot}$. We use sample-mean BHAR to approximate long-term average BHAR. Our sample-mean BHARs are derived from the $Chandra$ Deep Field-South 7 Ms observations, while the SFRs and $M_*$ have been estimated by the CANDELS team through SED fitting. The average BHAR is correlated positively with both SFR and $M_*$, and the BHAR-SFR and BHAR-$M_*$ relations can both be described acceptably by linear models with a slope of unity. However, BHAR appears to be correlated more strongly with $M_*$ than SFR. This result indicates that $M_*$ is the primary host-galaxy property related to black-hole growth, and the apparent BHAR-SFR relation is largely a secondary effect due to the star-forming main sequence. Among our sources, massive galaxies ($M_* \gtrsim 10^{10} M_{\odot}$) have significantly higher BHAR/SFR ratios than less-massive galaxies, indicating the former have higher black-hole fueling efficiency and/or higher SMBH occupation fraction than the latter. Our results can naturally explain the observed proportionality between $M_{\rm BH}$ and $M_*$ for local giant ellipticals, and suggest their $M_{\rm BH}/M_*$ is higher than that of local star-forming galaxies. Among local star-forming galaxies, massive systems might have higher $M_{\rm BH}/M_*$ compared to dwarfs.

We explore the evolution of the internal gas kinematics of star-forming galaxies from the peak of cosmic star-formation at $z\sim2$ to today. Measurements of galaxy rotation velocity $V_{rot}$, which quantify ordered motions, and gas velocity dispersion $\sigma_g$, which quantify disordered motions, are adopted from the DEEP2 and SIGMA surveys. This sample covers a continuous baseline in redshift from $z=2.5$ to $z=0.1$, spanning 10 Gyrs. At low redshift, nearly all sufficiently massive star-forming galaxies are rotationally supported ($V_{rot}>\sigma_g$). By $z=2$, the percentage of galaxies with rotational support has declined to 50$\%$ at low stellar mass ($10^{9}-10^{10}\,M_{\odot}$) and 70$\%$ at high stellar mass ($10^{10}-10^{11}M_{\odot}$). For $V_{rot}\,>\,3\,\sigma_g$, the percentage drops below 35$\%$ for all masses. From $z\,=\,2$ to now, galaxies exhibit remarkably smooth kinematic evolution on average. All galaxies tend towards rotational support with time, and it is reached earlier in higher mass systems. This is mostly due to an average decline in $\sigma_g$ by a factor of 3 since a redshift of 2, which is independent of mass. Over the same time period, $V_{rot}$ increases by a factor of 1.5 for low mass systems, but does not evolve for high mass systems. These trends in $V_{rot}$ and $\sigma_g$ with time are at a fixed stellar mass and should not be interpreted as evolutionary tracks for galaxy populations. When galaxy populations are linked in time with abundance matching, not only does $\sigma_g$ decline with time as before, but $V_{rot}$ strongly increases with time for all galaxy masses. This enhances the evolution in $V_{rot}/\sigma_g$. These results indicate that $z\,=\,2$ is a period of disk assembly, during which the strong rotational support present in today's massive disk galaxies is only just beginning to emerge.

Very few low-ionization broad absorption line (LoBAL) QSOs have been found at high redshifts to date. One high-redshift LoBAL QSO, J0122+1216, was recently discovered at the Lijiang 2.4-m Telescope with an initial redshift determination of 4.76. Aiming to investigate its physical properties, we carried out follow-up observations in the optical and near-IR spectroscopy. Near-IR spectra from UKIRT and P200 confirms that it is a LoBAL, with a new redshift determination of $4.82\pm0.01$ based on the \mgii~ emission-line. The new \mgii~ redshift determination reveals strong blueshifts and asymmetry of the high-ionization emission lines. We estimated a black hole mass of $\sim 2.3\times 10^9 M_\odot$ and Eddington ratio of $\sim 1.0$ according to the empirical \mgii-based single-epoch relation and bolometric correction factor. It is possible that strong outflows are the result of an extreme quasar environment driven by the high Eddington ratio. A lower limit on the outflowing kinetic power ($>0.9\% L_{Edd}$) was derived from both emission and absorption lines, indicating these outflows play a significant role in the feedback process to regulate the growth of its black hole as well as host galaxy evolution.

We present a sample of 10 low-mass active galactic nuclei (AGNs) selected from the 40-month NuSTAR serendipitous survey. The sample is selected to have robust NuSTAR detections at $3 - 24$~keV, to be at $z < 0.3$, and to have optical r-band magnitudes at least 0.5~mag fainter than an $L_\star$ galaxy at its redshift. The median values of absolute magnitude, stellar mass and 2--10 X-ray luminosity of our sample are $\langle M_r\rangle = -20.03$, $\langle M_\star\rangle = 4.6\times10^{9}M_\odot$, and $\langle L_{2-10\mathrm{keV}}\rangle = 3.1\times10^{42}$ erg s$^{-1}$, respectively. Five objects have detectable broad H$\alpha$ emission in their optical spectra, indicating black-hole masses of $(1.1-10.4)\times 10^6 M_\odot$. We find that $30^{+17}_{-10}\%$ of the galaxies in our sample do not show AGN-like optical narrow emission lines, and one of the ten galaxies in our sample, J115851+4243.2, shows evidence for heavy X-ray absorption. This result implies that a non-negligible fraction of low-mass galaxies might harbor accreting massive black holes that are missed by optical spectroscopic surveys and $<10$ keV X-ray surveys. The mid-IR colors of our sample also indicate these optically normal low-mass AGNs cannot be efficiently identified with typical AGN selection criteria based on WISE colors. While the hard ($>10$ keV) X-ray selected low-mass AGN sample size is still limited, our results show that sensitive NuSTAR observations are capable of probing faint hard X-ray emission originating from the nuclei of low-mass galaxies out to moderate redshift ($z<0.3$), thus providing a critical step in understanding AGN demographics in low-mass galaxies.

Quasar emission lines are often shifted from the systemic velocity due to various dynamical and radiative processes in the line-emitting region. The level of these velocity shifts depends both on the line species and on quasar properties. We study velocity shifts for the line peaks of various narrow and broad quasar emission lines relative to systemic using a sample of 849 quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. The coadded (from 32 epochs) spectra of individual quasars have sufficient signal-to-noise ratio (SNR) to measure stellar absorption lines to provide reliable systemic velocity estimates, as well as weak narrow emission lines. The sample also covers a large dynamic range in quasar luminosity (~2 dex), allowing us to explore potential luminosity dependence of the velocity shifts. We derive average line peak velocity shifts as a function of quasar luminosity for different lines, and quantify their intrinsic scatter. We further quantify how well the peak velocity can be measured for various lines as a function of continuum SNR, and demonstrate there is no systematic bias in the line peak measurements when the spectral quality is degraded to as low as SNR~3 per SDSS pixel. Based on the observed line shifts, we provide empirical guidelines on redshift estimation from [OII]3728, [OIII]5008, [NeV]3426, MgII, CIII], HeII1640, broad Hbeta, CIV, and SiIV, which are calibrated to provide unbiased systemic redshifts in the mean, but with increasing intrinsic uncertainties of 46, 56, 119, 205, 233, 242, 400, 415, and 477 km/s, in addition to the measurement uncertainties. These more realistic redshift uncertainties are generally much larger than the formal uncertainties reported by the redshift pipelines for spectroscopic quasar surveys, and demonstrate the infeasibility of measuring quasar redshifts to better than ~200 km/s with only broad lines.

Galaxy interactions are thought to be one of the main triggers of Active Galactic Nuclei (AGN), especially at high luminosities, where the accreted gas mass during the AGN lifetime is substantial. Evidence for a connection between mergers and AGN, however, remains mixed. Possible triggering mechanisms remain particularly poorly understood for luminous AGN, which are thought to require triggering by major mergers, rather than secular processes. We analyse the host galaxies of a sample of 20 optically and X-ray selected luminous AGN (log($L_{bol}$ [erg/s]) $>$ 45) at z $\sim$ 0.6 using HST WFC3 data in the F160W/H band. 15/20 sources have resolved host galaxies. We create a control sample of mock AGN by matching the AGN host galaxies to a control sample of non-AGN galaxies. Visual signs of disturbances are found in about 25% of sources in both the AGN hosts and control galaxies. Using both visual classification and quantitative morphology measures, we show that the levels of disturbance are not enhanced when compared to a matched control sample. We find no signs that major mergers play a dominant role in triggering AGN at high luminosities, suggesting that minor mergers and secular processes dominate AGN triggering up to the highest AGN luminosities. The upper limit on the enhanced fraction of major mergers is $\leqslant$20%. While major mergers might increase the incidence of (luminous AGN), they are not the prevalent triggering mechanism in the population of unobscured AGN.

We investigate the burstiness of star formation histories (SFHs) of galaxies at $0.4<z<1$ by using the ratio of star formation rates (SFRs) measured from H$\beta$ and FUV (1500 \AA) (H$\beta$--to--FUV ratio). Our sample contains 164 galaxies down to stellar mass (M*) of $10^{8.5} M_\odot$ in the CANDELS GOODS-N region, where Team Keck Redshift Survey DEIMOS spectroscopy and HST/WFC3 F275W images from CANDELS and Hubble Deep UV Legacy Survey are available. When the {\it ratio} of H$\beta$- and FUV-derived SFRs is measured, dust extinction correction is negligible (except for very dusty galaxies) with the Calzetti attenuation curve. The H$\beta$--to--FUV ratio of our sample increases with M* and SFR. The median ratio is $\sim$0.7 at M*$\sim10^{8.5} M_\odot$ (or SFR$\sim 0.5 M_\odot/yr$) and increases to $\sim$1 at M*$\sim10^{10} M_\odot$ (or SFR $\sim 10 M_\odot/yr$). At M*$<10^{9.5} M_\odot$, our median H$\beta$--to--FUV ratio is lower than that of local galaxies at the same M*, implying a redshift evolution. Bursty SFH on a timescale of a few tens of megayears on galactic scales provides a plausible explanation of our results, and the importance of the burstiness increases as M* decreases. Due to sample selection effects, our H$\beta$--to--FUV ratio may be an upper limit of the true value of a complete sample, which strengthens our conclusions. Other models, e.g., non-universal initial mass function or stochastic star formation on star cluster scales, are unable to plausibly explain our results.

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.

We exploit the 7 Ms \textit{Chandra} observations in the \chandra\,Deep Field-South (\mbox{CDF-S}), the deepest X-ray survey to date, coupled with CANDELS/GOODS-S data, to measure the total X-ray emission arising from 2076 galaxies at $3.5\leq z < 6.5$. This aim is achieved by stacking the \textit{Chandra} data at the positions of optically selected galaxies, reaching effective exposure times of $\geq10^9\mathrm{s}$. We detect significant ($>3.7\sigma$) X-ray emission from massive galaxies at $z\approx4$. We also report the detection of massive galaxies at $z\approx5$ at a $99.7\%$ confidence level ($2.7\sigma$), the highest significance ever obtained for X-ray emission from galaxies at such high redshifts. No significant signal is detected from galaxies at even higher redshifts. The stacking results place constraints on the BHAD associated with the known high-redshift galaxy samples, as well as on the SFRD at high redshift, assuming a range of prescriptions for X-ray emission due to X- ray binaries. We find that the X-ray emission from our sample is likely dominated by processes related to star formation. Our results show that low-rate mass accretion onto SMBHs in individually X-ray-undetected galaxies is negligible, compared with the BHAD measured for samples of X-ray detected AGN, for cosmic SMBH mass assembly at high redshift. We also place, for the first time, constraints on the faint-end of the AGN X-ray luminosity function ($\mathrm{logL_X\sim42}$) at $z>4$, with evidence for fairly flat slopes. The implications of all of these findings are discussed in the context of the evolution of the AGN population at high redshift.

We present spatially-resolved Atacama Large Millimeter/sub-millimeter Array (ALMA) 870 $\mu$m dust continuum maps of six massive, compact, dusty star-forming galaxies (SFGs) at $z\sim2.5$. These galaxies are selected for their small rest-frame optical sizes ($r_{\rm e, F160W}\sim1.6$ kpc) and high stellar-mass densities that suggest that they are direct progenitors of compact quiescent galaxies at $z\sim2$. The deep observations yield high far-infrared (FIR) luminosities of L$_{\rm IR}=10^{12.3-12.8}$ L$_{\odot}$ and star formation rates (SFRs) of SFR$=200-700$ M$_{\odot}$yr$^{-1}$, consistent with those of typical star-forming "main sequence" galaxies. The high-spatial resolution (FWHM$\sim$0.12"-0.18") ALMA and HST photometry are combined to construct deconvolved, mean radial profiles of their stellar mass and (UV+IR) SFR. We find that the dusty, nuclear IR-SFR overwhelmingly dominates the bolometric SFR up to $r\sim5$ kpc, by a factor of over 100$\times$ from the unobscured UV-SFR. Furthermore, the effective radius of the mean SFR profile ($r_{\rm e, SFR}\sim1$ kpc) is $\sim$30% smaller than that of the stellar mass profile. The implied structural evolution, if such nuclear starburst last for the estimated gas depletion time of $\Delta t=\pm100$ Myr, is a 4$\times$ increase of the stellar mass density within the central 1 kpc and a 1.6$\times$ decrease of the half-mass radius. This structural evolution fully supports dissipation-driven, formation scenarios in which strong nuclear starbursts transform larger, star-forming progenitors into compact quiescent galaxies.

We present results from a survey of the internal kinematics of 49 star-forming galaxies at z$\,\sim\,$2 in the CANDELS fields with the Keck/MOSFIRE spectrograph (SIGMA, Survey in the near-Infrared of Galaxies with Multiple position Angles). Kinematics (rotation velocity $V_{rot}$ and integrated gas velocity dispersion $\sigma_g$) are measured from nebular emission lines which trace the hot ionized gas surrounding star-forming regions. We find that by z$\,\sim\,$2, massive star-forming galaxies ($\log\,M_*/M_{\odot}\gtrsim10.2$) have assembled primitive disks: their kinematics are dominated by rotation, they are consistent with a marginally stable disk model, and they form a Tully-Fisher relation. These massive galaxies have values of $V_{rot}/\sigma_g$ which are factors of 2-5 lower than local well-ordered galaxies at similar masses. Such results are consistent with findings by other studies. We find that low mass galaxies ($\log\,M_*/M_{\odot}\lesssim10.2$) at this epoch are still in the early stages of disk assembly: their kinematics are often supported by gas velocity dispersion and they fall from the Tully-Fisher relation to significantly low values of $V_{rot}$. This "kinematic downsizing" implies that the process(es) responsible for disrupting disks at z$\,\sim\,$2 have a stronger effect and/or are more active in low mass systems. In conclusion, we find that the period of rapid stellar mass growth at z$\,\sim\,$2 is coincident with the nascent assembly of low mass disks and the assembly and settling of high mass disks.

We use the coadded spectra of 32 epochs of Sloan Digital Sky Survey (SDSS) Reverberation Mapping Project observations of 482 quasars with z>1.46 to highlight systematic biases in the SDSS- and BOSS-pipeline redshifts due to the natural diversity of quasar properties. We investigate the characteristics of this bias by comparing the BOSS-pipeline redshifts to an estimate from the centroid of HeII 1640. HeII has a low equivalent width but is often well-defined in high-S/N spectra, does not suffer from self-absorption, and has a narrow component that, when present (the case for about half of our sources), produces a redshift estimate that, on average, is consistent with that determined from [OII] to within 1-sigma of the quadrature sum of the HeII and [OII] centroid measurement uncertainties. The large redshift differences of ~1000 km/s, on average, between the BOSS-pipeline and HeII-centroid redshifts suggest there are significant biases in a portion of BOSS quasar redshift measurements. Adopting the HeII-based redshifts shows that CIV does not exhibit a ubiquitous blueshift for all quasars, given the precision probed by our measurements. Instead, we find a distribution of CIV centroid blueshifts across our sample, with a dynamic range that (i) is wider than that previously reported for this line, and (ii) spans CIV centroids from those consistent with the systemic redshift to those with significant blueshifts of thousands of kilometers per second. These results have significant implications for measurement and use of high-redshift quasar properties and redshifts and studies based thereon.

Differentiating between active galactic nuclei (AGN) activity and star formation in z ~ 2 galaxies is difficult because traditional methods, such as line ratio diagnostics, change with redshift while multi-wavelength methods (X-ray, radio, IR) are sensitive to only the brightest AGN. We have developed a new method for spatially resolving emission lines in HST/WFC3 G141 grism spectra and quantifying AGN activity through the spatial gradient of the [O III]/H$\beta$ line ratio. Through detailed simulations, we show that our novel line-ratio gradient approach identifies ~ sim 40% more low-mass and obscured AGN than obtained by classical methods. Based on our simulations, we developed a relationship that maps stellar mass, star formation rate, and measured [O III]/H$\beta$ gradient to AGN Eddington ratio. We apply our technique to previously studied stacked samples of galaxies at z ~2 and find that our results are consistent with these studies. Using this gradient method will also be able to inform other galaxy evolution science, such as inside-out quenching and metallicity gradients, and will be widely applicable to future spatially resolved JWST data.

We present results from the largest systematic investigation of broad absorption line (BAL) acceleration to date. We use spectra of 140 quasars from three Sloan Digital Sky Survey programs to search for global velocity offsets in BALs over timescales of ~2.5-5.5 years in the quasar rest frame. We carefully select acceleration candidates by requiring monolithic velocity shifts over the entire BAL trough, avoiding BALs with velocity shifts that might be caused by profile variability. The CIV BALs of two quasars show velocity shifts consistent with the expected signatures of BAL acceleration, and the BAL of one quasar shows a velocity-shift signature of deceleration. In our two acceleration candidates, we see evidence that the magnitude of the acceleration is not constant over time; the magnitudes of the change in acceleration for both acceleration candidates are difficult to produce with a standard disk-wind model or via geometric projection effects. We measure upper limits to acceleration and deceleration for 76 additional BAL troughs and find that the majority of BALs are stable to within about 3% of their mean velocities. The lack of widespread acceleration/deceleration could indicate that the gas producing most BALs is located at large radii from the central black hole and/or is not currently strongly interacting with ambient material within the host galaxy along our line of sight.

We present measurements of the evolution of normal-galaxy X-ray emission from $z \approx$ 0-7 using local galaxies and galaxy samples in the 6 Ms Chandra Deep Field-South (CDF-S) survey. The majority of the CDF-S galaxies are observed at rest-frame energies above 2 keV, where the emission is expected to be dominated by X-ray binary (XRB) populations; however, hot gas is expected to provide small contributions to the observed- frame < 1 keV emission at $z < 1$. We show that a single scaling relation between X-ray luminosity ($L_{\rm X}$) and star-formation rate (SFR) is insufficient for characterizing the average X-ray emission at all redshifts. We establish that scaling relations involving not only SFR, but also stellar mass ($M_\star$) and redshift, provide significantly improved characterizations of the average X-ray emission from normal galaxy populations at $z \approx$ 0-7. We further provide the first empirical constraints on the redshift evolution of X-ray emission from both low-mass XRB (LMXB) and high-mass XRB (HMXB) populations and their scalings with $M_\star$ and SFR, respectively. We find $L_{\rm 2-10 keV}({\rm LMXB})/M_\star \propto (1 + z)^{2-3}$ and $L_{\rm 2-10 keV}({\rm HMXB})$/SFR $\propto (1 + z)$, and show that these relations are consistent with XRB population-synthesis model predictions, which attribute the increase in LMXB and HMXB scaling relations with redshift as being due to declining host galaxy stellar ages and metallicities, respectively. We discuss how emission from XRBs could provide an important source of heating to the intergalactic medium in the early Universe, exceeding that of active galactic nuclei.

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 investigate the dependence on data quality of quasar properties measured from the CIV emission line region at high redshifts. Our measurements come from 32 epochs of Sloan Digital Sky Survey (SDSS) Reverberation Mapping Project spectroscopic observations of 482 z>1.46 quasars. We compare the differences between measurements made from the single-epoch and coadded spectra, focusing on the CIV emission line because of its importance for studies of high-redshift quasar demographics and physical properties, including black hole masses. In addition to statistical errors increasing (by factors of ~2-4), we find increasing systematic offsets with decreasing S/N. The systematic difference (measurement uncertainty) in our lowest S/N (<5) subsample between the single-epoch and coadded spectrum (i) CIV equivalent width is 17A (31A), (ii) centroid wavelength is <1A (2A), and fractional velocity widths, Delta V/V, characterized by (iii) the line dispersion, sigma_l, is 0.104 (0.12), and (iv) the mean absolute deviation (MAD) is 0.072 (0.11). These remain smaller than the 1-sigma measurement uncertainties for all subsamples considered. The MAD is found to be the most robust line-width characterization. Offsets in the CIV full-width at half maximum (FWHM) velocity width and the CIV profile characterized by FWHM/sigma_l are only smaller than the statistical uncertainties when S/N>10, although offsets in lower S/N spectra exceed the statistical uncertainties by only a factor of ~1.5. Characterizing the CIV line profile by the kurtosis is the least robust property investigated, as the median systematic coadded--single-epoch measurement differences are larger than the statistical uncertainties for all S/N subsamples.

We study a set of 3319 galaxies in the redshift interval 0.04 < z < 0.15 with far-infrared (FIR) coverage from the Herschel Stripe 82 survey (HerS), and emission-line measurements, redshifts, stellar masses and star-formation rates (SFRs) from the SDSS (DR7) MPA/JHU database. About 40% of the sample are detected in the Herschel/SPIRE 250 micron band. Total infrared (TIR) luminosities derived from HerS and ALLWISE photometry allow us to compare infrared and optical estimates of SFR with unprecedented statistics for diverse classes of galaxies. We find excellent agreement between TIR-derived and emission line-based SFRs for H II galaxies. Other classes, such as active galaxies and evolved galaxies, exhibit systematic discrepancies between optical and TIR SFRs. We demonstrate that these offsets are attributable primarily to survey biases and the large intrinsic uncertainties of the D4000- and colour-based optical calibrations used to estimate the SDSS SFRs of these galaxies. Using a classification scheme which expands upon popular emission-line methods, we demonstrate that emission-line galaxies with uncertain classifications include a population of massive, dusty, metal-rich star-forming systems that are frequently neglected in existing studies. We also study the capabilities of infrared selection of star-forming galaxies. FIR selection reveals a substantial population of galaxies dominated by cold dust which are missed by the long-wavelength WISE bands. Our results demonstrate that Herschel large-area surveys offer the means to construct large, relatively complete samples of local star-forming galaxies with accurate estimates of SFR that can be used to study the interplay between nuclear activity and star-formation.

Reverberation mapping (RM) measurements of broad-line region (BLR) lags in z>0.3 quasars are important for directly measuring black hole masses in these distant objects, but so far there have been limited attempts and success given the practical difficulties of RM in this regime. Here we report preliminary results of 15 BLR lag measurements from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, a dedicated RM program with multi-object spectroscopy designed for RM over a wide redshift range. The lags are based on the 2014 spectroscopic light curves alone (32 epochs over 6 months) and focus on the Hbeta and MgII broad lines in the 100 lowest-redshift (z<0.8) quasars included in SDSS-RM; they represent a small subset of the lags that SDSS-RM (including 849 quasars to z~4.5) is expected to deliver. The reported preliminary lag measurements are for intermediate-luminosity quasars at 0.3<~z<0.8, including 9 Hbeta lags and 6 MgII lags, for the first time extending RM results to this redshift-luminosity regime and providing direct quasar black hole mass estimates over ~ half of cosmic time. The MgII lags also increase the number of known MgII lags by several-fold, and start to explore the utility of MgII for RM at high redshift. The location of these new lags at higher redshifts on the observed BLR size-luminosity relationship is statistically consistent with previous Hbeta results at z<0.3. However, an independent constraint on the relationship slope at z>0.3 is not yet possible due to the limitations in our current sample. Our results demonstrate the general feasibility and potential of multi-object RM for z>0.3 quasars.

We compare the physical and morphological properties of z ~ 2 Lyman-alpha emitting galaxies (LAEs) identified in the HETDEX Pilot Survey and narrow band studies with those of z ~ 2 optical emission line galaxies (oELGs) identified via HST WFC3 infrared grism spectroscopy. Both sets of galaxies extend over the same range in stellar mass (7.5 < logM < 10.5), size (0.5 < R < 3.0 kpc), and star-formation rate (~1 < SFR < 100). Remarkably, a comparison of the most commonly used physical and morphological parameters -- stellar mass, half-light radius, UV slope, star formation rate, ellipticity, nearest neighbor distance, star formation surface density, specific star formation rate, [O III] luminosity, and [O III] equivalent width -- reveals no statistically significant differences between the populations. This suggests that the processes and conditions which regulate the escape of Ly-alpha from a z ~ 2 star-forming galaxy do not depend on these quantities. In particular, the lack of dependence on the UV slope suggests that Ly-alpha emission is not being significantly modulated by diffuse dust in the interstellar medium. We develop a simple model of Ly-alpha emission that connects LAEs to all high-redshift star forming galaxies where the escape of Ly-alpha depends on the sightline through the galaxy. Using this model, we find that mean solid angle for Ly-alpha escape is 2.4+/-0.8 steradians; this value is consistent with those calculated from other studies.

Astronomers in CANDELS outline changes for the academic system to promote a smooth transition for junior scientists from academia to industry.

We explore the evolution of the Stellar Mass-Star Formation Rate-Metallicity Relation using a set of 256 COSMOS and GOODS galaxies in the redshift range 1.90 < z < 2.35. We present the galaxies' rest-frame optical emission-line fluxes derived from IR-grism spectroscopy with the Hubble Space Telescope and combine these data with star formation rates and stellar masses obtained from deep, multi-wavelength (rest-frame UV to IR) photometry. We then compare these measurements to those for a local sample of galaxies carefully matched in stellar mass (7.5 < log(M*/Msol) < 10.5) and star formation rate (-0.5 < log(SFR) < 2.5 in Msol yr^-1). We find that the distribution of z ~ 2.1 galaxies in stellar mass-SFR-metallicity space is clearly different from that derived for our sample of similarly bright (L_H\b{eta} > 3 . 10^40 ergs s^-1) local galaxies, and this offset cannot be explained by simple systematic offsets in the derived quantities. At stellar masses above ~10^9 Msol and star formation rates above ~10 Msol yr^-1, the z ~ 2.1 galaxies have higher oxygen abundances than their local counterparts, while the opposite is true for lower-mass, lower-SFR systems.

We derive the mean wavelength dependence of stellar attenuation in a sample of 239 high redshift (1.90 < z < 2.35) galaxies selected via Hubble Space Telescope (HST) WFC3 IR grism observations of their rest-frame optical emission lines. Our analysis indicates that the average reddening law follows a form similar to that derived by Calzetti et al. for local starburst galaxies. However, over the mass range 7.2 < log M/Msolar < 10.2, the slope of the attenuation law in the UV is shallower than that seen locally, and the UV slope steepens as the mass increases. These trends are in qualitative agreement with Kriek & Conroy, who found that the wavelength dependence of attenuation varies with galaxy spectral type. However, we find no evidence of an extinction "bump" at 2175 A in any of the three stellar mass bins, or in the sample as a whole. We quantify the relation between the attenuation curve and stellar mass and discuss its implications.

We present the Team Keck Redshift Survey 2 (TKRS2), a near-infrared spectral observing program targeting selected galaxies within the CANDELS subsection of the GOODS-North Field. The TKRS2 program exploits the unique capabilities of MOSFIRE, an infrared multi-object spectrometer which entered service on the Keck I telescope in 2012 and contributes substantially to the study of galaxy spectral features at redshifts inaccessible to optical spectrographs. The TKRS2 project targets 97 galaxies drawn from samples that include z~2 emission-line galaxies with features observable in the JHK bands as well as lower-redshift targets with features in the Y band. We present a detailed measurement of MOSFIRE's sensitivity as a function of wavelength, including the effects of telluric features across the YJHK filters. The largest utility of our survey is in providing rest-frame-optical emission lines for z>1 galaxies, and we demonstrate that the ratios of strong, optical emission lines of z~2 galaxies suggest the presence of either higher N/O abundances than are found in z~0 galaxies or low-metallicity gas ionized by an active galactic nucleus. We have released all TKRS2 data products into the public domain to allow researchers access to representative raw and reduced MOSFIRE spectra.

We examine the host morphologies of heavily obscured active galactic nuclei (AGN) at $z\sim1$ to test whether obscured supermassive black hole growth at this epoch is preferentially linked to galaxy mergers. Our sample consists of 154 obscured AGN with $N_{\rm H}>10^{23.5}$ cm$^{-2}$ and $z<1.5$. Using visual classifications, we compare the morphologies of these AGN to control samples of moderately obscured ($10^{22}$ cm$^{-2}$ $<N_{\rm H}< 10^{23.5}$ cm$^{-2}$) and unobscured ($N_{\rm H}<10^{22}$ cm$^{-2}$) AGN. These control AGN are matched in redshift and intrinsic X-ray luminosity to our heavily obscured AGN. We find that heavily obscured AGN at z~1 are twice as likely to be hosted by late-type galaxies relative to unobscured AGN ($65.3^{+4.1}_{-4.6}\%$ vs $34.5^{+2.9}_{-2.7}\%$) and three times as likely to exhibit merger or interaction signatures ($21.5^{+4.2}_{-3.3}\%$ vs $7.8^{+1.9}_{-1.3}\%$). The increased merger fraction is significant at the 3.8$\sigma$ level. We also find that the incidence of point-like morphologies is inversely proportional to obscuration. If we exclude all point sources and consider only extended hosts, we find the correlation between merger fraction and obscuration is still evident, however at a reduced statistical significance ($2.5\sigma$). The fact that we observe a different disk/spheroid fraction versus obscuration indicates that viewing angle cannot be the only thing differentiating our three AGN samples, as a simple unification model would suggest. The increased fraction of disturbed morphologies with obscuration supports an evolutionary scenario, in which Compton-thick AGN are a distinct phase of obscured SMBH growth following a merger/interaction event. Our findings also suggest that some of the merger-triggered SMBH growth predicted by recent AGN fueling models may be hidden among the heavily obscured, Compton-thick population.

We use 317,000 emission-line galaxies from the Sloan Digital Sky Survey to investigate line-ratio selection of active galactic nuclei (AGNs). In particular, we demonstrate that "star formation dilution" by HII regions causes a significant bias against AGN selection in low-mass, blue, star-forming, disk-dominated galaxies. This bias is responsible for the observed preference of AGNs among high-mass, green, moderately star-forming, bulge-dominated hosts. We account for the bias and simulate the intrinsic population of emission-line AGNs using a physically-motivated Eddington ratio distribution, intrinsic AGN narrow line region line ratios, a luminosity-dependent Lbol/L[OIII] bolometric correction, and the observed Mbh-sigma relation. These simulations indicate that, in massive (log(M*/Msun) > 10) galaxies, AGN accretion is correlated with specific star formation rate but is otherwise uniform with stellar mass. There is some hint of lower black hole occupation in low-mass (log(M*/Msun) < 10) hosts, although our modeling is limited by uncertainties in measuring and interpreting the velocity dispersions of low-mass galaxies. The presence of star formation dilution means that AGNs contribute little to the observed strong optical emission lines (e.g., [OIII] and Ha) in low-mass and star-forming hosts. However the AGN population recovered by our modeling indicates that feedback by typical (low- to moderate-accretion) low-redshift AGNs has nearly uniform efficiency at all stellar masses, star formation rates, and morphologies. Taken together, our characterization of the observational bias and resultant AGN occupation function suggest that AGNs are unlikely to be the dominant source of star formation quenching in galaxies, but instead are fueled by the same gas which drives star formation activity.

We study the evolution of the core (r<1 kpc) and effective (r<r_e) stellar-mass surface densities, in star-forming and quiescent galaxies. Since z=3, both populations occupy distinct, linear relations in log(Sigma_e) and log(Sigma_1) vs. log(M). These structural relations exhibit slopes and scatter that remain almost constant with time while their normalizations decline. For SFGs, the normalization declines by less than a factor of 2 from z=3, in both Sigma_e and Sigma_1. Such mild declines suggest that SFGs build dense cores by growing along these relations. We define this evolution as the structural main sequence (Sigma-MS). Quiescent galaxies follow different relations (Sigma^Q_e, Sigma^Q_1) off the Sigma-MS by having higher densities than SFGs of the same mass and redshift. The normalization of Sigma^Q_e declines by a factor of 10 since z=3, but only a factor of 2 in Sigma^Q_1. Thus, the common denominator for quiescent galaxies at all redshifts is the presence of a dense stellar core, and the formation of such cores in SFGs is the main requirement for quenching. Expressed in 2D as deviations off the SFR-MS and off Sigma^Q_1 at each redshift, the distribution of massive galaxies forms a universal, L-shaped sequence that relates two fundamental physical processes: compaction and quenching. Compaction is a process of substantial core-growth in SFGs relative to that in the Sigma-MS. This process increases the core-to-total mass and Sersic index, thereby, making compact SFGs. Quenching occurs once compact SFGs reach a maximum central density above Sigma^Q_1 > 9.5 M_sun/kpc^2. This threshold provides the most effective selection criterion to identify the star-forming progenitors of quiescent galaxies at all redshifts.

Quasar host galaxies are key for understanding the relation between galaxies and the supermassive black holes (SMBHs) at their centers. We present a study of 191 broad-line quasars and their host galaxies at z < 1, using high signal-to-noise ratio (SNR) spectra produced by the Sloan Digital Sky Survey Reverberation Mapping project. Clear detection of stellar absorption lines allows a reliable decomposition of the observed spectra into nuclear and host components, using spectral models of quasar and stellar radiations as well as emission lines from the interstellar medium. We estimate age, mass (M*), and velocity dispersion (sigma*) of the host stars, the star formation rate (SFR), quasar luminosity, and SMBH mass (Mbh), for each object. The quasars are preferentially hosted by massive galaxies with M* ~ 10^{11} Msun characterized by stellar ages around a billion years, which coincides with the transition phase of normal galaxies from the blue cloud to the red sequence. The host galaxies have relatively low SFRs and fall below the main sequence of star-forming galaxies at similar redshifts. These facts suggest that the hosts have experienced an episode of major star formation sometime in the past billion years, which was subsequently quenched or suppressed. The derived Mbh - sigma* and Mbh - M* relations agree with our past measurements and are consistent with no evolution from the local Universe. The present analysis demonstrates that reliable measurements of stellar properties of quasar host galaxies are possible with high-SNR fiber spectra, which will be acquired in large numbers with future powerful instruments such as the Subaru Prime Focus Spectrograph.

We explore the variability of quasars in the MgII and Hbeta broad emission lines and UV/optical continuum emission using the Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM). This is the largest spectroscopic study of quasar variability to date: our study includes 29 spectroscopic epochs from SDSS-RM over $6$ months, containing 357 quasars with MgII and 41 quasars with Hbeta . On longer timescales, the study is also supplemented with two-epoch data from SDSS-I/II. The SDSS-I/II data include an additional $2854$ quasars with MgII and 572 quasars with Hbeta. The MgII emission line is significantly variable ($\Delta f/f$ 10% on 100-day timescales), a necessary prerequisite for its use for reverberation mapping studies. The data also confirm that continuum variability increases with timescale and decreases with luminosity, and the continuum light curves are consistent with a damped random-walk model on rest-frame timescales of $\gtrsim 5$ days. We compare the emission-line and continuum variability to investigate the structure of the broad-line region. Broad-line variability shows a shallower increase with timescale compared to the continuum emission, demonstrating that the broad-line transfer function is not a $\delta$-function. Hbeta is more variable than MgII (roughly by a factor of $1.5$), suggesting different excitation mechanisms, optical depths and/or geometrical configuration for each emission line. The ensemble spectroscopic variability measurements enabled by the SDSS-RM project have important consequences for future studies of reverberation mapping and black hole mass estimation of $1<z<2$ quasars.

We identify 885,503 type 1 quasar candidates to i<22 using the combination of optical and mid-IR photometry. Optical photometry is taken from the Sloan Digital Sky Survey-III: Baryon Oscillation Spectroscopic Survey (SDSS-III/BOSS), while mid-IR photometry comes from a combination of data from the Wide-Field Infrared Survey Explorer (WISE) "ALLWISE" data release and several large-area Spitzer Space Telescope fields. Selection is based on a Bayesian kernel density algorithm with a training sample of 157,701 spectroscopically-confirmed type-1 quasars with both optical and mid-IR data. Of the quasar candidates, 733,713 lack spectroscopic confirmation (and 305,623 are objects that we have not previously classified as photometric quasar candidates). These candidates include 7874 objects targeted as high probability potential quasars with 3.5<z<5 (of which 6779 are new photometric candidates). Our algorithm is more complete to z>3.5 than the traditional mid-IR selection "wedges" and to 2.2<z<3.5 quasars than the SDSS-III/BOSS project. Number counts and luminosity function analysis suggests that the resulting catalog is relatively complete to known quasars and is identifying new high-z quasars at z>3. This catalog paves the way for luminosity-dependent clustering investigations of large numbers of faint, high-redshift quasars and for further machine learning quasar selection using Spitzer and WISE data combined with other large-area optical imaging surveys.

We report the discovery of rapid variations of a high-velocity CIV broad absorption line trough in the quasar SDSS J141007.74+541203.3. This object was intensively observed in 2014 as a part of the Sloan Digital Sky Survey Reverberation Mapping Project, during which 32 epochs of spectroscopy were obtained with the Baryon Oscillation Spectroscopic Survey spectrograph. We observe significant (>4sigma) variability in the equivalent width of the broad (~4000 km/s wide) CIV trough on rest-frame timescales as short as 1.20 days (~29 hours), the shortest broad absorption line variability timescale yet reported. The equivalent width varied by ~10% on these short timescales, and by about a factor of two over the duration of the campaign. We evaluate several potential causes of the variability, concluding that the most likely cause is a rapid response to changes in the incident ionizing continuum. If the outflow is at a radius where the recombination rate is higher than the ionization rate, the timescale of variability places a lower limit on the density of the absorbing gas of n_e > 3.9 x 10^5 cm^-3. The broad absorption line variability characteristics of this quasar are consistent with those observed in previous studies of quasars, indicating that such short-term variability may in fact be common and thus can be used to learn about outflow characteristics and contributions to quasar/host-galaxy feedback scenarios.

We present host stellar velocity dispersion measurements for a sample of 88 broad-line quasars at 0.1<z<1 (46 at z>0.6) from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. High signal-to-noise ratio coadded spectra (average S/N~30 per 69 km/s pixel) from SDSS-RM allowed decomposition of the host and quasar spectra, and measurement of the host stellar velocity dispersions and black hole (BH) masses using the single-epoch (SE) virial method. The large sample size and dynamic range in luminosity (L5100=10^(43.2-44.7) erg/s) lead to the first clear detection of a correlation between SE virial BH mass and host stellar velocity dispersion far beyond the local universe. However, the observed correlation is significantly flatter than the local relation, suggesting that there are selection biases in high-z luminosity-threshold quasar samples for such studies. Our uniform sample and analysis enable an investigation of the redshift evolution of the M-sigma relation free of caveats by comparing different samples/analyses at disjoint redshifts. We do not observe evolution of the M-sigma relation in our sample, up to z~1, but there is an indication that the relation flattens towards higher redshifts. Coupled with the increasing threshold luminosity with redshift in our sample, this again suggests certain selection biases are at work, and simple simulations demonstrate that a constant M-sigma relation is favored to z~1. Our results highlight the scientific potential of deep coadded spectroscopy from quasar monitoring programs, and offer a new path to probe the co-evolution of BHs and galaxies at earlier times.

We present Keck-I MOSFIRE spectroscopy in the Y and H bands of GDN-8231, a massive, compact, star-forming galaxy (SFG) at a redshift $z\sim1.7$. Its spectrum reveals both H$_{\alpha}$ and [NII] emission lines and strong Balmer absorption lines. The H$_{\alpha}$ and Spitzer MIPS 24 $\mu$m fluxes are both weak, thus indicating a low star formation rate of SFR $\lesssim5-10$ M$_{\odot}$ yr$^{-1}$. This, added to a relatively young age of $\sim700$ Myr measured from the absorption lines, provides the first direct evidence for a distant galaxy being caught in the act of rapidly shutting down its star formation. Such quenching allows GDN-8231 to become a compact, quiescent galaxy, similar to 3 other galaxies in our sample, by $z\sim1.5$. Moreover, the color profile of GDN-8231 shows a bluer center, consistent with the predictions of recent simulations for an early phase of inside-out quenching. Its line-of-sight velocity dispersion for the gas, $\sigma^{\rm{gas}}_{\!_{\rm LOS}}=127\pm32$ km s$^{-1}$, is nearly 40% smaller than that of its stars, $\sigma^{\star}_{\!_{\rm LOS}}=215\pm35$ km s$^{-1}$. High-resolution hydro-simulations of galaxies explain such apparently colder gas kinematics of up to a factor of $\sim1.5$ with rotating disks being viewed at different inclinations and/or centrally concentrated star-forming regions. A clear prediction is that their compact, quiescent descendants preserve some remnant rotation from their star-forming progenitors.

We combine HST/WFC3 imaging and G141 grism observations from the CANDELS and 3D-HST surveys to produce a catalog of grism spectroscopic redshifts for galaxies in the CANDELS/GOODS-South field. The WFC3/G141 grism spectra cover a wavelength range of 1.1<lambda<1.7 microns with a resolving power of R~130 for point sources, thus providing rest-frame optical spectra for galaxies out to z~3.5. The catalog is selected in the H-band (F160W) and includes both galaxies with and without previously published spectroscopic redshifts. Grism spectra are extracted for all H-band detected galaxies with H<24 and a CANDELS photometric redshift z_phot > 0.6. The resulting spectra are visually inspected to identify emission lines and redshifts are determined using cross-correlation with empirical spectral templates. To establish the accuracy of our redshifts, we compare our results against high-quality spectroscopic redshifts from the literature. Using a sample of 411 control galaxies, this analysis yields a precision of sigma_NMAD=0.0028 for the grism-derived redshifts, which is consistent with the accuracy reported by the 3D-HST team. Our final catalog covers an area of 153 square arcmin and contains 1019 redshifts for galaxies in GOODS-S. Roughly 60% (608/1019) of these redshifts are for galaxies with no previously published spectroscopic redshift. These new redshifts span a range of 0.677 < z < 3.456 and have a median redshift of z=1.282. The catalog contains a total of 234 new redshifts for galaxies at z>1.5. In addition, we present 20 galaxy pair candidates identified for the first time using the grism redshifts in our catalog, including four new galaxy pairs at z~2, nearly doubling the number of such pairs previously identified.

We measure the location and evolutionary vectors of 69 Herschel-detected broad-line active galactic nuclei (BLAGNs) in the M_BH-M_* plane. BLAGNs are selected from the COSMOS and CDF-S fields, and span the redshift range 0.2< z<2.1. Black-hole masses are calculated using archival spectroscopy and single-epoch virial mass estimators, and galaxy total stellar masses are calculated by fitting the spectral energy distribution (subtracting the BLAGN component). The mass-growth rates of both the black hole and galaxy are calculated using Chandra/XMM-Newton X-ray and Herschel far-infrared data, reliable measures of the BLAGN accretion and galaxy star formation rates, respectively. We use Monte Carlo simulations to account for biases in our sample, due to both selection limits and the steep slope of the massive end of the galaxy stellar-mass distribution. We find our sample is consistent with no evolution in the M_BH-M_* relation from z~2 to z~0. BLAGNs and their host galaxies which lie off the black hole mass galaxy total stellar mass relation tend to have evolutionary vectors anti-correlated with their mass ratios: that is, galaxies with over-massive (under-massive) black holes tend to have a low (high) ratio of the specific accretion rate to the specific star formation rate. We also use the measured growth rates to estimate the preferred AGN duty cycle for our galaxies to evolve most consistently onto the local M_BH-M_Bul relation. Under reasonable assumptions of exponentially declining star formation histories, the data suggest a non-evolving (no more than a factor of a few) BLAGN duty cycle among star-forming galaxies of 10% (1sigma range of 1-42% at z<1 and 2-34% at z>1).

We present a new study investigating whether active galactic nuclei (AGN) beyond the local universe are preferentially fed via large-scale bars. Our investigation combines data from Chandra and Galaxy Zoo: Hubble (GZH) in the AEGIS, COSMOS, and GOODS-S surveys to create samples of face-on, disc galaxies at 0.2 < z < 1.0. We use a novel method to robustly compare a sample of 120 AGN host galaxies, defined to have 10^42 erg/s < L_X < 10^44 erg/s, with inactive control galaxies matched in stellar mass, rest-frame colour, size, Sersic index, and redshift. Using the GZH bar classifications of each sample, we demonstrate that AGN hosts show no statistically significant enhancement in bar fraction or average bar likelihood compared to closely-matched inactive galaxies. In detail, we find that the AGN bar fraction cannot be enhanced above the control bar fraction by more than a factor of two, at 99.7% confidence. We similarly find no significant difference in the AGN fraction among barred and non-barred galaxies. Thus we find no compelling evidence that large-scale bars directly fuel AGN at 0.2<z<1.0. This result, coupled with previous results at z=0, implies that moderate-luminosity AGN have not been preferentially fed by large-scale bars since z=1. Furthermore, given the low bar fractions at z>1, our findings suggest that large-scale bars have likely never directly been a dominant fueling mechanism for supermassive black hole growth.

Quenched galaxies at z>2 are nearly all very compact relative to z~0, suggesting a physical connection between high stellar density and efficient, rapid cessation of star-formation. We present restframe UV spectra of Lyman-break galaxies (LBGs) at z~3 selected to be candidate progenitors of quenched galaxies at z~2 based on their compact restframe optical sizes and high surface density of star-formation. We compare their UV properties to those of more extended LBGs of similar mass and star formation rate (non-candidates). We find that candidate progenitors have faster ISM gas velocities and higher equivalent widths of interstellar absorption lines, implying larger velocity spread among absorbing clouds. Candidates deviate from the relationship between equivalent widths of Lyman-alpha and interstellar absorption lines in that their Lyman-alpha emission remains strong despite high interstellar absorption, possibly indicating that the neutral HI fraction is patchy such that Lyman-alpha photons can escape. We detect stronger CIV P-Cygni features (emission and absorption) and HeII emission in candidates, indicative of larger populations of metal rich Wolf-Rayet stars compared to non-candidates. The faster bulk motions, broader spread of gas velocity, and Lyman-alpha properties of candidates are consistent with their ISM being subject to more energetic feedback than non-candidates. Together with their larger metallicity (implying more evolved star-formation activity) this leads us to propose, if speculatively, that they are likely to quench sooner than non-candidates, supporting the validity of selection criteria used to identify them as progenitors of z~2 passive galaxies. We propose that massive, compact galaxies undergo more rapid growth of stellar mass content, perhaps because the gas accretion mechanisms are different, and quench sooner than normally-sized LBGs at these early epochs.

Observations have shown that supermassive black holes in nearby elliptical galaxies correlate tightly with the stellar velocity dispersion (the $\MBH - \sigma$ relation) and the stellar mass (the $\MBH - \Mhost$ relation) of their host spheroids. However, the origin of these correlations remains ambiguous. In a previous paper by Zhu et al., we proposed a model which links the M-$\sigma$ relation to the the dynamical state of the system and the $\MBH - \Mhost$ relation to the self-regulation of galaxy growth. To test this model, we compile a sample of observed galaxies with different properties and examine the dependence of the above correlations on these parameters. We find that galaxies that satisfy the the $\MBH - \sigma$ correlation appear to have reached virial equilibrium, as indicated by the ratio between kinetic energy and gravitational potential, 2K/U $\sim$ 1. Furthermore, the ratio of black hole accretion rate to star formation rate remains nearly constant, BHAR /SFR $\sim$ $10^{-3}$, in active galaxies over a wide range of mass in the redshift range z=0 - 3. These results confirm our theoretical model that the observed correlations have different origins: the $\MBH - \sigma$ relation may result from galaxy relaxation, while the $\MBH$ - $\Mhost$ relation may be due to self-regulated black hole accretion and star formation in galaxies.

We use near-IR grism spectroscopy from the Hubble Space Telescope to examine the strength of [Ne~III] 3869 relative to H-beta, [O~II] 3727 and [O~III] 5007 in 236 low mass (7.5 < log (Mstar/Msolar) < 10.5) star-forming galaxies in the redshift range 1.90 < z < 2.35. By stacking the data by stellar mass, we show that the [Ne~III]/[O~II] ratios of the z ~ 2 universe are marginally higher than those seen in a comparable set of local SDSS galaxies, and that [Ne~III]/[O~III] is enhanced by ~0.2 dex. We consider the possible explanations for this ~4-sigma result, including higher oxygen depletion out of the gas-phase, denser H~II regions, higher production of Ne22 via Wolf-Rayet stars, and the existence of a larger population of X-ray obscured AGN at z ~ 2 compared to z ~ 0. None of these simple scenarios, alone, are favored to explain the observed line ratios. We conclude by suggesting several avenues of future observations to further explore the mystery of enhanced [Ne~III] emission.

We present photometric redshifts and associated probability distributions for all detected sources in the Extended Chandra Deep Field South (ECDFS). The work makes use of the most up-to-date data from the Cosmic Assembly Near-IR Deep Legacy Survey (CANDELS) and the Taiwan ECDFS Near-Infrared Survey (TENIS) in addition to other data. We also revisit multi-wavelength counterparts for published X-ray sources from the 4Ms-CDFS and 250ks-ECDFS surveys, finding reliable counterparts for 1207 out of 1259 sources ($\sim 96\%$). Data used for photometric redshifts include intermediate-band photometry deblended using the TFIT method, which is used for the first time in this work. Photometric redshifts for X-ray source counterparts are based on a new library of AGN/galaxy hybrid templates appropriate for the faint X-ray population in the CDFS. Photometric redshift accuracy for normal galaxies is 0.010 and for X-ray sources is 0.014, and outlier fractions are $4\%$ and $5.4\%$ respectively. The results within the CANDELS coverage area are even better as demonstrated both by spectroscopic comparison and by galaxy-pair statistics. Intermediate-band photometry, even if shallow, is valuable when combined with deep broad-band photometry. For best accuracy, templates must include emission lines.

We study the relationship between the structure and star-formation rate (SFR) of X-ray selected low and moderate luminosity active galactic nuclei (AGNs) in the two Chandra Deep Fields, using Hubble Space Telescope imaging from the Cosmic Assembly Near Infrared Extragalactic Legacy Survey (CANDELS) and deep far-infrared maps from the PEP+GOODS-Herschel survey. We derive detailed distributions of structural parameters and FIR luminosities from carefully constructed control samples of galaxies, which we then compare to those of the AGNs. At z~1, AGNs show slightly diskier light profiles than massive inactive (non-AGN) galaxies, as well as modestly higher levels of gross galaxy disturbance (as measured by visual signatures of interactions and clumpy structure). In contrast, at z~2, AGNs show similar levels of galaxy disturbance as inactive galaxies, but display a red central light enhancement, which may arise due to a more pronounced bulge in AGN hosts or due to extinguished nuclear light. We undertake a number of tests of these alternatives, but our results do not strongly favour one interpretation over the other. The mean SFR and its distribution among AGNs and inactive galaxies are similar at z>1.5. At z<1, however, clear and significant enhancements are seen in the SFRs of AGNs with bulge-dominated light profiles. These trends suggest an evolution in the relation between nuclear activity and host properties with redshift, towards a minor role for mergers and interactions at z>1.5.

The Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM) is a dedicated multi-object RM experiment that has spectroscopically monitored a sample of 849 broad-line quasars in a single 7 deg$^2$ field with the SDSS-III BOSS spectrograph. The RM quasar sample is flux-limited to i_psf=21.7 mag, and covers a redshift range of 0.1<z<4.5. Optical spectroscopy was performed during 2014 Jan-Jul dark/grey time, with an average cadence of ~4 days, totaling more than 30 epochs. Supporting photometric monitoring in the g and i bands was conducted at multiple facilities including the CFHT and the Steward Observatory Bok telescopes in 2014, with a cadence of ~2 days and covering all lunar phases. The RM field (RA, DEC=14:14:49.00, +53:05:00.0) lies within the CFHT-LS W3 field, and coincides with the Pan-STARRS 1 (PS1) Medium Deep Field MD07, with three prior years of multi-band PS1 light curves. The SDSS-RM 6-month baseline program aims to detect time lags between the quasar continuum and broad line region (BLR) variability on timescales of up to several months (in the observed frame) for ~10% of the sample, and to anchor the time baseline for continued monitoring in the future to detect lags on longer timescales and at higher redshift. SDSS-RM is the first major program to systematically explore the potential of RM for broad-line quasars at z>0.3, and will investigate the prospects of RM with all major broad lines covered in optical spectroscopy. SDSS-RM will provide guidance on future multi-object RM campaigns on larger scales, and is aiming to deliver more than tens of BLR lag detections for a homogeneous sample of quasars. We describe the motivation, design and implementation of this program, and outline the science impact expected from the resulting data for RM and general quasar science.

Limited studies have been performed on the radio-loud fraction in X-ray selected type 1 AGN samples. The consistency between various radio-loudness definitions also needs to be checked. We measure the radio-loudness of the 407 type 1 AGNs in the XMM-COSMOS quasar sample using nine criteria from the literature (six defined in the rest-frame and three defined in the observed frame): $R_L=\log(L_{5GHz}/L_B)$, $q_{24}=\log(L_{24\mu m}/L_{1.4GHz})$, $R_{uv}=\log(L_{5GHz}/L_{2500\AA})$, $R_{i}=\log(L_{1.4GHz}/L_i)$, $R_X=\log(\nu L_{\nu}(5GHz)/L_X)$, $P_{5GHz}=\log(P_{5GHz}(W/Hz/Sr))$, $R_{L,obs}=\log(f_{1.4GHz}/f_B)$ (observed frame), $R_{i,obs}=\log(f_{1.4GHz}/f_i)$ (observed frame), and $q_{24, obs}=\log(f_{24\mu m}/f_{1.4GHz})$ (observed frame). Using any single criterion defined in the rest-frame, we find a low radio-loud fraction of $\lesssim 5\%$ in the XMM-COSMOS type 1 AGN sample, except for $R_{uv}$. Requiring that any two criteria agree reduces the radio-loud fraction to $\lesssim 2\%$ for about 3/4 of the cases. The low radio-loud fraction cannot be simply explained by the contribution of the host galaxy luminosity and reddening. The $P_{5GHz}=\log(P_{5GHz}(W/Hz/Sr))$ gives the smallest radio-loud fraction. Two of the three radio-loud fractions from the criteria defined in the observed frame without k-correction ($R_{L,obs}$ and $R_{i,obs}$) are much larger than the radio-loud fractions from other criteria.

We use CANDELS imaging, 3D-HST spectroscopy, and Chandra X-ray data to investigate if active galactic nuclei (AGNs) are preferentially fueled by violent disk instabilities funneling gas into galaxy centers at 1.3<z<2.4. We select galaxies undergoing gravitational instabilities using the number of clumps and degree of patchiness as proxies. The CANDELS visual classification system is used to identify 44 clumpy disk galaxies, along with mass-matched comparison samples of smooth and intermediate morphology galaxies. We note that, despite being being mass-matched and having similar star formation rates, the smoother galaxies tend to be smaller disks with more prominent bulges compared to the clumpy galaxies. The lack of smooth extended disks is probably a general feature of the z~2 galaxy population, and means we cannot directly compare with the clumpy and smooth extended disks observed at lower redshift. We find that z~2 clumpy galaxies have slightly enhanced AGN fractions selected by integrated line ratios (in the mass-excitation method), but the spatially resolved line ratios indicate this is likely due to extended phenomena rather than nuclear AGNs. Meanwhile the X-ray data show that clumpy, smooth, and intermediate galaxies have nearly indistinguishable AGN fractions derived from both individual detections and stacked non-detections. The data demonstrate that AGN fueling modes at z~1.85 - whether violent disk instabilities or secular processes - are as efficient in smooth galaxies as they are in clumpy galaxies.

Post-starbursts are galaxies in transition from the blue cloud to the red sequence. Although they are rare today, integrated over time they may be an important pathway to the red sequence. This work uses SDSS, GALEX, and WISE observations to identify the evolutionary sequence from starbursts to fully quenched post-starbursts in the narrow mass range $\log M(M_\odot) = 10.3-10.7$, and identifies "transiting" post-starbursts which are intermediate between these two populations. In this mass range, $\sim 0.3\%$ of galaxies are starbursts, $\sim 0.1\%$ are quenched post-starbursts, and $\sim 0.5\%$ are the transiting types in between. The transiting post-starbursts have stellar properties that are predicted for fast-quenching starbursts and morphological characteristics that are already typical of early-type galaxies. The AGN fraction, as estimated from optical line ratios, of these post-starbursts is about 3 times higher ($\gtrsim 36 \pm 8 \%$) than that of normal star-forming galaxies of the same mass, but there is a significant delay between the starburst phase and the peak of nuclear optical AGN activity (median age difference of $\gtrsim 200 \pm 100$ Myr), in agreement with previous studies. The time delay is inferred by comparing the broad-band near NUV-to-optical photometry with stellar population synthesis models. We also find that starbursts and post-starbursts are significantly more dust-obscured than normal star-forming galaxies in the same mass range. About $20\%$ of the starbursts and $15\%$ of the transiting post-starbursts can be classified as the "Dust-Obscured Galaxies" (DOGs), while only $0.8\%$ of normal galaxies are DOGs.The time delay between the starburst phase and AGN activity suggests that AGN do not play a primary role in the original quenching of starbursts but may be responsible for quenching later low-level star formation during the post-starburst phase.

We consider how the profile and multi-year variability properties of a large sample of C IV Broad Absorption Line (BAL) troughs change when BALs from Si IV and/or Al III are present at corresponding velocities, indicating that the line-of-sight intercepts at least some lower ionization gas. We derive a number of observational results for C IV BALs separated according to the presence or absence of accompanying lower ionization transitions, including measurements of composite profile shapes, equivalent width (EW), characteristic velocities, composite variation profiles, and EW variability. We also measure the correlations between EW and fractional-EW variability for C IV, Si IV, and Al III. Our measurements reveal the basic correlated changes between ionization level, kinematics, and column density expected in accretion-disk wind models; e.g., lines-of-sight including lower ionization material generally show deeper and broader C IV troughs that have smaller minimum velocities and that are less variable. Many C IV BALs with no accompanying Si IV or Al III BALs may have only mild or no saturation.

We present near-infrared spectroscopy of a sample of 22 Extreme Emission Line Galaxies at redshifts 1.3 < z < 2.3, confirming that these are low-mass (M* = 10^8 - 10^9 M_sun) galaxies undergoing intense starburst episodes (M*/SFR ~ 10-100 Myr). The sample is selected by [O III] or H{\alpha} emission line flux and equivalent width using near-infrared grism spectroscopy from the 3D-HST survey. High-resolution NIR spectroscopy is obtained with LBT/LUCI and VLT/X-SHOOTER. The [O III]/H{\beta} line ratio is high (> 5) and [N II]/H{\alpha} is always significantly below unity, which suggests a low gas-phase metallicity. We are able to determine gas-phase metallicities for 7 of our objects using various strong-line methods, with values in the range 0.05-0.30 Z_sun and with a median of 0.15 Z_sun; for 3 of these objects we detect [O III]{\lambda}4363 which allows for a direct constraint on the metallicity. The velocity dispersion, as measured from the nebular emission lines, is typically ~50 km/s. Combined with the observed star-forming activity, the Jeans and Toomre stability criteria imply that the gas fraction must be large (> 2/3), consistent with the difference between our dynamical and stellar mass estimates. The implied gas depletion time scale (several hundred Myr) is substantially longer than the inferred mass-weighted ages (~50 Myr), which further supports the emerging picture that most stars in low-mass galaxies form in short, intense bursts of star formation.

We use Hubble Space Telescope near-IR grism spectroscopy to examine the H-beta line strengths of 260 star-forming galaxies in the redshift range 1.90 < z < 2.35. We show that at these epochs, the H-beta star formation rate (SFR) is a factor of ~1.8 higher than what would be expected from the systems' rest-frame UV flux density, suggesting a shift in the standard conversion between these quantities and star formation rate. We demonstrate that at least part of this shift can be attributed to metallicity, as H-beta is more greatly enhanced in systems with lower oxygen abundance. This offset must be considered when measuring the star formation rate history of the universe. We also show that the relation between stellar and nebular extinction in our z ~ 2 sample is consistent with that observed in the local universe.

We present Keck-I MOSFIRE near-infrared spectroscopy for a sample of 13 compact star-forming galaxies (SFGs) at redshift $2\leq z \leq2.5$ with star formation rates of SFR$\sim$100M$_{\odot}$ y$^{-1}$ and masses of log(M/M$_{\odot}$)$\sim10.8$. Their high integrated gas velocity dispersions of $\sigma_{\rm{int}}$=230$^{+40}_{-30}$ km s$^{-1}$, as measured from emission lines of H$_{\alpha}$ and [OIII], and the resultant M$_{\star}-\sigma_{\rm{int}}$ relation and M$_{\star}$$-$M$_{\rm{dyn}}$ all match well to those of compact quiescent galaxies at $z\sim2$, as measured from stellar absorption lines. Since log(M$_{\star}$/M$_{\rm{dyn}}$)$=-0.06\pm0.2$ dex, these compact SFGs appear to be dynamically relaxed and more evolved, i.e., more depleted in gas and dark matter ($<$13$^{+17}_{-13}$\%) than their non-compact SFG counterparts at the same epoch. Without infusion of external gas, depletion timescales are short, less than $\sim$300 Myr. This discovery adds another link to our new dynamical chain of evidence that compact SFGs at $z\gtrsim2$ are already losing gas to become the immediate progenitors of compact quiescent galaxies by $z\sim2$.

Emission line diagnostic diagrams probing the ionization sources in galaxies, such as the Baldwin-Phillips-Terlevich (BPT) diagram, have been used extensively to distinguish AGN from purely star-forming galaxies. Yet, they remain poorly understood at higher redshifts. We shed light on this issue with an empirical approach based on a z~0 reference sample built from ~300,000 SDSS galaxies, from which we mimic selection effects due to typical emission line detection limits at higher redshift. We combine this low-redshift reference sample with a simple prescription for luminosity evolution of the global galaxy population to predict the loci of high-redshift galaxies on the BPT and Mass-Excitation (MEx) diagnostic diagrams. The predicted bivariate distributions agree remarkably well with direct observations of galaxies out to z~1.5, including the observed stellar mass-metallicity (MZ) relation evolution. As a result, we infer that high-redshift star-forming galaxies are consistent with having "normal" ISM properties out to z~1.5, after accounting for selection effects and line luminosity evolution. Namely, their optical line ratios and gas-phase metallicities are comparable to that of low-redshift galaxies with equivalent emission-line luminosities. In contrast, AGN narrow-line regions may show a shift toward lower metallicities at higher redshift. While a physical evolution of the ISM conditions is not ruled out for purely star-forming galaxies, and may be more important starting at z>2, we find that reliably quantifying this evolution is hindered by selections effects. The recipes provided here may serve as a basis for future studies toward this goal. Code to predict the loci of galaxies on the BPT and MEx diagnostic diagrams, and the MZ relation as a function of emission line luminosity limits, is made publicly available.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) was a multi-cycle treasury program on the Hubble Space Telescope (HST) that surveyed a total area of ~0.25 deg^2 with ~900 HST orbits spread across 5 fields over 3 years. Within these survey images we discovered 65 supernovae (SN) of all types, out to z~2.5. We classify ~24 of these as Type Ia SN (SN Ia) based on host-galaxy redshifts and SN photometry (supplemented by grism spectroscopy of 6 SN). Here we present a measurement of the volumetric SN Ia rate as a function of redshift, reaching for the first time beyond z=2 and putting new constraints on SN Ia progenitor models. Our highest redshift bin includes detections of SN that exploded when the universe was only ~3 Gyr old and near the peak of the cosmic star-formation history. This gives the CANDELS high-redshift sample unique leverage for evaluating the fraction of SN Ia that explode promptly after formation (<500 Myr). Combining the CANDELS rates with all available SN Ia rate measurements in the literature we find that this prompt SN Ia fraction is fP=0.53 +0.09 -0.10 (stat) +0.10 -0.26 (sys), consistent with a delay time distribution that follows a simple t^{-1} power law for all times t>40 Myr. However, a mild tension is apparent between ground-based low-z surveys and space-based high-z surveys. In both CANDELS and the sister HST program CLASH, we find a low rate of SN Ia at z>1. This could be a hint that prompt progenitors are in fact relatively rare, accounting for only ~20% of all SN Ia explosions -- though further analysis and larger samples will be needed to examine that suggestion.

Apart from viewing-dependent obscuration, intrinsic broad-line emission from active galactic nuclei (AGNs) follows an evolutionary sequence: Type $1 \to 1.2/1.5 \to 1.8/1.9 \to 2$ as the accretion rate onto the central black hole is decreasing. This spectral evolution is controlled, at least in part, by the parameter $L_{\rm bol}/M^{2/3}$, where $L_{\rm bol}$ is the AGN bolometric luminosity and $M$ is the black hole mass. Both this dependence and the double-peaked profiles that emerge along the sequence arise naturally in the disk-wind scenario for the AGN broad-line region.

We took advantage of the observations carried out by XMM in the COSMOS field during 3.5 years, to study the long term variability of a large sample of AGN (638 sources), in a wide range of redshift (0.1<z<3.5) and X-ray luminosity ($10^{41}<$L(2-10)$<10^{45.5}$). Both a simple statistical method to asses the significance of variability, and the Normalized Excess Variance ($\sigma^{2}_{rms}$) parameter, where used to obtain a quantitative measurement of the variability. Variability is found to be prevalent in most AGN, whenever we have good statistic to measure it, and no significant differences between type-1 and type-2 AGN were found. A flat (slope -0.23+/-0.03) anti-correlation between $\sigma^{2}_{rms}$ and X-ray luminosity is found, when significantly variable sources are considered all together. When divided in three redshift bins, the anti-correlation becomes stronger and evolving with z, with higher redshift AGN being more variable. We prove however that this effect is due to the pre-selection of variable sources: considering all the sources with available $\sigma^{2}_{rms}$ measurement, the evolution in redshift disappears. For the first time we were also able to study the long term X-ray variability as a function of $M_{\rm BH}$ and Eddington ratio, for a large sample of AGN spanning a wide range of redshift. An anti-correlation between $\sigma^{2}_{rms}$ and $M_{\rm BH}$ is found, with the same slope of the anti-correlation between $\sigma^{2}_{rms}$ and X-ray luminosity, suggesting that the latter can be a byproduct of the former one. No clear correlation is found between $\sigma^{2}_{rms}$ and the Eddington ratio in our sample. Finally, no correlation is found between the X-ray $\sigma^{2}_{rms}$ and the optical variability.

The mid-infrared to ultraviolet (0.1 -- 10 $\mu m$) spectral energy distribution (SED) shapes of 407 X-ray-selected radio-quiet type 1 AGN in the wide-field ``Cosmic Evolution Survey" (COSMOS) have been studied for signs of evolution. For a sub-sample of 200 radio-quiet quasars with black hole mass estimates and host galaxy corrections, we studied their mean SEDs as a function of a broad range of redshift, bolometric luminosity, black hole mass and Eddington ratio, and compared them with the Elvis et al. (1994, E94) type 1 AGN mean SED. We found that the mean SEDs in each bin are closely similar to each other, showing no statistical significant evidence of dependence on any of the analyzed parameters. We also measured the SED dispersion as a function of these four parameters, and found no significant dependencies. The dispersion of the XMM-COSMOS SEDs is generally larger than E94 SED dispersion in the ultraviolet, which might be due to the broader ``window function'' for COSMOS quasars, and their X-ray based selection.

We analyze the star-forming and structural properties of 45 massive (log(M/Msun)>10) compact star-forming galaxies (SFGs) at 2<z<3 to explore whether they are progenitors of compact quiescent galaxies at z~2. The optical/NIR and far-IR Spitzer/Herschel colors indicate that most compact SFGs are heavily obscured. Nearly half (47%) host an X-ray bright AGN. In contrast, only about 10% of other massive galaxies at that time host AGNs. Compact SFGs have centrally-concentrated light profiles and spheroidal morphologies similar to quiescent galaxies, and are thus strikingly different from other SFGs. Most compact SFGs lie either within the SFR-M main sequence (65%) or below (30%), on the expected evolutionary path towards quiescent galaxies. These results show conclusively that galaxies become more compact before they lose their gas and dust, quenching star formation. Using extensive HST photometry from CANDELS and grism spectroscopy from the 3D-HST survey, we model their stellar populations with either exponentially declining (tau) star formation histories (SFHs) or physically-motivated SFHs drawn from semi-analytic models (SAMs). SAMs predict longer formation timescales and older ages ~2 Gyr, which are nearly twice as old as the estimates of the tau models. While both models yield good SED fits, SAM SFHs better match the observed slope and zero point of the SFR-M main sequence. Some low-mass compact SFGs (log(M/Msun)=10-10.6) have younger ages but lower sSFRs than that of more massive galaxies, suggesting that the low-mass galaxies reach the red sequence faster. If the progenitors of compact SFGs are extended SFGs, state-of-the-art SAMs show that mergers and disk instabilities are both able to shrink galaxies, but disk instabilities are more frequent (60% versus 40%) and form more concentrated galaxies. We confirm this result via high-resolution hydrodynamic simulations.

Spectroscopic observations from the Large Binocular Telescope and the Very Large Telescope reveal kinematically narrow lines (~50 km/s) for a sample of 14 Extreme Emission Line Galaxies (EELGs) at redshifts 1.4 < z < 2.3. These measurements imply that the total dynamical masses of these systems are low (< 3x10^9 M_sun). Their large [O III] 5007 equivalent widths (500-1100 Angstroms) and faint blue continuum emission imply young ages of 10-100 Myr and stellar masses of 10^8-10^9 M_sun, confirming the presence of a violent starburst. The dynamical masses represent the first such determinations for low-mass galaxies at z > 1. The stellar mass formed in this vigorous starburst phase represents a large fraction of the total (dynamical) mass, without a significantly massive underlying population of older stars. The occurrence of such intense events in shallow potentials strongly suggests that supernova-driven winds must be of critical importance in the subsequent evolution of these systems.

Quasi-stellar objects (QSOs) occur in galaxies in which supermassive black holes (SMBHs) are growing substantially through rapid accretion of gas. Many popular models of the co-evolutionary growth of galaxies and SMBHs predict that QSOs are also sites of substantial recent star formation, mediated by important processes, such as major mergers, which rapidly transform the nature of galaxies. A detailed study of the star-forming properties of QSOs is a critical test of such models. We present a far-infrared Herschel/PACS study of the mean star formation rate (SFR) of a sample of spectroscopically observed QSOs to z~2 from the COSMOS extragalactic survey. This is the largest sample to date of moderately luminous AGNs studied using uniform, deep far-infrared photometry. We study trends of the mean SFR with redshift, black hole mass, nuclear bolometric luminosity and specific accretion rate (Eddington ratio). To minimize systematics, we have undertaken a uniform determination of SMBH properties, as well as an analysis of important selection effects within spectroscopic QSO samples that influence the interpretation of SFR trends. We find that the mean SFRs of these QSOs are consistent with those of normal massive star-forming galaxies with a fixed scaling between SMBH and galaxy mass at all redshifts. No strong enhancement in SFR is found even among the most rapidly accreting systems, at odds with several co-evolutionary models. Finally, we consider the qualitative effects on mean SFR trends from different assumptions about the star-forming properties of QSO hosts and redshift evolution of the SMBH-galaxy relationship. While limited currently by uncertainties, valuable constraints on AGN-galaxy co-evolution can emerge from our approach.

We define a quasar-galaxy mixing diagram using the slopes of their spectral energy distributions (SEDs) from 1\mu m to 3000\AA\ and from 1\mu m to 3\mu m in the rest frame. The mixing diagram can easily distinguish among quasar-dominated, galaxy-dominated and reddening-dominated SED shapes. By studying the position of the 413 XMM selected Type 1 AGN in the wide-field "Cosmic Evolution Survey" (COSMOS) in the mixing diagram, we find that a combination of the Elvis et al. (1994, hereafter E94) quasar SED with various contributions from galaxy emission and some dust reddening is remarkably effective in describing the SED shape from 0.3-3\mu m for large ranges of redshift, luminosity, black hole mass and Eddington ratio of type 1 AGN. In particular, the location in the mixing diagram of the highest luminosity AGN is very close (within 1\sigma) to that of the E94 SED. The mixing diagram can also be used to estimate the host galaxy fraction and reddening in quasar. We also show examples of some outliers which might be AGN in different evolutionary stages compared to the majority of AGN in the quasar-host galaxy co-evolution cycle.

We have made a serendipitous discovery of a massive cD galaxy at z=1.096 in a candidate rich cluster in the HUDF area of GOODS-South. This brightest cluster galaxy is the most distant cD galaxy confirmed to date. Ultra-deep HST/WFC3 images reveal an extended envelope starting from ~10 kpc and reaching ~70 kpc in radius along the semi-major axis. The spectral energy distributions indicate that both its inner component and outer envelope are composed of an old, passively-evolving stellar population. The cD galaxy lies on the same mass-size relation as the bulk of quiescent galaxies at similar redshifts. The cD galaxy has a higher stellar mass surface density but a similar velocity dispersion to those of more-massive, nearby cDs. If the cD galaxy is one of the progenitors of today's more massive cDs, its size and stellar mass have had to increase on average by factors of $3.4\pm1.1$ and $3.3\pm1.3$ over the past ~8 Gyrs, respectively. Such increases in size and stellar mass without being accompanied by significant increases in velocity dispersion are consistent with evolutionary scenarios driven by both major and minor dry mergers. If such cD envelopes originate from dry mergers, our discovery of even one example proves that some BCGs entered the dry merger phase at epochs earlier than z=1. Our data match theoretical models which predict that the continuance of dry mergers at z<1 can result in structures similar to those of massive cD galaxies seen today. Moreover, our discovery is a surprise given that the extreme depth of the HUDF is essential to reveal such an extended cD envelope at z>1 and, yet, the HUDF covers only a minuscule region of sky. Adding that cDs are rare, Our serendipitous discovery hints that such cDs may be more common than expected. [Abridged]

We characterize the incidence of active galactic nuclei (AGNs) is 0.3 < z < 1 star-forming galaxies by applying multi-wavelength AGN diagnostics (X-ray, optical, mid-infrared, radio) to a sample of galaxies selected at 70-micron from the Far-Infrared Deep Extragalactic Legacy survey (FIDEL). Given the depth of FIDEL, we detect "normal" galaxies on the specific star formation rate (sSFR) sequence as well as starbursting systems with elevated sSFR. We find an overall high occurrence of AGN of 37+/-3%, more than twice as high as in previous studies of galaxies with comparable infrared luminosities and redshifts but in good agreement with the AGN fraction of nearby (0.05 < z < 0.1) galaxies of similar infrared luminosities. The more complete census of AGNs comes from using the recently developed Mass-Excitation (MEx) diagnostic diagram. This optical diagnostic is also sensitive to X-ray weak AGNs and X-ray absorbed AGNs, and reveals that absorbed active nuclei reside almost exclusively in infrared-luminous hosts. The fraction of galaxies hosting an AGN appears to be independent of sSFR and remains elevated both on the sSFR sequence and above. In contrast, the fraction of AGNs that are X-ray absorbed increases substantially with increasing sSFR, possibly due to an increased gas fraction and/or gas density in the host galaxies.

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

We present some of the first science data with the new Keck/MOSFIRE instrument to test the effectiveness of different AGN/SF diagnostics at z~1.5. MOSFIRE spectra were obtained in three H-band multi-slit masks in the GOODS-S field, resulting in two hour exposures of 36 emission-line galaxies. We compare X-ray data with the traditional emission-line ratio diagnostics and the alternative mass-excitation and color-excitation diagrams, combining new MOSFIRE infrared data with previous HST/WFC3 infrared spectra (from the 3D-HST survey) and multiwavelength photometry. We demonstrate that a high [OIII]/Hb ratio is insufficient as an AGN indicator at z>1. For the four X-ray detected galaxies, the classic diagnostics ([OIII]/Hb vs. [NII]/Ha and [SII]/Ha) remain consistent with X-ray AGN/SF classification. The X-ray data also suggest that "composite" galaxies (with intermediate AGN/SF classification) host bona-fide AGNs. Nearly 2/3 of the z~1.5 emission-line galaxies have nuclear activity detected by either X-rays or the classic diagnostics. Compared to the X-ray and line ratio classifications, the mass-excitation method remains effective at z>1, but we show that the color-excitation method requires a new calibration to successfully identify AGNs at these redshifts.

We present the first quantified, statistical map of broad-line active galactic nucleus (AGN) frequency with host galaxy color and stellar mass in nearby (0.01 < z < 0.11) galaxies. Aperture photometry and z-band concentration measurements from the Sloan Digital Sky Survey (SDSS) are used to dis- entangle AGN and galaxy emission, resulting in estimates of uncontaminated galaxy rest-frame color, luminosity, and stellar mass. Broad-line AGNs are distributed throughout the blue cloud and green valley at a given stellar mass, and are much rarer in quiescent (red sequence) galaxies. This is in contrast to the published host galaxy properties of weaker narrow-line AGNs, indicating that broad-line AGNs occur during a different phase in galaxy evolution. More luminous broad-line AGNs have bluer host galaxies, even at fixed mass, suggesting that the same processes that fuel nuclear activity also efficiently form stars. The data favor processes that simultaneously fuel both star formation activity and rapid supermassive black hole accretion. If AGNs cause feedback on their host galaxies in the nearby universe, the evidence of galaxy-wide quenching must be delayed until after the broad-line AGN phase.

The "Cosmic Evolution Survey" (COSMOS) enables the study of the Spectral Energy Distributions (SEDs) of Active Galactic Nuclei (AGN) because of the deep coverage and rich sampling of frequencies from X-ray to radio. Here we present a SED catalog of 413 X-ray (\xmm) selected type 1 (emission line FWHM$>2000$ km s$^{-1}$) AGN with Magellan, SDSS or VLT spectrum. The SEDs are corrected for the Galactic extinction, for broad emission line contributions, constrained variability, and for host galaxy contribution. We present the mean SED and the dispersion SEDs after the above corrections in the rest frame 1.4 GHz to 40 keV, and show examples of the variety of SEDs encountered. In the near-infrared to optical (rest frame $\sim 8\mu m$-- 4000\AA), the photometry is complete for the whole sample and the mean SED is derived from detections only. Reddening and host galaxy contamination could account for a large fraction of the observed SED variety. The SEDs are all available on-line.

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

We combine high-resolution HST/WFC3 images with multi-wavelength photometry to track the evolution of structure and activity of massive (log(M*) > 10) galaxies at redshifts z = 1.4 - 3 in two fields of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). We detect compact, star-forming galaxies (cSFGs) whose number densities, masses, sizes, and star formation rates qualify them as likely progenitors of compact, quiescent, massive galaxies (cQGs) at z = 1.5 - 3. At z > 2 most cSFGs have specific star-formation rates (sSFR = 10^-9 yr^-1) half that of typical, massive SFGs at the same epoch, and host X-ray luminous AGN 30 times (~30%) more frequently. These properties suggest that cSFGs are formed by gas-rich processes (mergers or disk-instabilities) that induce a compact starburst and feed an AGN, which, in turn, quench the star formation on dynamical timescales (few 10^8 yr). The cSFGs are continuously being formed at z = 2 - 3 and fade to cQGs by z = 1.5. After this epoch, cSFGs are rare, thereby truncating the formation of new cQGs. Meanwhile, down to z = 1, existing cQGs continue to enlarge to match local QGs in size, while less-gas-rich mergers and other secular mechanisms shepherd (larger) SFGs as later arrivals to the red sequence. In summary, we propose two evolutionary scenarios of QG formation: an early (z > 2), fast-formation path of rapidly-quenched cSFGs that evolve into cQGs that later enlarge within the quiescent phase, and a slow, late-arrival (z < 2) path for SFGs to form QGs without passing through a compact state.

Bolometric luminosities and Eddington ratios of both X-ray selected broad-line (Type-1) and narrow-line (Type-2) AGN from the XMM-Newton survey in the COSMOS field are presented. The sample is composed by 929 AGN (382 Type-1 AGN and 547 Type-2 AGN) and it covers a wide range of redshifts, X-ray luminosities and absorbing column densities. About 65% of the sources are spectroscopically identified as either Type-1 or Type-2 AGN (83% and 52% respectively), while accurate photometric redshifts are available for the rest of the sample. The study of such a large sample of X-ray selected AGN with a high quality multi-wavelength coverage from the far-infrared (now with the inclusion of Herschel data at 100 micron and 160 micron) to the optical-UV allows us to obtain accurate estimates of bolometric luminosities, bolometric corrections and Eddington ratios. The kbol-Lbol relations derived in the present work are calibrated for the first time against a sizable AGN sample, and rely on observed redshifts, X-ray luminosities and column density distributions. We find that kbol is significantly lower at high Lbol with respect to previous estimates by Marconi et al. (2004) and Hopkins et al. (2007). Black hole masses and Eddington ratios are available for 170 Type-1 AGN, while black hole masses for Type-2 AGN are computed for 481 objects using the black hole mass-stellar mass relation and the morphological information. We confirm a trend between kbol and lambda_Edd, with lower hard X-ray bolometric corrections at lower Eddington ratios for both Type-1 and Type-2 AGN. We find that, on average, Eddington ratio increases with redshift for all Types of AGN at any given Mbh, while no clear evolution with redshift is seen at any given Lbol.

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

We present Hubble Space Telescope Wide Field Camera 3 slitless grism spectroscopy of 28 emission-line galaxies at z~2, in the GOODS-S region of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). The high sensitivity of these grism observations, with 1-sigma detections of emission lines to f > 2.5x10^{-18} erg/s/cm^2, means that the galaxies in the sample are typically ~7 times less massive (median M_* = 10^{9.5} M_sun) than previously studied z~2 emission-line galaxies. Despite their lower mass, the galaxies have OIII/Hb ratios which are very similar to previously studied z~2 galaxies and much higher than the typical emission-line ratios of local galaxies. The WFC3 grism allows for unique studies of spatial gradients in emission lines, and we stack the two-dimensional spectra of the galaxies for this purpose. In the stacked data the OIII emission line is more spatially concentrated than the Hb emission line with 98.1 confidence. We additionally stack the X-ray data (all sources are individually undetected), and find that the average L(OIII)/L(0.5-10 keV) ratio is intermediate between typical z~0 obscured active galaxies and star-forming galaxies. Together the compactness of the stacked OIII spatial profile and the stacked X-ray data suggest that at least some of these low-mass, low-metallicity galaxies harbor weak active galactic nuclei.

Spitzer IRAC selection is a powerful tool for identifying luminous AGN. For deep IRAC data, however, the AGN selection wedges currently in use are heavily contaminated by star-forming galaxies, especially at high redshift. Using the large samples of luminous AGN and high-redshift star-forming galaxies in COSMOS, we redefine the AGN selection criteria for use in deep IRAC surveys. The new IRAC criteria are designed to be both highly complete and reliable, and incorporate the best aspects of the current AGN selection wedges and of infrared power-law selection while excluding high redshift star-forming galaxies selected via the BzK, DRG, LBG, and SMG criteria. At QSO-luminosities of log L(2-10 keV) (ergs/s) > 44, the new IRAC criteria recover 75% of the hard X-ray and IRAC-detected XMM-COSMOS sample, yet only 38% of the IRAC AGN candidates have X-ray counterparts, a fraction that rises to 52% in regions with Chandra exposures of 50-160 ks. X-ray stacking of the individually X-ray non-detected AGN candidates leads to a hard X-ray signal indicative of heavily obscured to mildly Compton-thick obscuration (log N_H (cm^-2) = 23.5 +/- 0.4). While IRAC selection recovers a substantial fraction of luminous unobscured and obscured AGN, it is incomplete to low-luminosity and host-dominated AGN.

We use a sample of active galaxies from the Cosmic Evolution Survey to show that host galaxy morphology is tied to the accretion rate and X-ray obscuration of its active galactic nucleus (AGN). Unobscured and rapidly accreting broad-line AGNs are more likely to be in spheroid-dominated hosts than weak or obscured AGNs, and obscured AGNs are more likely to have disturbed host galaxies. Much of the disagreement in previous work on the AGN-merger connection is likely due to each study probing AGNs with different obscuration and accretion properties. Only obscured AGNs seem to merger-driven, while weak AGNs are fed by stochastic processes in disks, and rapidly-accreting broad-line AGNs require massive bulges. Our observed "unified model" for AGN hosts fits with theoretical models for merger-driven AGN evolution, but is also consistent with steady-state AGN activity.

In order to understand environmental effects on star formation in high-redshift galaxies, we investigate the physical relationships between the star formation activity, stellar mass, and environment for z ~1.2 galaxies in the 2 deg^2 COSMOS field. We estimate star formation using the [OII] emission line and environment from the local galaxy density. Our analysis shows that for massive galaxies M_*>10^10 M_sun, the fraction of [OII] emitters in high-density environments is 1.7 times higher than in low-density environments, while the [OII] emitter fraction does not depend on environment for low-mass M_* < 10^10 M_sun galaxies. In order to understand what drives these trends, we investigate the role of companion galaxies in our sample. We find that the fraction of [OII] emitters in galaxies with companions is 2.4 times as high as that in galaxies without companions at M_* > 10^10 M_sun. In addition, massive galaxies are more likely to have companions in high-density environments. However, although the "number" of star forming galaxies increases for massive galaxies with close companions and in dense environments, the "average" star formation rate of star forming galaxies at a given mass is independent of environment and the presence/absence of a close companion. These results suggest that interactions and/or mergers in high-density environment could induce star formation in massive galaxies at z~1.2, increasing the fraction of star-forming galaxies with M_* > 10^10 M_sun.

We present a rest-frame spectral stacking analysis of ~1000 X-ray sources detected in the XMM-COSMOS field in order to investigate the iron K line properties of active galaxies beyond redshift z~1. In Type I AGN that have a typical X-ray luminosity of Lx~1.5e44 erg/s and z~1.6, the cold Fe K at 6.4 keV is weak (EW~0.05keV), in agreement with the known trend. In contrast, high-ionization lines of Fe XXV and Fe XXVI are pronounced. These high-ionization Fe K lines appear to have a connection with high accretion rates. While no broad Fe emission is detected in the total spectrum, it might be present, albeit at low significance, when the X-ray luminosity is restricted to the range below 3e44 erg/s, or when an intermediate range of Eddington ratio around 0.1 is selected. In Type II AGN, both cold and high-ionzation lines become weak with increasing X-ray luminosity. However, strong high-ionization Fe K (EW~0.3 keV) is detected in the spectrum of objects at z>2, while no 6.4 keV line is found. It is then found that the primary source of the high-ionization Fe K emission is those objects detected with Spitzer-MIPS at 24 micron. Given their median redshift of z=2.5, their bolometric luminosity is likely to reach 10^13 Lsun and the MIPS-detected emission most likely originates from hot dust heated by embedded AGN, probably accreting at high Eddington ratio. These properties match those of rapidly growing black holes in ultra-luminous infrared galaxies at the interesting epoch (z=2-3) of galaxy formation.

The Lick AGN Monitoring Project 2011 observing campaign was carried out over the course of 11 weeks in Spring 2011. Here we present the first results from this program, a measurement of the broad-line reverberation lag in the Seyfert 1 galaxy Mrk 50. Combining our data with supplemental observations obtained prior to the start of the main observing campaign, our dataset covers a total duration of 4.5 months. During this time, Mrk 50 was highly variable, exhibiting a maximum variability amplitude of a factor of 4 in the U-band continuum and a factor of 2 in the H-beta line. Using standard cross-correlation techniques, we find that H-beta and H-gamma lag the V-band continuum by tau_cen = 10.64(-0.93,+0.82) and 8.43(-1.28,+1.30) days, respectively, while the lag of He II 4686 is unresolved. The H-beta line exhibits a symmetric velocity-resolved reverberation signature with shorter lags in the high-velocity wings than in the line core, consistent with an origin in a broad-line region dominated by orbital motion rather than infall or outflow. Assuming a virial normalization factor of f=5.25, the virial estimate of the black hole mass is (3.2+-0.5)*10^7 solar masses. These observations demonstrate that Mrk 50 is among the most promising nearby active galaxies for detailed investigations of broad-line region structure and dynamics.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5\sigma\ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5\sigma\ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

We identify an abundant population of extreme emission line galaxies (EELGs) at redshift z~1.7 in the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) imaging from Hubble Space Telescope/Wide Field Camera 3 (HST/WFC3). 69 EELG candidates are selected by the large contribution of exceptionally bright emission lines to their near-infrared broad-band magnitudes. Supported by spectroscopic confirmation of strong [OIII] emission lines -- with rest-frame equivalent widths ~1000\AA -- in the four candidates that have HST/WFC3 grism observations, we conclude that these objects are galaxies with 10^8 Msol in stellar mass, undergoing an enormous starburst phase with M_*/(dM_*/dt) of only ~15 Myr. These bursts may cause outflows that are strong enough to produce cored dark matter profiles in low-mass galaxies. The individual star formation rates and the co-moving number density (3.7x10^-4 Mpc^-3) can produce in ~4 Gyr much of the stellar mass density that is presently contained in 10^8-10^9 Msol dwarf galaxies. Therefore, our observations provide a strong indication that many or even most of the stars in present-day dwarf galaxies formed in strong, short-lived bursts, mostly at z>1.

We present new constraints on the ratio of black hole (BH) mass to total galaxy stellar mass at 0.3<z<0.9 for a sample of 32 type-1 active galactic nuclei (AGNs) from the XMM-COSMOS survey covering the range M_BH/M_sun~10^(7.2--8.7). Virial M_BH estimates based on H_beta are available from the COSMOS Magellan/IMACS survey. We use high-resolution Hubble Space Telescope (HST) imaging to decompose the light of each type-1 AGN and host galaxy, and employ a specially-built mass-to-light ratio to estimate the stellar masses (M_*). The M_BH-M_* ratio shows a zero offset with respect to the local relation for galactic bulge masses, and we also find no evolution in the mass ratio M_BH/M_*=(1+z)^{0.02+-0.34} up to z~0.9. Interestingly, at the high-M_BH end there is a positive offset from the z=0 relation, which can be fully explained by a mass function bias with a cosmic scatter of 0.3, reaffirming that the intrinsic distribution is consistent with zero evolution. From our results we conclude that since z~0.9 no substantial addition of stellar mass is required: the decline in star formation rates and merger activity at z<1 support this scenario. Nevertheless, given that a significant fraction of these galaxies show a disk component, their bulges are indeed undermassive. This is a direct indication that for the last 7 Gyr the only essential mechanism required in order that these galaxies obey the z=0 relation is a redistribution of stellar mass to the bulge, likely driven by secular processes, i.e., internal instabilities and minor merging.

We explore the connection between black hole growth at the center of obscured quasars selected from the XMM-COSMOS survey and the physical properties of their host galaxies. We study a bolometric regime (<Lbol > 8 x 10^45 erg/s) where several theoretical models invoke major galaxy mergers as the main fueling channel for black hole accretion. We confirm that obscured quasars mainly reside in massive galaxies (Mstar>10^10 Msun) and that the fraction of galaxies hosting such powerful quasars monotonically increases with the stellar mass. We stress the limitation of the use of rest-frame color-magnitude diagrams as a diagnostic tool for studying galaxy evolution and inferring the influence that AGN activity can have on such a process. We instead use the correlation between star-formation rate and stellar mass found for star-forming galaxies to discuss the physical properties of the hosts. We find that at z ~1, ~62% of Type-2 QSOs hosts are actively forming stars and that their rates are comparable to those measured for normal star-forming galaxies. The fraction of star-forming hosts increases with redshift: ~71% at z ~2, and 100% at z ~3. We also find that the the evolution from z ~1 to z ~3 of the specific star-formation rate of the Type-2 QSO hosts is in excellent agreement with that measured for star-forming galaxies. From the morphological analysis, we conclude that most of the objects are bulge-dominated galaxies, and that only a few of them exhibit signs of recent mergers or disks. Finally, bulge-dominated galaxies tend to host Type-2 QSOs with low Eddington ratios (lambda<0.1), while disk-dominated or merging galaxies have at their centers BHs accreting at high Eddington ratios (lambda > 0.1).

Using HST/WFC3 imaging taken as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS), we examine the role that major galaxy mergers play in triggering active galactic nuclei (AGN) activity at z~2. Our sample consists of 72 moderate-luminosity (Lx ~ 1E42-1E44 erg/s) AGN at 1.5<z<2.5 that are selected using the 4 Msec Chandra observations in the Chandra Deep Field South, the deepest X-ray observations to date. Employing visual classifications, we have analyzed the rest-frame optical morphologies of the AGN host galaxies and compared them to a mass-matched control sample of 216 non-active galaxies at the same redshift. We find that most of the AGN reside in disk galaxies (51.4%), while a smaller percentage are found in spheroids (27.8%). Roughly 16.7% of the AGN hosts have highly disturbed morphologies and appear to be involved in a major merger or interaction, while most of the hosts (55.6%) appear relatively relaxed and undisturbed. These fractions are statistically consistent with the fraction of control galaxies that show similar morphological disturbances. These results suggest that the hosts of moderate-luminosity AGN are no more likely to be involved in an ongoing merger or interaction relative to non-active galaxies of similar mass at z~2. The high disk fraction observed among the AGN hosts also appears to be at odds with predictions that merger-driven accretion should be the dominant AGN fueling mode at z~2, even at moderate X-ray luminosities. Although we cannot rule out that minor mergers are responsible for triggering these systems, the presence of a large population of relatively undisturbed disk-like hosts suggests that secular processes play a greater role in fueling AGN activity at z~2 than previously thought.

We explore the evolution with redshift of the rest-frame colours and space densities of AGN hosts (relative to normal galaxies) to shed light on the dominant mechanism that triggers accretion onto supermassive black holes as a function of cosmic time. Data from serendipitous wide-area XMM surveys of the SDSS footprint (XMM/SDSS, Needles in the Haystack survey) are combined with Chandra deep observations in the AEGIS, GOODS-North and GOODS-South to compile uniformly selected samples of moderate luminosity X-ray AGN [L_X(2-10keV) = 1e41-1e44erg/s] at redshifts 0.1, 0.3 and 0.8. It is found that the fraction of AGN hosted by red versus blue galaxies does not change with redshift. Also, the X-ray luminosity density associated with either red or blue AGN hosts remains nearly constant since z=0.8. X-ray AGN represent a roughly fixed fraction of the space density of galaxies of given optical luminosity at all redshifts probed by our samples. In contrast the fraction of X-ray AGN among galaxies of a given stellar mass decreases with decreasing redshift. These findings suggest that the same process or combination of processes for fueling supermassive black holes are in operation in the last 5 Gyrs of cosmic time. The data are consistent with a picture in which the drop of the accretion power during that period (1dex since z=0.8) is related to the decline of the space density of available AGN hosts, as a result of the evolution of the specific star-formation rate of the overall galaxy population. Scenarios which attribute the evolution of moderate luminosity AGN since z \approx 1 to changes in the suppermassive black hole accretion mode are not favored by our results.

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

We present a study of the multi-wavelength properties, from the mid-infrared to the hard X-rays, of a sample of 255 spectroscopically identified X-ray selected Type-2 AGN from the XMM-COSMOS survey. Most of them are obscured the X-ray absorbing column density is determined by either X-ray spectral analyses (for the 45% of the sample), or from hardness ratios. Spectral Energy Distributions (SEDs) are computed for all sources in the sample. The average SEDs in the optical band is dominated by the host-galaxy light, especially at low X-ray luminosities and redshifts. There is also a trend between X-ray and mid-infrared luminosity: the AGN contribution in the infrared is higher at higher X-ray luminosities. We calculate bolometric luminosities, bolometric corrections, stellar masses and star formation rates (SFRs) for these sources using a multi-component modeling to properly disentangle the emission associated to stellar light from that due to black hole accretion. For 90% of the sample we also have the morphological classifications obtained with an upgraded version of the Zurich Estimator of Structural Types (ZEST+). We find that on average Type-2 AGN have lower bolometric corrections than Type-1 AGN. Moreover, we confirm that the morphologies of AGN host-galaxies indicate that there is a preference for these Type-2 AGN to be hosted in bulge-dominated galaxies with stellar masses greater than 10^10 solar masses.

We show how accretion rate governs the physical properties of a sample of unobscured broad-line, narrow-line, and lineless active galactic nuclei (AGNs). We avoid the systematic errors plaguing previous studies of AGN accretion rate by using accurate accretion luminosities (L_int) from well-sampled multiwavelength SEDs from the Cosmic Evolution Survey (COSMOS), and accurate black hole masses derived from virial scaling relations (for broad-line AGNs) or host-AGN relations (for narrow-line and lineless AGNs). In general, broad emission lines are present only at the highest accretion rates (L_int/L_Edd > 0.01), and these rapidly accreting AGNs are observed as broad-line AGNs or possibly as obscured narrow-line AGNs. Narrow-line and lineless AGNs at lower specific accretion rates (L_int/L_Edd < 0.01) are unobscured and yet lack a broad line region. The disappearance of the broad emission lines is caused by an expanding radiatively inefficient accretion flow (RIAF) at the inner radius of the accretion disk. The presence of the RIAF also drives L_int/L_Edd < 10^-2 narrow-line and lineless AGNs to 10 times higher ratios of radio to optical/UV emission than L_int/L_Edd > 0.01 broad-line AGNs, since the unbound nature of the RIAF means it is easier to form a radio outflow. The IR torus signature also tends to become weaker or disappear from L_int/L_Edd < 0.01 AGNs, although there may be additional mid-IR synchrotron emission associated with the RIAF. Together these results suggest that specific accretion rate is an important physical "axis" of AGN unification, described by a simple model.

We present Subaru/FOCAS spectropolarimetry of two active galaxies in the Cosmic Evolution Survey. These objects were selected to be optically dull, with the bright X-ray emission of an AGN but missing optical emission lines in our previous spectroscopy. Our new observations show that one target has very weak emission lines consistent with an optically dull AGN, while the other object has strong emission lines typical of a host-diluted Type 2 Seyfert galaxy. In neither source do we observe polarized emission lines, with 3-sigma upper limits of P_BLR < 2%. This means that the missing broad emission lines (and weaker narrow emission lines) are not due to simple anisotropic obscuration, e.g., by the canonical AGN torus. The weak-lined optically dull AGN exhibits a blue polarized continuum with P = 0.78 +/- 0.07% at 4400 A < lambda_rest < 7200 A (P = 1.37 +/- 0.16% at 4400 A < lambda_rest < 5050 A). The wavelength dependence of this polarized flux is similar to that of an unobscured AGN continuum and represents the intrinsic AGN emission, either as synchrotron emission or the outer part of an accretion disk reflected by a clumpy dust scatterer. Because this intrinsic AGN emission lacks emission lines, this source is likely to have a radiatively inefficient accretion flow.

We searched for quasars that are ~ 3 mag fainter than the SDSS quasars in the redshift range 3.7 < z < 4.7 in the COSMOS field to constrain the faint end of the quasar luminosity function. Using optical photometric data, we selected 31 quasar candidates with 22 < i' < 24 at z ~ 4. We obtained optical spectra for most of these candidates using FOCAS on the Subaru telescope, and identified 8 low-luminosity quasars at z ~ 4. In order to derive the quasar luminosity function (QLF) based on our spectroscopic follow-up campaign, we estimated the photometric completeness of our quasar survey through detailed Monte Carlo simulations. Our QLF at z ~ 4 has a much shallower faint-end slope beta = -1.67^{+0.11}_{-0.17} than that obtained by other recent surveys in the same redshift. Our result is consistent with the scenario of downsizing evolution of active galactic nuclei inferred by recent optical and X-ray quasar surveys at lower redshifts.

What is the relevance of major mergers and interactions as triggering mechanisms for active galactic nuclei (AGN) activity? To answer this longstanding question, we analyze 140 XMM-selected AGN host galaxies and a matched control sample of 1264 inactive galaxies over z~0.3-1.0 and log(M_*/M_sun)<11.7 with high-resolution HST/ACS imaging from the COSMOS field. The visual analysis of their morphologies by 10 independent human classifiers yields a measure of the fraction of distorted morphologies in the AGN and control samples, i.e. quantifying the signature of recent mergers which might potentially be responsible for fueling/triggering the AGN. We find that (1) the vast majority (>85%) of the AGN host galaxies do not show strong distortions, and (2) there is no significant difference in the distortion fractions between active and inactive galaxies. Our findings provide the best direct evidence that, since z~1, the bulk of black hole accretion has not been triggered by major galaxy mergers, therefore arguing that the alternative mechanisms, i.e., secular processes and minor interactions, are the leading triggers for the episodes of major black hole growth. We also exclude an alternative interpretation of our results: a significant time lag between merging and the observability of the AGN phase could wash out the most significant merging signatures, explaining the lack of enhancement of strong distortions on the AGN hosts. We show that this alternative scenario is unlikely due to: (1) recent major mergers being ruled out for the majority of sources due to the high fraction of disk-hosted AGN, (2) the lack of a significant X-ray signal in merging inactive galaxies as a signature of a potential buried AGN, and (3) the low levels of soft X-ray obscuration for AGN hosted by interacting galaxies, in contrast to model predictions.

We report a sizable class of type 1 active galactic nuclei (AGNs) with unusually weak near-infrared (1-3{\mu}m) emission in the XMM-COSMOS type 1 AGN sample. The fraction of these "hot-dust-poor" AGNs increases with redshift from 6% at lowredshift (z < 2) to 20% at moderate high redshift (2 < z < 3.5). There is no clear trend of the fraction with other parameters: bolometric luminosity, Eddington ratio, black hole mass, and X-ray luminosity. The 3{\mu}m emission relative to the 1{\mu}m emission is a factor of 2-4 smaller than the typical Elvis et al. AGN spectral energy distribution (SED), which indicates a "torus" covering factor of 2%-29%, a factor of 3-40 smaller than required by unified models. The weak hot dust emission seems to expose an extension of the accretion disk continuum in some of the source SEDs. We estimate the outer edge of their accretion disks to lie at $(0.3-2.0) /times 10^4$ Schwarzschild radii, ~10-23 times the gravitational stability radii. Formation scenarios for these sources are discussed.

Based on the recent very deep near-infrared imaging of the Hubble Ultra Deep Field with WFC3 on the Hubble Space Telescope, five groups published most probable samples of galaxies at z~8, selected by the so-called dropout method or photometric redshift; e.g., Y_105-dropouts (Y_105-J_125 > 0.8). These studies are highly useful for investigating both the early star formation history of galaxies and the sources of cosmic re-ionization. In order to better understand these issues, we carefully examine if there are low-$z$ interlopers in the samples of z~8 galaxy candidates. We focus on the strong emission-line galaxies at z~2 in this paper. Such galaxies may be selected as Y_105-dropouts since the [OIII] lambda 5007 emission line is redshifted into the J_125-band. We have found that the contamination from such low-$z$ interlopers is negligibly small. Therefore, all objects found by the five groups are free from this type of contamination. However, it remains difficult to extract real z~8 galaxies because all the sources are very faint and the different groups have found different candidates. With this in mind, we construct a robust sample of eight galaxies at z~8 from the objects found by the five groups: each of these eight objects has been selected by at least two groups. Using this sample, we discuss their UV continuum slope. We also discuss the escape fraction of ionizing photons adopting various metallicities. Our analysis suggests that massive stars forming in low-metallicity gas (Z~5 \times 10^-4 Z_sun) can be responsible for the completion of cosmic re-ionization if the escape fraction of ionizing continuum from galaxies is as large as 0.5, and this is consistent with the observed blue UV continua.

We analyze the morphological properties of a large sample of 1503 70 micron selected galaxies in the COSMOS field spanning the redshift range 0.01<z< 3.5 with a median redshift of 0.5 and an infrared luminosity range of 10^8<L_IR<10^14L_sun with a median luminosity of 10^11.4 L_sun. In general these galaxies are massive, with a stellar mass range of 10^10-10^12 M_sun, and luminous, with -25<M_K<-20. We find a strong correlation between the fraction of major mergers and L_IR, with the fraction at the highest luminosity being up to 50%. We also find that the fraction of spirals drops dramatically with L_IR. Minor mergers likely play a role in boosting the infrared luminosity for sources with low luminosities. The precise fraction of mergers in any given L_IR bin varies by redshift due to sources at z>1 being difficult to classify and subject to the effects of band pass shifting, therefore, these numbers can only be considered lower limits. At z<1, where the morphological classifications are most robust, major mergers clearly dominate the ULIRG population (50-80%) and are important for the LIRG population (25-40%). At z>1 the fraction of major mergers is at least 30-40% for ULIRGs. Although the general morphological trends agree with what has been observed for local (U)LIRGs, the fraction of major mergers is slightly lower than seen locally. This is in part due to the difficulty of identifying merger signatures at high redshift. We argue that given the number of major gas-rich mergers observed and the relatively short timescale that they would be observable in the (U)LIRG phase that it is plausible for the observed red sequence of massive ellipticals (<10^12 M_sun) to have been formed entirely by gas-rich major mergers.

We report the final optical identifications of the medium-depth (~60 ksec), contiguous (2 deg^2) XMM-Newton survey of the COSMOS field. XMM-Newton has detected ~800 X-ray sources down to limiting fluxes of ~5x10^{-16}, ~3x10^{-15}, and ~7x10^{-15} erg/cm2/s in the 0.5-2 keV, 2-10 keV and 5-10 keV bands, respectively. The work is complemented by an extensive collection of multi-wavelength data from 24 micron to UV, available from the COSMOS survey, for each of the X-ray sources, including spectroscopic redshifts for ~50% of the sample, and high-quality photometric redshifts for the rest. The XMM and multiwavelength flux limits are well matched: 1760 (98%) of the X-ray sources have optical counterparts, 1711 (~95%) have IRAC counterparts, and 1394 (~78%) have MIPS 24micron detections. Thanks to the redshift completeness (almost 100%) we were able to constrain the high-luminosity tail of the X-ray luminosity function confirming that the peak of the number density of logL_X>44.5 AGN is at z~2. Spectroscopically-identified obscured and unobscured AGN, as well as normal and starforming galaxies, present well-defined optical and infrared properties. We devised a robust method to identify a sample of ~150 high redshift (z>1), obscured AGN candidates for which optical spectroscopy is not available. We were able to determine that the fraction of the obscured AGN population at the highest (L_X>10^{44} erg s^{-1}) X-ray luminosity is ~15-30% when selection effects are taken into account, providing an important observational constraint for X-ray background synthesis. We studied in detail the optical spectrum and the overall spectral energy distribution of a prototypical Type 2 QSO, caught in a stage transitioning from being starburst dominated to AGN dominated, which was possible to isolate only thanks to the combination of X-ray and infrared observations.

We present a detailed study of a peculiar source in the COSMOS survey at z=0.359. Source CXOCJ100043.1+020637 (CID-42) presents two compact optical sources embedded in the same galaxy. The distance between the 2, measured in the HST/ACS image, is 0.495" that, at the redshift of the source, corresponds to a projected separation of 2.46 kpc. A large (~1200 km/s) velocity offset between the narrow and broad components of Hbeta has been measured in three different optical spectra from the VLT/VIMOS and Magellan/IMACS instruments. CID-42 is also the only X-ray source having in its X-ray spectra a strong redshifted broad absorption iron line, and an iron emission line, drawing an inverted P-Cygni profile. The Chandra and XMM data show that the absorption line is variable in energy by 500 eV over 4 years and that the absorber has to be highly ionized, in order not to leave a signature in the soft X-ray spectrum. That these features occur in the same source is unlikely to be a coincidence. We envisage two possible explanations: (1) a gravitational wave recoiling black hole (BH), caught 1-10 Myr after merging, (2) a Type 1/ Type 2 system in the same galaxy where the Type 1 is recoiling due to slingshot effect produced by a triple BH system. The first possibility gives us a candidate gravitational waves recoiling BH with both spectroscopic and imaging signatures. In the second case, the X-ray absorption line can be explained as a BAL-like outflow from the foreground nucleus (a Type 2 AGN) at the rearer one (a Type 1 AGN), which illuminates the otherwise undetectable wind, giving us the first opportunity to show that fast winds are present in obscured AGN.

Using Chandra observations we have identified a sample of seven off-nuclear X-ray sources, in the redshift range z=0.072-0.283, located within optically bright galaxies in the COSMOS Survey. Using the multi-wavelength coverage available in the COSMOS field, we study the properties of the host galaxies of these ULXs. In detail, we derived their star formation rate from H_alpha measurements and their stellar masses using SED fitting techniques with the aim to compute the probability to have an off-nuclear source based on the host galaxy properties. We divide the host galaxies in different morphological classes using the available ACS/HST imaging. We find that our ULXs candidates are located in regions of the SFR versus M$_star$ plane where one or more off-nuclear detectable sources are expected. From a morphological analysis of the ACS imaging and the use of rest-frame colours, we find that our ULXs are hosted both in late and early type galaxies. Finally, we find that the fraction of galaxies hosting a ULX ranges from ~0.5% to ~0.2% going from L[0.5-2 keV]=3 x 10^39 erg s^-1 to L[0.5-2 keV]= 2 x 10^40 erg s^-1.

We present a study of the X-ray to optical properties of a sample of 545 X-ray selected Type 1 AGN, from the XMM-COSMOS survey, over a wide range of redshifts ($0.04<\z<4.25$) and X-ray luminosities ($40.6 \leq \Log \Lhard \leq 45.3$). About 60% of them are spectroscopically identified Type 1 AGN, while the others have a reliable photometric redshift and are classified as Type 1 AGN on the basis of their multi-band Spectral Energy Distributions. We discuss the relationship between UV and X-ray luminosity, as parameterized by the $\alphaox$ spectral slope, and its dependence on redshift and luminosity. We compare our findings with previous investigations of optically selected broad-line AGN (mostly from SDSS). A highly significant correlation between $\alphaox$ and $\lo$ is found, in agreement with previous investigations of optically selected samples. We calculate bolometric corrections, $\kbol$, for the whole sample using hard X-ray luminosities ($\Lhard$), and the Eddington ratios for a subsample of 150 objects for which black hole mass estimates are available. We confirm the trend of increasing bolometric correction with increasing Eddington ratio as proposed in previous works. A tight correlation is found between $\alphaox$ and $\kbol$, which can be used to estimate accurate bolometric corrections using only optical and X-ray data. We find a significant correlation between $\alphaox$ and Eddington ratio, in which $\alphaox$ increases for increasing Eddington ratios.

We present a large robust sample of 1503 reliable and unconfused 70microm selected sources from the multiwavelength data set of the Cosmic Evolution Survey (COSMOS). Using the Spitzer IRAC and MIPS photometry, we estimate the total infrared luminosity, L_IR (8--1000 microns), by finding the best fit template from several different template libraries. The long wavelength 70 and 160 micron data allow us to obtain a reliable estimate of L_IR, accurate to within 0.2 and 0.05 dex, respectively. The 70 micron data point enables a significant improvement over the luminosity estimates possible with only a 24 micron detection. The full sample spans a wide range in L_IR, L_IR ~ 10^8-10^14 L_sun, with a median luminosity of 10^11.4 L_sun. We identify a total of 687 luminous, 303 ultraluminous, and 31 hyperluminous infrared galaxies (LIRGs, ULIRGs, and HyLIRGs) over the redshift range 0.01<z<3.5 with a median redshift of 0.5. Presented here are the full spectral energy distributions for each of the sources compiled from the extensive multiwavelength data set from the ultraviolet (UV) to the far-infrared (FIR). Using SED fits we find possible evidence for a subset of cooler ultraluminous objects than observed locally. However, until direct observations at longer wavelengths are obtained, the peak of emission and the dust temperature cannot be well constrained. We use these SEDs, along with the deep radio and X-ray coverage of the field, to identify a large sample of candidate active galactic nuclei (AGN). We find that the fraction of AGN increases strongly with L_IR, as it does in the local universe, and that nearly 70% of ULIRGs and all HyLIRGs likely host a powerful AGN.

We constrain the ratio of black hole (BH) mass to total stellar mass of type-1 AGN in the COSMOS survey at 1<z<2. For 10 AGN at mean redshift z~1.4 with both HST/ACS and HST/NICMOS imaging data we are able to compute total stellar mass M_(*,total), based on restframe UV-to-optical host galaxy colors which constrain mass-to-light ratios. All objects have virial BH mass-estimates available from the COSMOS Magellan/IMACS and zCOSMOS surveys. We find zero difference between the M_BH--M_(*,total)-relation at z~1.4 and the M_BH--M_(*,bulge)-relation in the local Universe. Our interpretation is: (a) If our objects were purely bulge-dominated, the M_BH--M_(*,bulge)-relation has not evolved since z~1.4. However, (b) since we have evidence for substantial disk components, the bulges of massive galaxies (logM_(*,total)=11.1+-0.25 or logM_BH~8.3+-0.2) must have grown over the last 9 Gyrs predominantly by redistribution of disk- into bulge-mass. Since all necessary stellar mass exists in the galaxy at z=1.4, no star-formation or addition of external stellar material is required, only a redistribution e.g. induced by minor and major merging or through disk instabilities. Merging, in addition to redistributing mass in the galaxy, will add both BH and stellar/bulge mass, but does not change the overall final M_BH/M_(*,bulge) ratio. Since the overall cosmic stellar and BH mass buildup trace each other tightly over time, our scenario of bulge-formation in massive galaxies is independent of any strong BH-feedback and means that the mechanism coupling BH and bulge mass until the present is very indirect.

(Abriged) We report on the measurement of the rest frame K-band luminosity and total stellar mass of the hosts of 89 broad line Active Galactic Nuclei detected in the zCOSMOS survey in the redshift range 1<z<2.2. The unprecedented multiwavelength coverage of the survey field allows us to disentangle the emission of the host galaxy from that of the nuclear black hole in their Spectral Energy Distributions. We derive an estimate of black hole masses through the analysis of the broad Mg II emission lines observed in the medium-resolution spectra taken with VIMOS/VLT as part of the zCOSMOS project. We found that, as compared to the local value, the average black hole to host galaxy mass ratio appears to evolve positively with redshift, with a best fit evolution of the form (1+z)^{0.68 \pm0.12 +0.6 -0.3}, where the large asymmetric systematic errors stem from the uncertainties in the choice of IMF, in the calibration of the virial relation used to estimate BH masses and in the mean QSO SED adopted. A thorough analysis of observational biases induced by intrinsic scatter in the scaling relations reinforces the conclusion that an evolution of the MBH-M* relation must ensue for actively growing black holes at early times: either its overall normalization, or its intrinsic scatter (or both) appear to increase with redshift. This can be interpreted as signature of either a more rapid growth of supermassive black holes at high redshift, a change of structural properties of AGN hosts at earlier times, or a significant mismatch between the typical growth times of nuclear black holes and host galaxies.

We present infrared, optical, and X-ray data of 48 X-ray bright, optically dull AGNs in the COSMOS field. These objects exhibit the X-ray luminosity of an active galactic nucleus (AGN) but lack broad and narrow emission lines in their optical spectrum. We show that despite the lack of optical emission lines, most of these optically dull AGNs are not well-described by a typical passive red galaxy spectrum: instead they exhibit weak but significant blue emission like an unobscured AGN. Photometric observations over several years additionally show significant variability in the blue emission of four optically dull AGNs. The nature of the blue and infrared emission suggest that the optically inactive appearance of these AGNs cannot be caused by obscuration intrinsic to the AGNs. Instead, up to ~70% of optically dull AGNs are diluted by their hosts, with bright or simply edge-on hosts lying preferentially within the spectroscopic aperture. The remaining ~30% of optically dull AGNs have anomalously high f_x/f_o ratios and are intrinsically weak, not obscured, in the optical. These optically dull AGNs are best described as a weakly accreting AGN with a truncated accretion disk from a radiatively inefficient accretion flow.

We present detailed morphological properties of Lyman alpha emitters (LAEs) at z~ 5.7 in the COSMOS field, based on {\it Hubble Space Telescope} Advanced Camera for Surveys (ACS) data. The ACS imaging in the F814W filter covered 85 LAEs of the 119 LAEs identified in the full two square degree field, and 47 LAEs of them are detected in the ACS images. Nearly half of them are spatially extended with a size larger than 0.15 arcsec (~0.88 kpc at z=5.7) up to 0.4 arcsec (~2.5 kpc at z=5.7). The others are nearly unresolved compact objects. Two LAEs show double-component structures, indicating interaction or merging of building components to form more massive galaxies. By stacking the ACS images of all the detected sources, we obtain a Sersic parameter of n~0.7 with a half-light radius of 0.13 arcsec (0.76 kpc), suggesting that the majority of ACS detected LAEs have not spheroidal-like but disk-like or irregular light profiles. Comparing ACS F814W magnitudes (I_814) with Subaru/Suprime-Cam magnitudes in the NB816, i', and z' bands, we find that the ACS imaging in the F814W band mainly probes UV continuum rather than Lyman alpha line emission. UV continuum sizes tend to be larger for LAEs with larger Ly\alpha emission regions as traced by the NB816 imaging. The non-detection of 38 LAEs in the ACS images is likely due to the fact that their surface brightness is even too low both in the UV continuum and Lyalpha emission. Estimating I_814 for the LAEs with ACS non-detection from the z' and NB816 magnitudes, we find that 16 of these are probably LAEs with a size larger than 0.15 arcsec in UV continuum. All these results suggest that our LAE sample contains systematically larger LAEs in UV continuum size than those previously studied at z~6.

We investigate the relation between the star-formation activity in galaxies and environment at z=1.2 in the COSMOS field, using the fraction of [OII] emitters and the local galaxy density. The fraction of [OII] emitters appears to be almost constant over the surface density of galaxies between 0.2 and 10 Mpc^-2. This trend is different from that seen in the local universe where the star-formation activity is weaker in higher density regions. To understand this difference between z~1 and z~0, we study the fraction of non-isolated galaxies as a function of local galaxy density. We find that the fraction of non-isolated galaxies increases with increasing density. Our results suggest that the star formation in galaxies at z~1 is triggered by galaxy interaction and/or mergers.

We present black hole masses and accretion rates for 182 Type 1 AGN in COSMOS. We estimate masses using the scaling relations for the broad Hb, MgII, and CIV emission lines in the redshift ranges 0.16<z<0.88, 1<z<2.4, and 2.7<z<4.9. We estimate the accretion rate using an Eddington ratio L_I/L_Edd estimated from optical and X-ray data. We find that very few Type 1 AGN accrete below L_I/L_Edd ~ 0.01, despite simulations of synthetic spectra which show that the survey is sensitive to such Type 1 AGN. At lower accretion rates the BLR may become obscured, diluted or nonexistent. We find evidence that Type 1 AGN at higher accretion rates have higher optical luminosities, as more of their emission comes from the cool (optical) accretion disk with respect to shorter wavelengths. We measure a larger range in accretion rate than previous works, suggesting that COSMOS is more efficient at finding low accretion rate Type 1 AGN. However the measured range in accretion rate is still comparable to the intrinsic scatter from the scaling relations, suggesting that Type 1 AGN accrete at a narrow range of Eddington ratio, with L_I/L_Edd ~ 0.1.

We present results of a survey for Ly alpha emitters at z=4.86 based on optical narrowband (lambda_c=7126 angstrom, Delta lambda=73 angstrom) and broadband (B, V, r', i', and z') observations of the Cosmic Evolution Survey (COSMOS) field using Suprime-Cam on the Subaru Telescope. We find 79 LAE candidates at z=4.86 over a contiguous survey area of 1.83 deg^2, down to the Ly alpha line flux of 1.47 x 10^-17 ergs s^-1 cm^-2. We obtain the Ly alpha luminosity function with a best-fit Schechter parameters of log L^*=42.9^+0.5_-0.3 ergs s^-1 and phi^* = 1.2^+8.0_-1.1 x 10^-4 Mpc^-3 for alpha=-1.5 (fixed). The two-point correlation function for our LAE sample is xi(r) = (r/4.4^+5.7_-2.9 Mpc)^-1.90+/-0.22. In order to investigate the field-to-field variations of the properties of Ly alpha emitters, we divide the survey area into nine tiles of 0.5^circ x 0.5^circ each. We find that the number density varies with a factor of ~ 2 from field to field with high statistical significance. However, we find no significant field-to-field variance when we divide the field into four tiles with 0.7^circ x 0.7^circ each. We conclude that at least 0.5 deg^2 survey area is required to derive averaged properties of LAEs at z~5, and our survey field is wide enough to overcome the cosmic variance.

We study the spatial clustering of 538 X-ray selected AGN in the 2 deg^2 XMM-COSMOS field that are spectroscopically identified to I_{AB}<23 and span the redshift range z=0.2-3.0. The median redshift and luminosity of the sample are z = 0.98 and L_{0.5-10}=6.3 x 10^{43} erg/s, respectively. A strong clustering signal is detected at ~18sigma level, which is the most significant measurement obtained to date for clustering of X-ray selected AGN. By fitting the projected correlation function w(r_p) with a power law on scales of r_p=0.3-40 Mpc/h, we derive a best fit comoving correlation length of r_0 = 8.6 +- 0.5 Mpc/h and slope of gamma=1.88 +- 0.07 (Poissonian errors; bootstrap errors are about a factor of 2 larger). An excess signal is observed in the range r_p~5-15 Mpc/h, which is due to a large scale structure at z ~ 0.36 containing about 40 AGN. When removing the z ~ 0.36 structure, or computing w(r_p) in a narrower range around the peak of the redshift distribution (e.g. z=0.4-1.6), the correlation length decreases to r_0 ~ 5-6 Mpc/h, which is consistent with that observed for bright optical QSOs at the same redshift. We investigate the clustering properties of obscured and unobscured AGN separately. Within the statistical uncertainties, we do not find evidence that AGN with broad optical lines (BLAGN) cluster differently from AGN without broad optical lines (non-BLAGN). The correlation length measured for XMM-COSMOS AGN at z~1 is similar to that of massive galaxies (stellar mass M_*> 3 x 10^{10} M_sun) at the same redshift. This suggests that AGN at z~1 are preferentially hosted by massive galaxies, as observed both in the local and in the distant (z~2) Universe. (shortened)

We present optical spectroscopy for an X-ray and optical flux-limited sample of 677 XMM-Newton selected targets covering the 2 deg^2 COSMOS field, with a yield of 485 high-confidence redshifts. The majority of the spectra were obtained over three seasons (2005-2007) with the IMACS instrument on the Magellan (Baade) telescope. We also include in the sample previously published Sloan Digital Sky Survey spectra and supplemental observations with MMT/Hectospec. We detail the observations and classification analyses. The survey is 90% complete to flux limits of f_{0.5-10 keV}>8 x 10^-16 erg cm^-2 s^-1 and i_AB+<22, where over 90% of targets have high-confidence redshifts. Making simple corrections for incompleteness due to redshift and spectral type allows for a description of the complete population to $i_AB+<23. The corrected sample includes 57% broad emission line (Type 1, unobscured) AGN at 0.13<z<4.26, 25% narrow emission line (Type 2, obscured) AGN at 0.07<z<1.29, and 18% absorption line (host-dominated, obscured) AGN at 0<z<1.22 (excluding the stars that made up 4% of the X-ray targets). We show that the survey's limits in X-ray and optical flux include nearly all X-ray AGN (defined by L_{0.5-10 keV}>3 x 10^42 erg s^-1) to z<1, of both optically obscured and unobscured types. We find statistically significant evidence that the obscured to unobscured AGN ratio at z<1 increases with redshift and decreases with luminosity.

(abridged) We take advantage of the deep Chandra and Spitzer coverage of a large area (more than 10 times the area covered by the Chandra deep fields, CDFs in the COSMOS field, to extend the search of highly obscured, Compton-thick active nuclei to higher luminosity. These sources have low surface density and large samples can be provided only through large area surveys, like the COSMOS survey. We analyze the X-ray properties of COSMOS MIPS sources with 24$\mu$m fluxes higher than 550$\mu$Jy. For the MIPS sources not directly detected in the Chandra images we produce stacked images in soft and hard X-rays bands. To estimate the fraction of Compton-thick AGN in the MIPS source population we compare the observed stacked count rates and hardness ratios to those predicted by detailed Monte Carlo simulations including both obscured AGN and star-forming galaxies. The density of lower luminosity Compton-thick AGN (logL(2-10keV)=43.5-44) at z=0.7--1.2 is $(3.7\pm1.1) \times10^{-5}$ Mpc$^{-3}$, corresponding to $\sim67%$ of that of X-ray selected AGN. The comparison between the fraction of infrared selected, Compton thick AGN to the X-ray selected, unobscured and moderately obscured AGN at high and low luminosity suggests that Compton-thick AGN follow a luminosity dependence similar to that discovered for Compton-thin AGN, becoming relatively rarer at high luminosities. We estimate that the fraction of AGN (unobscured, moderately obscured and Compton thick) to the total MIPS source population is $49\pm10%$, a value significantly higher than that previously estimated at similar 24$\mu$m fluxes. We discuss how our findings can constrain AGN feedback models.

We present a new measurement of the space density of high redshift (3.0<z<4.5), X-ray selected QSOs obtained by exploiting the deep and uniform multiwavelength coverage of the COSMOS survey. We have assembled a statistically large (40 objects), X-ray selected (F_{0.5-2 keV} >10^{-15} cgs), homogeneous sample of z>3 QSOs for which spectroscopic (22) or photometric (18) redshifts are available. We present the optical (color-color diagrams) and X-ray properties, the number counts and space densities of the z>3 X-ray selected quasars population and compare our findings with previous works and model predictions. We find that the optical properties of X-ray selected quasars are not significantly different from those of optically selected samples. There is evidence for substantial X-ray absorption (logN_H>23 cm^{-2}) in about 20% of the sources in the sample. The comoving space density of luminous (L_X >10^{44} erg s^-1) QSOs declines exponentially (by an e--folding per unit redshift) in the z=3.0-4.5 range, with a behavior similar to that observed for optically bright unobscured QSOs selected in large area optical surveys. Prospects for future, large and deep X-ray surveys are also discussed.

We present photometric redshifts and spectral energy distribution (SED) classifications for a sample of 1542 optically identified sources detected with XMM in the COSMOS field. Our template fitting classifies 46 sources as stars and 464 as non-active galaxies, while the remaining 1032 require templates with an AGN contribution. High accuracy in the derived photometric redshifts was accomplished as the result of 1) photometry in up to 30 bands with high significance detections, 2) a new set of SED templates including 18 hybrids covering the far-UV to mid-infrared, which have been constructed by the combination of AGN and non-active galaxies templates, and 3) multi-epoch observations that have been used to correct for variability (most important for type 1 AGN). The reliability of the photometric redshifts is evaluated using the sub-sample of 442 sources with measured spectroscopic redshifts. We achieved an accuracy of $\sigma_{\Delta z/(1+z_{spec})} = 0.014$ for i$_{AB}^*<$22.5 ($\sigma_{\Delta z/(1+z_{spec})} \sim0.015$ for i$_{AB}^*<$24.5). The high accuracies were accomplished for both type 2 (where the SED is often dominated by the host galaxy) and type 1 AGN and QSOs out to $z=4.5$. The number of outliers is a large improvement over previous photometric redshift estimates for X-ray selected sources (4.0% and 4.8% outliers for i$_{AB}^*<$22.5 and i$_{AB}^*<$24.5, respectively). We show that the intermediate band photometry is vital to achieving accurate photometric redshifts for AGN, whereas the broad SED coverage provided by mid infrared (Spitzer/IRAC) bands is important to reduce the number of outliers for normal galaxies.

We use HST/ACS images and a photometric catalog of the COSMOS field to analyze morphologies of the host galaxies of approximately 400 AGN candidates at redshifts 0.3 < z < 1.0. We compare the AGN hosts with a sample of non-active galaxies drawn from the COSMOS field to match the magnitude and redshift distribution of the AGN hosts. We perform 2-D surface brightness modeling with GALFIT to yield host galaxy and nuclear point source magnitudes. X-ray selected AGN host galaxy morphologies span a substantial range that peaks between those of early-type, bulge-dominated and late-type, disk-dominated systems. We also measure the asymmetry and concentration of the host galaxies. Unaccounted for, the nuclear point source can significantly bias results of these measured structural parameters, so we subtract the best-fit point source component to obtain images of the underlying host galaxies. Our concentration measurements reinforce the findings of our 2-D morphology fits, placing X-ray AGN hosts between early- and late-type inactive galaxies. AGN host asymmetry distributions are consistent with those of control galaxies. Combined with a lack of excess companion galaxies around AGN, the asymmetry distributions indicate that strong interactions are no more prevalent among AGN than normal galaxies. In light of recent work, these results suggest that the host galaxies of AGN at these X-ray luminosities may be in a transition from disk-dominated to bulge-dominated, but that this transition is not typically triggered by major mergers.

We explore the properties of the submillijansky radio population at 20 cm by applying a newly developed optical color-based method to separate star forming (SF) from AGN galaxies at intermediate redshifts (z<1.3). Although optical rest-frame colors are used, our separation method is shown to be efficient, and not biased against dusty starburst galaxies. This classification method has been calibrated and tested on a local radio selected optical sample. Given accurate multi-band photometry and redshifts, it carries the potential to be generally applicable to any galaxy sample where SF and AGN galaxies are the two dominant populations. In order to quantify the properties of the submillijansky radio population, we have analyzed ~2,400 radio sources, detected at 20 cm in the VLA-COSMOS survey. 90% of these have submillijansky flux densities. We classify the objects into 1) star candidates, 2) quasi stellar objects, 3) AGN, 4) SF, and 5) high redshift (z>1.3) galaxies. We find, for the composition of the submillijansky radio population, that SF galaxies are not the dominant population at submillijansky flux levels, as previously often assumed, but that they make up an approximately constant fraction of 30-40% in the flux density range of ~50 microJy to 0.7 mJy. In summary, based on the entire VLA-COSMOS radio population at 20 cm, we find that the radio population at these flux densities is a mixture of roughly 30-40% of SF and 50-60% of AGN galaxies, with a minor contribution (~10%) of QSOs.

We outline a strategy to select faint (i<24.5) type 1 AGN candidates down to the Seyfert/QSO boundary for spectroscopic targeting in the COSMOS field, picking candidates by their nonstellar colors in broadband ground-based photometry and morphological properties extracted from HST-ACS. AGN optical color selection has not been applied to such faint magnitudes in such a large continuous part of the sky. Hot stars are known to be the dominant contaminant for bright AGN candidate selection at z<2, but we anticipate the highest color contamination at all redshifts to be from faint starburst and compact galaxies. Morphological selection via the Gini Coefficient separates most potential AGN from these faint blue galaxies. Recent models of the quasar luminosity function are used to estimate quasar surface densities, and studies of stellar populations in the COSMOS field infer stellar contamination. We use 292 spectroscopically confirmed type 1 AGN and quasar templates to predict AGN colors with redshift, and contrast those predictions with the colors of known contaminating populations. The motivation of this study and subsequent spectroscopic follow-up is to populate and refine the faint end of the QLF where the population of type 1 AGN is presently not well known. The anticipated AGN observations will add to the ~300 already known AGN in the COSMOS field, making COSMOS a densely packed field of quasars to be used to understand supermassive black holes and probe the structure of the intergalactic medium in the intervening volume.

In this work we use a sample of 318 radio-quiet quasars (RQQ) to investigate the dependence of the ratio of optical/UV flux to X-ray flux, alpha_ox, and the X-ray photon index, Gamma_X, on black hole mass, UV luminosity relative to Eddington, and X-ray luminosity relative to Eddington. Our sample is drawn from the SDSS, with X-ray data from ROSAT and Chandra, and optical data mostly from the SDSS; 153 of these sources have estimates of Gamma_X from Chandra. We estimate M_BH using standard estimates derived from the Hbeta, Mg II, and C IV broad emission lines. Our sample spans a broad range in black hole mass (10^6 < M_BH / M_Sun < 10^10) and redshift (z < 4.8). We find that alpha_ox increases with increasing M_BH and L_UV / L_Edd, and decreases with increasing L_X / L_Edd. In addition, we confirm the correlation seen in previous studies between Gamma_X and M_BH and both L_UV / L_Edd and L_X / L_Edd; however, we also find evidence that the dependence of Gamma_X of these quantities is not monotonic, changing sign at M_BH ~ 3 x 10^8 M_Sun. We argue that the alpha_ox correlations imply that the fraction of bolometric luminosity emitted by the accretion disk, as compared to the corona, increases with increasing accretion rate relative to Eddington. In addition, we argue that the Gamma_X trends are caused by a dependence of X-ray spectral index on accretion rate. We discuss our results within the context of accretion models with comptonizing corona, and discuss the implications of the alpha_ox correlations for quasar feedback. To date, this is the largest study of the dependence of RQQ X-ray parameters on black hole mass and related quantities, and the first to attempt to correct for the large statistical uncertainty in the broad line mass estimates.

We present a detailed spectral analysis of point-like X-ray sources in the XMM-COSMOS field. Our sample of 135 sources only includes those that have more than 100 net counts in the 0.3-10 keV energy band and have been identified through optical spectroscopy. The majority of the sources are well described by a simple power-law model with either no absorption (76%) or a significant intrinsic, absorbing column (20%).As expected, the distribution of intrinsic absorbing column densities is markedly different between AGN with or without broad optical emission lines. We find within our sample four Type-2 QSOs candidates (L_X > 10^44 erg/s, N_H > 10^22 cm^-2), with a spectral energy distribution well reproduced by a composite Seyfert-2 spectrum, that demonstrates the strength of the wide field XMM/COSMOS survey to detect these rare and underrepresented sources.

[ABRIGED] We present the optical identification of a sample of 695 X-ray sources detected in the first 1.3 deg^2 of the XMM-COSMOS survey, down to a 0.5-2 keV (2-10 keV) limiting flux of ~10^-15 erg cm-2 s-1 (~5x10^-15 erg cm^-2 s-1). We were able to associate a candidate optical counterpart to ~90% (626) of the X-ray sources, while for the remaining ~10% of the sources we were not able to provide a unique optical association due to the faintness of the possible optical counterparts (I_AB>25) or to the presence of multiple optical sources within the XMM-Newton error circles. We also cross-correlated the candidate optical counterparts with the Subaru multicolor and ACS catalogs and with the Magellan/IMACS, zCOSMOS and literature spectroscopic data; the spectroscopic sample comprises 248 objects (~40% of the full sample). Our analysis reveals that for ~80% of the counterparts there is a very good agreement between the spectroscopic classification, the morphological parameters as derived from ACS data, and the optical to near infrared colors. About 20% of the sources show an apparent mismatch between the morphological and spectroscopic classifications. All the ``extended'' BL AGN lie at redshift <1.5, while the redshift distribution of the full BL AGN population peaks at z~1.5. Our analysis also suggests that the Type 2/Type 1 ratio decreases towards high luminosities, in qualitative agreement with the results from X-ray spectral analysis and the most recent modeling of the X-ray luminosity function evolution.

We present the results of a search for galaxy clusters in the first 36 XMM-Newton pointings on the COSMOS field. We reach a depth for a total cluster flux in the 0.5-2 keV band of 3x10-15 ergs cm-2 s-1, having one of the widest XMM-Newton contiguous raster surveys, covering an area of 2.1 square degrees. Cluster candidates are identified through a wavelet detection of extended X-ray emission. Verification of the cluster candidates is done based on a galaxy concentration analysis in redshift slices of thickness of 0.1-0.2 in redshift, using the multi-band photometric catalog of the COSMOS field and restricting the search to z<1.3 and i_AB < 25. We identify 72 clusters and derive their properties based on the X-ray cluster scaling relations. A statistical description of the survey in terms of the cumulative log(N>S)-lg(S) distribution compares well with previous results, although yielding a somewhat higher number of clusters at similar fluxes. The X-ray luminosity function of COSMOS clusters matches well the results of nearby surveys, providing a comparably tight constraint on the faint end slope of alpha=1.93+/-0.04. For the probed luminosity range of 8x10+42 - 2x10+44 ergs s-1, our survey is in agreement with and adds significantly to the existing data on the cluster luminosity function at high redshifts and implies no substantial evolution at these luminosities to z=1.3.

We have identified a complex galaxy cluster system in the COSMOS field via a wide angle tail (WAT) radio galaxy consistent with the idea that WAT galaxies can be used as tracers of clusters. The WAT galaxy, CWAT-01, is coincident with an elliptical galaxy resolved in the HST-ACS image. Using the COSMOS multiwavelength data set, we derive the radio properties of CWAT-01 and use the optical and X-ray data to investigate its host environment. The cluster hosting CWAT-01 is part of a larger assembly consisting of a minimum of four X-ray luminous clusters within ~2 Mpc distance. We apply hydrodynamical models that combine ram pressure and buoyancy forces on CWAT-01. These models explain the shape of the radio jets only if the galaxy's velocity relative to the intra-cluster medium (ICM) is in the range of about 300-550 km/s which is higher than expected for brightest cluster galaxies (BCGs) in relaxed systems. This indicates that the CWAT-01 host cluster is not relaxed, but is possibly dynamically young. We argue that such a velocity could have been induced through subcluster merger within the CWAT-01 parent cluster and/or cluster-cluster interactions. Our results strongly indicate that we are witnessing the formation of a large cluster from an assembly of multiple clusters, consistent with the hierarchical scenario of structure formation. We estimate the total mass of the final cluster to be approximately 20% of the mass of the Coma cluster.

We present spectroscopic redshifts for the first 466 X-ray and radio-selected AGN targets in the 2 deg^2 COSMOS field. Spectra were obtained with the IMACS instrument on the Magellan (Baade) telescope, using the nod-and-shuffle technique. We identify a variety of Type 1 and Type 2 AGN, as well as red galaxies with no emission lines. Our redshift yield is 72% down to i_AB=24, although the yield is >90% for i_AB<22. We expect the completeness to increase as the survey continues. When our survey is complete and additional redshifts from the zCOSMOS project are included, we anticipate ~1100 AGN with redshifts over the entire COSMOS field. Our redshift survey is consistent with an obscured AGN population that peaks at z~0.7, although further work is necessary to disentangle the selection effects.

We present a total of 4784 unique broad absorption line quasars from the Sloan Digital Sky Survey Third Data Release. An automated algorithm was used to match a continuum to each quasar and to identify regions of flux at least 10% below the continuum over a velocity range of at least 1000 km/s in the CIV and MgII absorption regions. The model continuum was selected as the best-fit match from a set of template quasar spectra binned in luminosity, emission line width, and redshift, with the power-law spectral index and amount of dust reddening as additional free parameters. We characterize our sample through the traditional ``balnicity'' index and a revised absorption index, as well as through parameters such as the width, outflow velocity, fractional depth and number of troughs. From a sample of 16883 quasars at 1.7<z<4.38, we identify 4386 (26.0%) quasars with broad CIV absorption, of which 1756 (10.4%) satisfy traditional selection criteria. From a sample of 34973 quasars at 0.5<z<2.15, we identify 457 (1.31%) quasars with broad MgII absorption, 191 (0.55%) of which satisfy traditional selection criteria. We also provide a supplementary list of 39 visually identified z>4.38 quasars with broad CIV absorption. We find that broad absorption line quasars may have broader emission lines on average than other quasars.

We study 16,707 quasar spectra from the SDSS using the Karhunen-Lo\`eve (KL) transform (or Principal Component Analysis, PCA). The quasar eigenspectra of the full catalog reveal the following: 1st order - the mean spectrum; 2nd order - a host-galaxy component; 3rd order - the UV-optical continuum slope; 4th order - the correlations of Balmer emission lines. We find that the spectral classification of quasars is redshift and luminosity dependent, as such there does not exist a compact set (< 10 modes) of eigenspectra (covering 900 to 8000 angstrom) which can describe most variations (> 95%) of the entire catalog. We therefore construct several sets of eigenspectra in different redshift and luminosity bins. From these eigenspectra we find that quasar spectra can be classified (by the first two eigenspectra) into a sequence that is defined by a simple progression in the steepness of the slope of the continuum. We also find a dependence on redshift and luminosity in the eigencoefficients. The dominant redshift effect is a result of the evolution of the blended FeII emission (optical) and the Balmer continuum (the ``small bump''). A luminosity dependence is also present in the eigencoefficients and is related to the Baldwin effect. The eigenspectra constructed from the wavelength-selected SDSS spectra are found to agree with the principal components by Francis et al. (1992) and the well-known ``Eigenvector-1'' (Boroson and Green 1992). (Abridged)

We present the second edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog. The catalog consists of 16713 objects in the SDSS First Data Release (DR1) that have luminosities larger than M_i=-22 (H_0=70 km/s, Omega_M=0.3, Omega_Lambda=0.7), have at least one emission line with FWHM larger than 1000 km/s, and have highly reliable redshifts. The area covered is ~1360 deg^2 and the redshifts range from 0.08 to 5.41, with a median value of 1.43. Each object has positions accurate to better than 0.2" rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method, in addition to radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large area surveys. Calibrated digital spectra covering 3800-9200 Angstroms at a spectral resolution 1800-2100, are available. This publication supersedes the first SDSS Quasar Catalog, which was based on material from the SDSS Early Data Release. A summary of corrections to current quasar databases is provided. The majority of the objects were found in SDSS commissioning data using a multicolor selection technique. Since the quasar selection algorithm was undergoing testing during the entire DR1 observational period, care must be taken when assembling samples for use in statistical studies. A total of 15786 objects (94%) in the catalog were discovered by the SDSS; 12,173 of which are reported here for the first time, including five quasars brighter than i=16.0 and 17 quasars with redshifts larger than 4.5.