Quasars are the most luminous of active galactic nuclei (AGN), and are perhaps responsible for quenching star formation in their hosts. The Stripe 82X catalog covers 31.3 deg$^2$ of the Stripe 82 field, of which the 15.6 deg$^2$ covered with XMM-Newton is also covered by Herschel/SPIRE. We have 2500 X-ray detected sources with multi-wavelength counterparts, and 30% of these are unobscured quasars, with $L_X > 10^{44}\,$erg/s and $M_B < -23$. We define a new population of quasars which are unobscured, have X-ray luminosities in excess of $10^{44}\,$erg/s, have broad emission lines, and yet are also bright in the far-infrared, with a 250$\mu$m flux density of $S_{\rm 250}>30$mJy. We refer to these Herschel-detected, unobscured quasars as "Cold Quasars". A mere 4% (21) of the X-ray- and optically-selected unobscured quasars in Stripe 82X are detected at 250$\mu$m. These Cold Quasars lie at $z\sim1-3$, have $L_{\rm IR}>10^{12}\,L_\odot$, and have star formation rates of $\sim200-1400\,M_\odot$/yr. Cold Quasars are bluer in the mid-IR than the full quasar population, and 72% of our Cold Quasars have WISE W3 $<$ 11.5 [Vega], while only 19% of the full quasar sample meets this criteria. Crucially, Cold Quasars have on average $\sim9\times$ as much star formation as the main sequence of star forming galaxies at similar redshifts. Although dust-rich, unobscured quasars have occasionally been noted in the literature before, we argue that they should be considered as a separate class of quasars due to their high star formation rates. This phase is likely short-lived, as the central engine and immense star formation consume the gas reservoir. Cold Quasars are type-1 blue quasars that reside in starburst galaxies.

The clustering of active galactic nuclei (AGN) sheds light on their typical large (Mpc-scale) environments, which can constrain the growth and evolution of supermassive black holes. Here we measure the clustering of luminous X-ray-selected AGN in the Stripe 82X and XMM-XXL-North surveys around the peak epoch of black hole growth, in order to investigate the dependence of luminosity on large-scale AGN environment. We compute the auto-correlation function of AGN in two luminosity bins, $10^{43}\leq L_X<10^{44.5}$ erg s$^{-1}$ at $z\sim 0.8$ and $L_X\geq 10^{44.5}$ erg s$^{-1}$ at $z\sim 1.8$, and calculate the AGN bias taking into account the redshift distribution of the sources using three different methods. Our results show that while the less luminous sample has an inferred typical halo mass that is smaller than for the more luminous AGN, the host halo mass may be less dependent on luminosity than suggested in previous work. Focusing on the luminous sample, we calculate a typical host halo mass of $\sim 10^{13}$ M$_{\odot}~h^{-1}$, which is similar to previous measurements of moderate-luminosity X-ray AGN and significantly larger than the values found for optical quasars of similar luminosities and redshifts. We suggest that the clustering differences between different AGN selection techniques are dominated by selection biases, and not due to a dependence on AGN luminosity. We discuss the limitations of inferring AGN triggering mechanisms from halo masses derived by large-scale bias.

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized.

We constrain X-ray spectral shapes for the ensemble of AGN based on the shape of the Cosmic X-ray Background (CXB). Specifically, we rule out regions of X-ray spectral parameter space that do not reproduce the CXB in the energy range 1-100 keV. The key X-ray spectral parameters are the photon index, {\Gamma}; the cutoff energy, Ecutoff; and the reflection scaling factor, R. Assuming each parameter follows a Gaussian distribution, we first explore the parameter space using a Bayesian approach and a fixed X-ray luminosity function (XLF). For {\sigma}_E = 36 keV and {\sigma}_R = 0.14, fixed at the observed values from the Swift-BAT 70-month sample, we allow <R>, <Ecutoff > and <{\Gamma}> to vary subject to reproducing the CXB. We report results for {\sigma}_{\Gamma} = 0.1-0.5. In an alternative approach, we define the parameter distributions, then forward model to fit the CXB by perturbing the XLF using a neural network. This approach allows us to rule out parameter combinations that cannot reproduce the CXB for any XLF. The marginalized conditional probabilities for the four free parameters are: <R> = 0.99^{+0.11}_{-0.26}, <Ecutoff> = 118^{+24}_{-23}, {\sigma}_{\Gamma} = 0.101^{+0.097}_{-0.001} and <{\Gamma}> = 1.9^{+0.08}_{-0.09}. We provide an interactive online tool for users to explore any combination of <Ecutoff>, {\sigma}_E, <{\Gamma}>, {\sigma}_{\Gamma}, <R> and {\sigma}_R including different distributions for each absorption bin, subject to the integral CXB constraint. The distributions observed in many AGN samples can be ruled out by our analysis, meaning these samples can not be representative of the full AGN population. The few samples that fall within the acceptable parameter space are hard X-ray-selected, commensurate with their having fewer selection biases.

With the launch of eROSITA (extended Roentgen Survey with an Imaging Telescope Array), successfully occurred on 2019 July 13, we are facing the challenge of computing reliable photometric redshifts for 3 million of active galactic nuclei (AGNs) over the entire sky, having available only patchy and inhomogeneous ancillary data. While we have a good understanding of the photo-z quality obtainable for AGN using spectral energy distribution (SED)-fitting technique, we tested the capability of machine learning (ML), usually reliable in computing photo-z for QSO in wide and shallow areas with rich spectroscopic samples. Using MLPQNA as example of ML, we computed photo-z for the X-ray-selected sources in Stripe 82X, using the publicly available photometric and spectroscopic catalogues. Stripe 82X is at least as deep as eROSITA will be and wide enough to include also rare and bright AGNs. In addition, the availability of ancillary data mimics what can be available in the whole sky. We found that when optical, and near- and mid-infrared data are available, ML and SED fitting perform comparably well in terms of overall accuracy, realistic redshift probability density functions, and fraction of outliers, although they are not the same for the two methods. The results could further improve if the photometry available is accurate and including morphological information. Assuming that we can gather sufficient spectroscopy to build a representative training sample, with the current photometry coverage we can obtain reliable photo-z for a large fraction of sources in the Southern hemisphere well before the spectroscopic follow-up, thus timely enabling the eROSITA science return. The photo-z catalogue is released here.

We present spectral energy distributions (SEDs) of 41 active galactic nuclei, derived from multiwavelength photometry and archival spectroscopy. All of the SEDs span at least 0.09 to 30 micron, but in some instances wavelength coverage extends into the X-ray, far-infrared and radio. For some AGNs we have fitted the measured far-infrared photometry with greybody models, while radio flux density measurements have been approximated by power-laws or polynomials. We have been able to fill some of the gaps in the spectral coverage using interpolation or extrapolation of simple models. In addition to the 41 individual AGN SEDs, we have produced 72 Seyfert SEDs by mixing SEDs of the central regions of Seyferts with galaxy SEDs. Relative to the literature, our templates have broader wavelength coverage and/or higher spectral resolution. We have tested the utility of our SEDs by using them to generate photometric redshifts for 0 < z < 6.12 AGNs in the Bootes field (selected with X-ray, IR and optical criteria) and, relative to SEDs from the literature, they produce comparable or better photometric redshifts with reduced flux density residuals.

The cosmic X-ray background (CXB), which peaks at an energy of ~30 keV, is produced primarily by emission from accreting supermassive black holes (SMBHs). The CXB therefore serves as a constraint on the integrated SMBH growth in the Universe and the accretion physics and obscuration in active galactic nuclei (AGNs). This paper gives an overview of recent progress in understanding the high-energy (>~10 keV) X-ray emission from AGNs and the synthesis of the CXB, with an emphasis on results from NASA's NuSTAR hard X-ray mission. We then discuss remaining challenges and open questions regarding the nature of AGN obscuration and AGN physics. Finally, we highlight the exciting opportunities for a next-generation, high-resolution hard X-ray mission to achieve the long-standing goal of resolving and characterizing the vast majority of the accreting SMBHs that produce the CXB.

We report the results of a Sloan Digital Sky Survey-IV eBOSS program to target X-ray sources and mid-infrared-selected WISE AGN candidates in a 36.8 deg$^2$ region of Stripe 82. About half this survey (15.6 deg$^2$) covers the largest contiguous portion of the Stripe 82 X-ray survey. This program represents the largest spectroscopic survey of AGN candidates selected solely by their WISE colors. We combine this sample with X-ray and WISE AGN in the field identified via other sources of spectroscopy, producing a catalog of 4847 sources that is 82% complete to $r\sim22$. Based on X-ray luminosities or WISE colors, 4730 of these sources are AGN, with a median sample redshift of $z\sim1$. About 30% of the AGN are optically obscured (i.e., lack broad lines in their optical spectra). BPT analysis, however, indicates that 50% of the WISE AGN at $z<0.5$ have emission line ratios consistent with star-forming galaxies, so whether they are buried AGN or star-forming galaxy contaminants is currently unclear. We find that 61% of X-ray AGN are not selected as MIR AGN, with 22% of X-ray AGN undetected by WISE. Most of these latter AGN have high X-ray luminosities ($L_{\rm x} > 10^{44}$ erg s$^{-1}$), indicating that MIR selection misses a sizable fraction of the highest luminosity AGN, as well as lower luminosity sources where AGN heated dust is not dominating the MIR emission. Conversely, $\sim$58% of WISE AGN are undetected by X-rays, though we do not find that they are preferentially redder than the X-ray detected WISE AGN.

As matter accretes onto the central supermassive black holes in active galactic nuclei (AGN), X-rays are emitted. We present a population synthesis model that accounts for the summed X-ray emission from growing black holes; modulo the efficiency of converting mass to X-rays, this is effectively a record of the accreted mass. We need this population synthesis model to reproduce observed constraints from X-ray surveys: the X-ray number counts, the observed fraction of Compton-thick AGN [$\log$ (N$_{\rm H}$/cm$^{\rm -2}$) $>$ 24] and the spectrum of the cosmic X-ray background (CXB), after accounting for selection biases. Over the past decade, X-ray surveys by {\it XMM-Newton}, {\it Chandra}, \textit{NuSTAR} and \textit{Swift}-BAT have provided greatly improved observational constraints. We find that no existing X-ray luminosity function (XLF) consistently reproduces all these observations. We take the uncertainty in AGN spectra into account, and use a neural network to compute an XLF that fits all observed constraints, including observed Compton-thick number counts and fractions. This new population synthesis model suggests that, intrinsically, 50$\pm$9\% (56$\pm$7\%) of all AGN within z $\simeq$ 0.1 (1.0) are Compton-thick.

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

Multi-wavelength surveys covering large sky volumes are necessary to obtain an accurate census of rare objects such as high luminosity and/or high redshift active galactic nuclei (AGN). Stripe 82X is a 31.3 deg$^2$ X-ray survey with $Chandra$ and $XMM$-Newton observations overlapping the legacy Sloan Digital Sky Survey (SDSS) Stripe 82 field, which has a rich investment of multi-wavelength coverage from the ultraviolet to the radio. The wide-area nature of this survey presents new challenges for photometric redshifts for AGN compared to previous work on narrow-deep fields because it probes different populations of objects that need to be identified and represented in the library of templates. Here we present an updated X-ray plus multi-wavelength matched catalog, including $Spitzer$ counterparts, and estimated photometric redshifts for 5961 (96% of a total of 6181) X-ray sources, which have a normalized median absolute deviation, $\sigma_{\rm nmad}$ = 0.06 and an outlier fraction, $\eta$ = 13.7%. The populations found in this survey, and the template libraries used for photometric redshifts, provide important guiding principles for upcoming large-area surveys such as $eROSITA$ and 3$XMM$ (in X-ray) and the Large Synoptic Survey Telescope (LSST; optical).

We present results of a ground-based near-infrared campaign with Palomar TripleSpec, Keck NIRSPEC, and Gemini GNIRS to target two samples of reddened active galactic nucleus (AGN) candidates from the 31 deg$^2$ Stripe 82 X-ray survey. One sample, which is $\sim$89\% complete to $K<16$ (Vega), consists of eight confirmed AGNs, four of which were identified with our follow-up program, and is selected to have red $R-K$ colors ($>4$, Vega). The fainter sample ($K>17$, Vega) represents a pilot program to follow-up four sources from a parent sample of 34 that are not detected in the single-epoch SDSS catalog and have {\it WISE} quasar colors. All twelve sources are broad-line AGNs (at least one permitted emission line has a FWHM exceeding 1300 km s$^{-1}$) and span a redshift range $0.59 < z < 2.5$. Half the ($R-K$)-selected AGNs have features in their spectra suggestive of outflows. When comparing these sources to a matched sample of blue Type 1 AGNs, we find the reddened AGNs are more distant ($z > 0.5$) and a greater percentage have high X-ray luminosities ($L_{\rm X,full} > 10^{44}$ erg s$^{-1}$). Such outflows and high luminosities may be consistent with the paradigm that reddened broad-line AGNs represent a transitory phase in AGN evolution as described by the major merger model for black hole growth. Results from our pilot program demonstrate proof-of-concept that our selection technique is successful in discovering reddened quasars at $z > 1$ missed by optical surveys.

The wide-area XMM-XXL X-ray survey is used to explore the fraction of obscured AGN at high accretion luminosities, $L_X (\rm 2-10 \, keV) > 10^{44} \, erg \,s ^{-1}$, and out to redshift $z\approx1.5$. The sample covers an area of about $\rm14\,deg^2$ and provides constraints on the space density of powerful AGN over a wide range of neutral hydrogen column densities extending beyond the Compton-thick limit, $\rm N_H\approx10^{24}\,cm^{-2}$. The fraction of obscured Compton-thin ($\rm N_H=10^{22}-10^{24}\,cm^{-2}$) AGN is estimated to be $\approx0.35$ for luminosities $L_X(\rm 2-10\,keV)>10^{44}\,erg\,s^{-1}$ independent of redshift. For less luminous sources the fraction of obscured Compton-thin AGN increases from $0.45\pm0.10$ at $z=0.25$ to $0.75\pm0.05$ at $z=1.25$. Studies that select AGN in the infrared via template fits to the observed Spectral Energy Distribution of extragalactic sources estimate space densities at high accretion luminosities consistent with the XMM-XXL constraints. There is no evidence for a large population of AGN (e.g. heavily obscured) identified in the infrared and missed at X-ray wavelengths. We further explore the mid-infrared colours of XMM-XXL AGN as a function of accretion luminosity, column density and redshift. The fraction of XMM-XXL sources that lie within the mid-infrared colour wedges defined in the literature to select AGN is primarily a function of redshift. This fraction increases from about 20-30% at z=0.25 to about 50-70% at $z=1.5$.

Using {\em Chandra} observations in the 2.15 deg$^{2}$ COSMOS legacy field, we present one of the most accurate measurements of the Cosmic X-ray Background (CXB) spectrum to date in the [0.3-7] keV energy band. The CXB has three distinct components: contributions from two Galactic collisional thermal plasmas at kT$\sim$0.27 and 0.07 keV and an extragalactic power-law with photon spectral index $\Gamma$=1.45$\pm{0.02}$. The 1 keV normalization of the extragalactic component is 10.91$\pm{0.16}$ keV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ keV$^{-1}$. Removing all X-ray detected sources, the remaining unresolved CXB is best-fit by a power-law with normalization 4.18$\pm{0.26}$ keV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ keV$^{-1}$ and photon spectral index $\Gamma$=1.57$\pm{0.10}$. Removing faint galaxies down to i$_{AB}\sim$27-28 leaves a hard spectrum with $\Gamma\sim$1.25 and a 1 keV normalization of $\sim$1.37 keV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ keV$^{-1}$. This means that $\sim$91\% of the observed CXB is resolved into detected X-ray sources and undetected galaxies. Unresolved sources that contribute $\sim 8-9\%$ of the total CXB show a marginal evidence of being harder and possibly more obscured than resolved sources. Another $\sim$1\% of the CXB can be attributed to still undetected star forming galaxies and absorbed AGN. According to these limits, we investigate a scenario where early black holes totally account for non source CXB fraction and constrain some of their properties. In order to not exceed the remaining CXB and the $z\sim$6 accreted mass density, such a population of black holes must grow in Compton-thick envelopes with N$_{H}>$1.6$\times$10$^{25}$ cm$^{-2}$ and form in extremely low metallicity environments $(Z_\odot)\sim10^{-3}$.

We release the next installment of the Stripe 82 X-ray survey point-source catalog, which currently covers 31.3 deg$^2$ of the Sloan Digital Sky Survey (SDSS) Stripe 82 Legacy field. In total, 6181 unique X-ray sources are significantly detected with {\it XMM-Newton} ($>5\sigma$) and {\it Chandra} ($>4.5\sigma$). This catalog release includes data from {\it XMM-Newton} cycle AO 13, which approximately doubled the Stripe 82X survey area. The flux limits of the Stripe 82X survey are $8.7\times10^{-16}$ erg s$^{-1}$ cm$^{-2}$, $4.7\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$, and $2.1\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$ in the soft (0.5-2 keV), hard (2-10 keV), and full bands (0.5-10 keV), respectively, with approximate half-area survey flux limits of $5.4\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$, $2.9\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$, and $1.7\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$. We matched the X-ray source lists to available multi-wavelength catalogs, including updated matches to the previous release of the Stripe 82X survey; 88\% of the sample is matched to a multi-wavelength counterpart. Due to the wide area of Stripe 82X and rich ancillary multi-wavelength data, including coadded SDSS photometry, mid-infrared {\it WISE} coverage, near-infrared coverage from UKIDSS and VHS, ultraviolet coverage from {\it GALEX}, radio coverage from FIRST, and far-infrared coverage from {\it Herschel}, as well as existing $\sim$30\% optical spectroscopic completeness, we are beginning to uncover rare objects, such as obscured high-luminosity AGN at high-redshift. The Stripe 82X point source catalog is a valuable dataset for constraining how this population grows and evolves, as well as for studying how they interact with the galaxies in which they live.

We present an extensive search for RR Lyrae stars in and around the ultra-faint Milky Way companions Segue 2 and Segue 3. The former (M_V = -2.5, Belokurov et al. 2009) appears to be an extremely faint dwarf galaxy companion of the Milky Way. The latter (M_V = 0.0, Fadely et al. 2011) is among the faintest star clusters known. We use B and V band time-series imaging obtained at the WIYN 0.9 meter telescope at Kitt Peak National Observatory to search for RR Lyrae in these objects. In our Segue 2 observations, we present a previously unknown fundamental mode (RRab) RR Lyrae star with a period of P_ab = 0.748 days. With this measurement, we revisit the inverse correlation between <P_ab> and <[Fe/H]> established in the literature for Milky Way dwarf galaxies and their RR Lyrae. In this context, the long period of Segue 2's RRab star as well as the known significant spread in metallicity in this dwarf galaxy are consistent with the observed trend in <P_ab> and <[Fe/H]>. We derive the first robust distance to Segue 2, using both its RRab star and spectroscopically confirmed blue horizontal branch stars. Using [Fe/H] = -2.16 and -2.44 dex, we find d_RRL = 36.6 +2.5/-2.4 and 37.7 +/- 2.7 kpc; assuming [Fe/H] = -2.257 dex, we find d_BHB = 34.4 +/- 2.6 kpc. Although no RR Lyrae were present in the Segue 3 field, we found a candidate eclipsing binary star system.

The Hubble Space Telescope (HST) Treasury Program on the Orion Nebula Cluster has used 104 orbits of HST time to image the Great Orion Nebula region with the Advanced Camera for Surveys (ACS), the Wide-Field/Planetary Camera 2 (WFPC2) and the Near Infrared Camera and Multi Object Spectrograph (NICMOS) instruments in 11 filters ranging from the U-band to the H-band equivalent of HST. The program has been intended to perform the definitive study of the stellar component of the ONC at visible wavelengths, addressing key questions like the cluster IMF, age spread, mass accretion, binarity and cirumstellar disk evolution. The scanning pattern allowed to cover a contiguous field of approximately 600 square arcminutes with both ACS and WFPC2, with a typical exposure time of approximately 11 minutes per ACS filter, corresponding to a point source depth AB(F435W) = 25.8 and AB(F775W)=25.2 with 0.2 magnitudes of photometric error. We describe the observations, data reduction and data products, including images, source catalogs and tools for quick look preview. In particular, we provide ACS photometry for 3399 stars, most of them detected at multiple epochs, WFPC2 photometry for 1643 stars, 1021 of them detected in the U-band, and NICMOS JH photometry for 2116 stars. We summarize the early science results that have been presented in a number of papers. The final set of images and the photometric catalogs are publicly available through the archive as High Level Science Products at the STScI Multimission Archive hosted by the Space Telescope Science Institute.