We report on observations of redshifted CO(1-0) line emission and observed-frame $\rm\sim$ 30GHz radio continuum emission from five ultra-luminous, mid-IR selected hot, Dust-Obscured Galaxies (Hot DOGs) at $z\rm\gtrsim$ 3 using the Karl G. Jansky Very Large Array. We detect CO(1-0) line emission in all five Hot DOGs, with one of them at high signal to noise. We analyse FIR-radio spectral energy distributions, including dust, free-free and synchrotron emission for the galaxies. We find that most of the 115 GHz rest-frame continuum is mostly due to synchrotron or free-free emission, with only a potentially small contribution from thermal emission. We see a deficit in the rest-frame 115 GHz continuum emission compared to dusty star-forming galaxies (DSFGs) and sub-millimetre galaxies (SMGs) at high redshift, suggesting that Hot DOGs do not have similar cold gas reserves compared with star-forming galaxies. One target, W2305-0039, is detected in the FIRST 1.4 GHz survey, and is likely to possess compact radio jets. We compare to the FIR-radio correlation, and find that at least half of the Hot DOGs in our sample are radio-quiet with respect to normal galaxies. These findings suggest that Hot DOGs have comparably less cold molecular gas than star-forming galaxies at lower, $z\rm\sim$ 2 redshifts, and are dominated by powerful, yet radio-quiet AGN.

Hot Dust-Obscured Galaxies (Hot DOGs) are among the most luminous galaxies in the Universe. Powered by highly obscured, possibly Compton-thick, active galactic nuclei (AGNs), Hot DOGs are characterized by SEDs that are very red in the mid-IR yet dominated by the host galaxy stellar emission in the UV and optical. An earlier study identified a sub-sample of Hot DOGs with significantly enhanced UV emission. One target, W0204-0506, was studied in detail and, based on Chandra observations, it was concluded that the enhanced emission was most likely due to either extreme unobscured star-formation (${\rm SFR}>1000~M_{\odot}~\rm yr^{-1}$) or to light from the highly obscured AGN scattered by gas or dust into our line of sight. Here, we present a follow-up study of W0204-0506 as well as two more Hot DOGs with excess UV emission. For the two new objects we obtained Chandra/ACIS-S observations, and for all three targets we obtained HST/WFC3 F555W and F160W imaging. We conclude that the excess UV emission is primarily dominated by light from the central highly obscured, hyper-luminous AGN that has been scattered into our line of sight. We cannot rule out, however, that star-formation may significantly contribute to the UV excess of W0204-0506.

The intergalactic medium (IGM) at $z\sim$ 5 to 6 is largely ionized, and yet the main source for the IGM ionization in the early universe is uncertain. Of the possible contributors are faint quasars with $-26 \lesssim M_{\rm 1450} \lesssim -23$, but their number density is poorly constrained at $z\sim5$. In this paper, we present our survey of faint quasars at $z\sim5$ in the European Large-Area {\it ISO} Survey-North 1 (ELAIS-N1) field over a survey area of 6.51 deg$^2$ and examine if such quasars can be the dominant source of the IGM ionization. We use the deep optical/near-infrared data of the ELAIS-N1 field as well as the additional medium-band observations to find $z \sim 5$ quasars through a two-step approach using the broadband color selection, and SED fitting with the medium-band information included. Adopting Bayesian information criterion, we identify ten promising quasar candidates. Spectra of three of the candidates are obtained, confirming all of them to be quasars at $z\sim5$ and supporting the reliability of the quasar selection. Using the promising candidates, we derive the $z\sim5$ quasar luminosity function at $-26 \lesssim M_{\rm 1450} \lesssim -23$. The number density of faint $z\sim5$ quasars in the ELAIS-N1 field is consistent with several previous results that quasars are not the main contributors to the IGM-ionizing photons at $z\sim5$

We present analysis of Chandra X-ray observations of seven quasars that were identified as candidate sub-parsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey (CRTS) based on apparent periodicity in their optical light curves. Simulations predict close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs, including harder or softer X-ray spectra, ripple-like profiles in the Fe K-$\alpha$ line, and distinct peaks in the spectrum due to the separation of the accretion disk into a circumbinary disk and mini-disks around each SMBH. We obtained Chandra observations to test these models and assess whether these quasars could contain binary SMBHs. We instead find that the quasar spectra are all well fit by simple absorbed power law models, with the rest frame 2-10 keV photon indices, $\Gamma$, and the X-ray-to-optical power slopes, $\alpha_{\rm OX}$, indistinguishable from the larger quasar population. This may indicate that these seven quasars are not truly sub-parsec binary SMBH systems, or it may simply reflect that our sample size was too small to robustly detect any differences. Alternatively, the X-ray spectral changes might only be evident at higher energies than probed by Chandra. Given the available models and current data, no firm conclusions are drawn. These observations will help motivate and direct further work on theoretical models of binary SMBH systems, such as modeling systems with thinner accretion disks and larger binary separations.

We present VLT/XSHOOTER rest-frame UV-optical spectra of 10 Hot Dust-Obscured Galaxies (Hot DOGs) at $z\sim2$ to investigate AGN diagnostics and to assess the presence and effect of ionized gas outflows. Most Hot DOGs in this sample are narrow-line dominated AGN (type 1.8 or higher), and have higher Balmer decrements than typical type 2 quasars. Almost all (8/9) sources show evidence for ionized gas outflows in the form of broad and blueshifted [O III] profiles, and some sources have such profiles in H$\alpha$ (5/7) or [O II] (3/6). Combined with the literature, these results support additional sources of obscuration beyond the simple torus invoked by AGN unification models. Outflow rates derived from the broad [O III] line ($\rm \gtrsim10^{3}\,M_{\odot}\,yr^{-1}$) are greater than the black hole accretion and star formation rates, with feedback efficiencies ($\sim0.1-1\%$) consistent with negative feedback to the host galaxy's star formation in merger-driven quasar activity scenarios. We find the broad emission lines in luminous, obscured quasars are often better explained by outflows within the narrow line region, and caution that black hole mass estimates for such sources in the literature may have substantial uncertainty. Regardless, we find lower bounds on the Eddington ratio for Hot DOGs near unity.

The discovery of luminous quasars at redshifts up to 7.5 demonstrates the existence of several billion M_sun supermassive black holes (SMBHs) less than a billion years after the Big Bang. They are accompanied by intense star formation in their host galaxies, pinpointing sites of massive galaxy assembly in the early universe, while their absorption spectra reveal an increasing neutral intergalactic medium (IGM) at the epoch of reionization. Extrapolating from the rapid evolution of the quasar density at z=5-7, we expect that there is only one luminous quasar powered by a billion M_sun SMBH in the entire observable universe at z~9. In the next decade, new wide-field, deep near-infrared (NIR) sky surveys will push the redshift frontier to the first luminous quasars at z~9-10; the combination with new deep X-ray surveys will probe fainter quasar populations that trace earlier phases of SMBH growth. The identification of these record-breaking quasars, and the measurements of their BH masses and accretion properties require sensitive spectroscopic observations with next generation of ground-based and space telescopes at NIR wavelengths. High-resolution integral-field spectroscopy at NIR, and observations at millimeter and radio wavelengths, will together provide a panchromatic view of the quasar host galaxies and their galactic environment at cosmic dawn, connecting SMBH growth with the rise of the earliest massive galaxies. Systematic surveys and multiwavelength follow-up observations of the earliest luminous quasars will strongly constrain the seeding and growth of the first SMBHs in the universe, and provide the best lines of sight to study the history of reionization.

Black hole mass scaling relations suggest that extremely massive black holes (EMBHs) with $M_\mathrm{BH}\ge10^{9.4}\,M_{\odot}$ are found in the most massive galaxies with $M_\mathrm{star}\ge10^{11.6}\,M_{\odot}$, which are commonly found in dense environments, like galaxy clusters. Therefore, one can expect that there is a close connection between active EMBHs and dense environments. Here, we study the environments of 9461 galaxies and 2943 quasars at $0.24 \le z \le 0.40$, among which 52 are extremely massive quasars with $\log(M_\mathrm{BH}/M_{\odot}) \ge 9.4$, using Sloan Digital Sky Survey and MMT Hectospec data. We find that, on average, both massive quasars and massive galaxies reside in environments more than $\sim2$ times as dense as those of their less massive counterparts with $\log(M_\mathrm{BH}/M_{\odot}) \le 9.0$. However, massive quasars reside in environments about half as dense as inactive galaxies with $\log(M_\mathrm{BH}/M_{\odot}) \ge 9.4$, and only about one third of massive quasars are found in galaxy clusters, while about two thirds of massive galaxies reside in such clusters. This indicates that massive galaxies are a much better signpost for galaxy clusters than massive quasars. The prevalence of massive quasars in moderate to low density environments is puzzling, considering that several simulation results show that these quasars appear to prefer dense environments. Several possible reasons for this discrepancy are discussed, although further investigation is needed to obtain a definite explanation.

Recently Assef et al.(2018) presented two catalogs of AGN candidates over 30,093 deg^2 selected from the Wide-field Infrared Survey Explorer (WISE) observations. From their most reliable sample, Assef et al. (2018) identified 45 AGN candidates with the highest variability levels in the AllWISE catalog, but that are not blazars. Here we present new spectroscopic observations of some of these targets to further constrain their nature. We also study their optical lightcurves using observations from CRTS, and find that only seven show significant optical variability, and that five of those seven are spectroscopically classified as AGN. In one of them, WISEA J094806.56+031801.7 (W0948+0318), we identify a transient event in the CRTS lightcurve. We present a detailed analysis of this transient, and characterize it through its CRTS lightcurve and its multi-wavelength spectral energy distribution obtained from GALEX, Pan-STARRS and WISE observations. We find that the most likely source of the transient is a super-luminous supernova (SLSN) in W0948+0318. We estimate the total radiated energy to be E=1.6\pm 0.3 x 10^52 erg, making it one of the most energetic SLSN observed. Based on the lack of change in mid-IR color throughout and after the transient event, we speculate that the location of the SLSN is within the torus of the AGN. We identify 9 possible analogs to W0948+0318 based on their WISE lightcurves. None show optically detected transients and hence suggest significant dust obscuration. Finally, we estimate a rate of >2x10^-7 yr^-1 per AGN for these transients under the conservative assumption that none of the identified analogs have a common origin with the transient in W0948+0318.

The faint quasars with $M_{1450}>-24$ mag are known to hold the key to the determination of the ultraviolet emissivity for the cosmic re-ionization. But only a few have been identified so far because of the limitations on the survey data. Here, we present the first results of the $z\sim5$ faint quasar survey with the Infrared Medium-deep Survey (IMS), which covers $\sim100$ deg$^{2}$ areas in $J$-band to the depths of $J_{\rm AB}\sim23$ mag. To improve selection methods, the medium-band follow-up imaging has been carried out using the SED camera for QUasars in Early uNiverse (SQUEAN) on the Otto Struve 2.1 m Telescope. The optical spectra of the candidates were obtained with 8-m class telescopes. We newly discovered 10 quasars with $-25<M_{1450}<-23$ at $z\sim5$, among which three have been missed in a previous survey using the same optical data over the same area, implying the necessity for improvements in high redshift faint quasars selection. We derived photometric redshifts from the medium-band data, and find that they have high accuracies of $\langle|\Delta z|/(1+z)\rangle=0.016$. The medium-band-based approach allows us to rule out many of the interlopers that contaminate $\gtrsim20~\%$ of the broad-band-selected quasar candidates. These results suggest that the medium-band-based approach is a powerful way to identify $z\sim5$ quasars and measure their redshifts at high accuracy (1-2 %). It is also a cost-effective way to understand the contribution of quasars to the cosmic re-ionization history.

Galaxy mergers and gas accretion from the cosmic web drove the growth of galaxies and their central black holes at early epochs. We report spectroscopic imaging of a multiple merger event in the most luminous known galaxy, WISE J224607.56-052634.9 (W2246-0526), a dust-obscured quasar at redshift 4.6, 1.3 Gyr after the Big Bang. Far-infrared dust continuum observations show three galaxy companions around W2246-0526 with disturbed morphologies, connected by streams of dust likely produced by the dynamical interaction. The detection of tidal dusty bridges shows that W2246-0526 is accreting its neighbors, suggesting merger activity may be a dominant mechanism through which the most luminous galaxies simultaneously obscure and feed their central supermassive black holes.

We use optical and near-infrared spectroscopy to observe rest-UV emission lines and estimate the black hole mass of WISEA J224607.56-052634.9 (W2246-0526) at z = 4.601, the most luminous hot dust-obscured galaxy yet discovered by WISE. From the broad component of the MgII-2799A emission line, we measure a black hole mass of log (M_BH/M_sun) = 9.6 +- 0.4. The broad CIV-1549A line is asymmetric and significantly blue-shifted. The derived M_BH from the blueshift-corrected broad CIV line width agrees with the MgII result. From direct measurement using a well-sampled SED, the bolometric luminosity is 3.6 * 10^14 L_sun. The corresponding Eddington ratio for W2246-0526 is lambda_Edd = L_AGN / L_Edd = 2.8. This high Eddington ratio may reach the level where the luminosity is saturating due to photon trapping in the accretion flow, and be insensitive to the mass accretion rate. In this case, the M_BH growth rate in W2246-0526 would exceed the apparent accretion rate derived from the observed luminosity.

Changing-look quasars are a recently identified class of active galaxies in which the strong UV continuum and/or broad optical hydrogen emission lines associated with unobscured quasars either appear or disappear on timescales of months to years. The physical processes responsible for this behaviour are still debated, but changes in the black hole accretion rate or accretion disk structure appear more likely than changes in obscuration. Here we report on four epochs of spectroscopy of SDSS J110057.70-005304.5, a quasar at a redshift of $z=0.378$ whose UV continuum and broad hydrogen emission lines have faded, and then returned over the past $\approx$20 years. The change in this quasar was initially identified in the infrared, and an archival spectrum from 2010 shows an intermediate phase of the transition during which the flux below rest-frame $\approx$3400\AA\ has decreased by close to an order of magnitude. This combination is unique compared to previously published examples of changing-look quasars, and is best explained by dramatic changes in the innermost regions of the accretion disk. The optical continuum has been rising since mid-2016, leading to a prediction of a rise in hydrogen emission line flux in the next year. Increases in the infrared flux are beginning to follow, delayed by a $\sim$3 year observed timescale. If our model is confirmed, the physics of changing-look quasars are governed by processes at the innermost stable circular orbit (ISCO) around the black hole, and the structure of the innermost disk. The easily identifiable and monitored changing-look quasars would then provide a new probe and laboratory of the nuclear central engine.

Active Galactic Nucleus (AGN) variability can be used to study the physics of the region in the vicinity of the central black hole. In this paper, we investigated intra-night optical variability of AGN in the COSMOS field in order to understand the AGN instability at the smallest scale. Observations were performed using the KMTNet on three separate nights for 2.5-5 hour at a cadence of 20-30 min. We find that the observation enables the detection of the short-term variability as small as $\sim$ 0.02 and 0.1 mag for $R \sim$ 18 and 20 mag sources, respectively. Using four selection methods (X-rays, mid-infrared, radio, and matching with SDSS quasars), 394 AGNs are detected in the 4 deg$^2$ field of view. After differential photometry and $\chi^2-$test, we classify intra-night variable AGNs. But the fraction of variable AGNs (0-8 %) is consistent with a statistical fluctuation from null result. Eight out of 394 AGNs are found to be intra-night variable in two filters or two nights with a variability level of 0.1 mag, suggesting that they are strong candidates for intra-night variable AGNs. Still they represent a small population (2 %). There is no sub-category of AGNs that shows a statistically significant intra-night variability.

SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$\mu$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$\sigma$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop.

We report a new changing-look quasar, WISE~J105203.55+151929.5 at $z=0.303$, found by identifying highly mid-IR variable quasars in the WISE/NEOWISE data stream. Compared to multi-epoch mid-IR photometry of a large sample of SDSS-confirmed quasars, WISE J1052+1519 is an extreme photometric outlier, fading by more than a factor of two at $3.4$ and $4.6 \mu$m since 2009. Swift target-of-opportunity observations in 2017 show even stronger fading in the soft X-rays compared to the ROSAT detection of this source in 1995, with at least a factor of fifteen decrease. We obtained second-epoch spectroscopy with the Palomar telescope in 2017 which, when compared with the 2006 archival SDSS spectrum, reveals that the broad H$\beta$ emission has vanished and that the quasar has become significantly redder. The two most likely interpretations for this dramatic change are source fading or obscuration, where the latter is strongly disfavored by the mid-IR data. We discuss various physical scenarios that could cause such changes in the quasar luminosity over this timescale, and favor changes in the innermost regions of the accretion disk that occur on the thermal and heating/cooling front timescales. We discuss possible physical triggers that could cause these changes, and predict the multiwavelength signatures that could distinguish these physical scenarios.

To date, most of the luminous quasars known at $z\sim6$ have been found to be in maximal accretion with the Eddington ratios, $\lambda_{\rm{Edd}}\sim1$, suggesting enhanced nuclear activities in the early universe. However, this may not be the whole picture of supermassive black hole (SMBH) growth since previous studies have not reached on faint quasars that are more likely to harbor SMBHs with low $\lambda_{\rm{Edd}}$. To gain a better understanding on the accretion activities in quasars in the early universe, we obtained a deep near-infrared (NIR) spectrum of a quasar, IMS J220417.92+011144.8 (hereafter IMS J2204+0112), one of the faintest quasars that have been identified at $z\sim6$. From the redshifted C IV $\lambda 1549$ emission line in the NIR spectrum, we find that IMS J2204+0112 harbors a SMBH with about a billion solar mass and $\lambda_{\rm{Edd}} \sim 0.1$, but with a large uncertainty in both quantities (0.41 dex). IMS J2204+0112 has one of the lowest Eddington ratios among quasars at $z\sim6$, but a common value among quasars at $z\sim2$. Its low $\lambda_{\rm{Edd}}$ can be explained with two scenarios; the SMBH growth from a stellar mass black hole through short-duration super-Eddington accretion events or from a massive black hole seed ($\sim10^{5}\,M_{\odot}$) with Eddington-limited accretion. NIR spectra of more faint quasars are needed to better understand the accretion activities of SMBHs at $z \sim 6$.

Hot, Dust-Obscured Galaxies, or "Hot DOGs", are a rare, dusty, hyperluminous galaxy population discovered by the WISE mission. Predominantly at redshifts 2-3, they include the most luminous known galaxies in the universe. Their high luminosities likely come from accretion onto highly obscured super massive black holes (SMBHs). We have conducted a pilot survey to measure the SMBH masses of five z~2 Hot DOGs via broad H_alpha emission lines, using Keck/MOSFIRE and Gemini/FLAMINGOS-2. We detect broad H_alpha emission in all five Hot DOGs. We find substantial corresponding SMBH masses for these Hot DOGs (~ 10^{9} M_sun), and their derived Eddington ratios are close to unity. These z~2 Hot DOGs are the most luminous AGNs at given BH masses, suggesting they are accreting at the maximum rates for their BHs. A similar property is found for known z~6 quasars. Our results are consistent with scenarios in which Hot DOGs represent a transitional, high-accretion phase between obscured and unobscured quasars. Hot DOGs may mark a special evolutionary stage before the red quasar and optical quasar phases, and they may be present at other cosmic epochs.

We present two large catalogs of AGN candidates identified across ~75% of the sky from the Wide-field Infrared Survey Explorer's AllWISE Data Release. Both catalogs, some of the largest such catalogs published to date, are selected purely on the basis of mid-IR photometry in the WISE W1 and W2 bands. The catalogs are designed to be appropriate for a broad range of scientific investigations, with one catalog emphasizing reliability while the other emphasizes completeness. Specifically, the R90 catalog consists of 4,543,530 AGN candidates with 90% reliability, while the C75 catalog consists of 20,907,127 AGN candidates with 75% completeness. We provide a detailed discussion of potential artifacts, and excise portions of the sky close to the Galactic Center, Galactic Plane, nearby galaxies, and other expected contaminating sources. Our final catalogs cover 30,093 deg^2 of extragalactic sky. These catalogs are expected to enable a broad range of science, and we present a few simple illustrative cases. From the R90 sample we identify 45 highly variable AGN lacking radio counterparts in the FIRST survey, implying they are unlikely to be blazars. One of these sources, WISEA J142846.71+172353.1, is a mid-IR-identified changing-look quasar at z=0.104. We characterize our catalogs by comparing them to large, wide-area AGN catalogs in the literature, specifically UV-to-near-IR quasar selections from SDSS and XDQSOz, mid-IR selection from Secrest et al. (2015) and X-ray selection from ROSAT. From the latter work, we identify four ROSAT X-ray sources that each are matched to three WISE-selected AGN in the R90 sample within 30". Palomar spectroscopy reveals one of these systems, 2RXS J150158.6+691029, to consist of a triplet of quasars at z=1.133 +/- 0.004, suggestive of a rich group or forming galaxy cluster.(Abridged)

We analyze the evolution of the dual AGN in ESO 509-IG066, a galaxy pair located at z=0.034 whose nuclei are separated by 11 kpc. Previous observations with XMM-Newton on this dual AGN found evidence for two moderately obscured ($N_H\sim10^{22}$ cm$^{-2}$) X-ray luminous ($L_X\sim10^{43}$ erg/s) nuclear sources. We present an analysis of subsequent Chandra, NuSTAR and Swift/XRT observations that show one source has dropped in flux by a factor of 10 between 2004 and 2011, which could be explained by either an increase in the absorbing column or an intrinsic fading of the central engine possibly due to a decrease in mass accretion. Both of these scenarios are predicted by galaxy merger simulations. The source which has dropped in flux is not detected by NuSTAR, which argues against absorption, unless it is extreme. However, new Keck/LRIS optical spectroscopy reveals a previously unreported broad H-alpha line which is highly unlikely to be visible under the extreme absorption scenario. We therefore conclude that the black hole in this nucleus has undergone a dramatic drop in accretion rate. From AO-assisted near-infrared integral-field spectroscopy of the other nucleus, we find evidence that the galaxy merger is having a direct effect on the kinematics of the gas close to the nucleus of the galaxy, providing a direct observational link between the galaxy merger and the mass accretion rate on to the black hole.

We present our first results of the survey for high redshift quasars at $5 \lesssim {\rm z} \lesssim 5.7$. The search for quasars in this redshift range has been known to be challenging due to limitations of filter sets used in previous studies. We conducted a quasar survey for two specific redshift ranges, 4.60 $\leq$ z $\leq$ 5.40 and 5.50 $\leq$ z $\leq$ 6.05, using multi-wavelength data that include observations using custom-designed filters, $is$ and $iz$. Using these filters and a new selection technique, we were able to reduce the fraction of interlopers. Through optical spectroscopy, we confirmed seven quasars at 4.7 $\leq$ z $\leq$ 5.4 with $-27.4 < M_{1450} < -26.4$ which were discovered independently by another group recently. We estimated black hole masses and Eddington ratios of four of these quasars from optical and near-infrared spectra, and found that these quasars are undergoing nearly Eddington-limited accretion which is consistent with the rapid growth of supermassive black holes in luminous quasars at z $\sim$ 5.

CRTS J084133.15+200525.8 is an optically bright quasar at z=2.345 that has shown extreme spectral variability over the past decade. Photometrically, the source had a visual magnitude of V~17.3 between 2002 and 2008. Then, over the following five years, the source slowly brightened by approximately one magnitude, to V~16.2. Only ~1 in 10,000 quasars show such extreme variability, as quantified by the extreme parameters derived for this quasar assuming a damped random walk model. A combination of archival and newly acquired spectra reveal the source to be an iron low-ionization broad absorption line (FeLoBAL) quasar with extreme changes in its absorption spectrum. Some absorption features completely disappear over the 9 years of optical spectra, while other features remain essentially unchanged. We report the first definitive redshift for this source, based on the detection of broad H-alpha in a Keck/MOSFIRE spectrum. Absorption systems separated by several 1000 km/s in velocity show coordinated weakening in the depths of their troughs as the continuum flux increases. We interpret the broad absorption line variability to be due to changes in photoionization, rather than due to motion of material along our line of sight. This source highlights one sort of rare transition object that astronomy will now be finding through dedicated time-domain surveys.

We obtained near-infrared spectra of 26 SDSS quasars at $0.7<z<2.5$ with reported rest-frame ultraviolet $M_{\rm BH} \sim 10^{10}M_{\odot}$ to critically examine the systematic effects involved with their mass estimations. We find that AGNs heavier than $10^{10}M_{\odot}$ often display double-peaked H$\alpha$ emission, extremely broad FeII complex emission around MgII, and highly blueshifted and broadened CIV emission. The weight of this evidence, combined with previous studies, cautions against the use of $M_{\rm BH}$ values based on any emission line with a width over 8000 km/s. Also, the $M_{\rm BH}$ estimations are not positively biased along the presence of ionized narrow line outflows, anisotropic radiation, or the use of line FWHM instead of $\sigma$ for our sample, and unbiased with variability, scatter in broad line equivalent width, or obscuration for general type-1 quasars. Removing the systematically uncertain $M_{\rm BH}$ values, $\sim10^{10}M_{\odot}$ BHs in $1\lesssim z \lesssim 2$ AGNs can still be explained by anisotropic motion of the broad line region from $\sim10^{9.5}M_{\odot}$ BHs, although current observations support they are intrinsically most massive, and overmassive to the host's bulge mass.

We present the first full catalog and science results for the NuSTAR serendipitous survey. The catalog incorporates data taken during the first 40 months of NuSTAR operation, which provide ~20Ms of effective exposure time over 331 fields, with an areal coverage of 13 sq deg, and 497 sources detected in total over the 3-24 keV energy range. There are 276 sources with spectroscopic redshifts and classifications, largely resulting from our extensive campaign of ground-based spectroscopic followup. We characterize the overall sample in terms of the X-ray, optical, and infrared source properties. The sample is primarily comprised of active galactic nuclei (AGNs), detected over a large range in redshift from z = 0.002 - 3.4 (median of <z> = 0.56), but also includes 16 spectroscopically confirmed Galactic sources. There is a large range in X-ray flux, from log( f_3-24keV / erg s^-1 cm^-2 ) ~ -14 to -11, and in rest-frame 10-40 keV luminosity, from log( L_10-40keV / erg s^-1 ) ~ 39 to 46, with a median of 44.1. Approximately 79% of the NuSTAR sources have lower energy (<10 keV) X-ray counterparts from XMM-Newton, Chandra, and Swift/XRT. The mid-infrared (MIR) analysis, using WISE all-sky survey data, shows that MIR AGN color selections miss a large fraction of the NuSTAR-selected AGN population, from ~15% at the highest luminosities (Lx > 10^44 erg s^-1) to ~80% at the lowest luminosities (Lx < 10^43 erg s^-1). Our optical spectroscopic analysis finds that the observed fraction of optically obscured AGNs (i.e., the Type 2 fraction) is F_Type2 = 53(+14-15)%, for a well-defined subset of the 8-24 keV selected sample. This is higher, albeit at a low significance level, than the Type 2 fraction measured for redshift- and luminosity-matched AGNs selected by <10 keV X-ray missions.

Hot, Dust-Obscured Galaxies (Hot DOGs), selected from the WISE all sky infrared survey, host some of the most powerful Active Galactic Nuclei (AGN) known, and might represent an important stage in the evolution of galaxies. Most known Hot DOGs are at $z> 1.5$, due in part to a strong bias against identifying them at lower redshift related to the selection criteria. We present a new selection method that identifies 153 Hot DOG candidates at $z\sim 1$, where they are significantly brighter and easier to study. We validate this approach by measuring a redshift $z=1.009$, and an SED similar to higher redshift Hot DOGs for one of these objects, WISE J1036+0449 ($L_{\rm\,Bol}\simeq 8\times 10^{46}\rm\,erg\,s^{-1}$), using data from Keck/LRIS and NIRSPEC, SDSS, and CSO. We find evidence of a broadened component in MgII, which, if due to the gravitational potential of the supermassive black hole, would imply a black hole mass of $M_{\rm\,BH}\simeq 2 \times 10^8 M_{\odot}$, and an Eddington ratio of $\lambda_{\rm\,Edd}\simeq 2.7$. WISE J1036+0449 is the first Hot DOG detected by NuSTAR, and the observations show that the source is heavily obscured, with a column density of $N_{\rm\,H}\simeq(2-15)\times10^{23}\rm\,cm^{-2}$. The source has an intrinsic 2-10 keV luminosity of $\sim 6\times 10^{44}\rm\,erg\,s^{-1}$, a value significantly lower than that expected from the mid-infrared/X-ray correlation. We also find that the other Hot DOGs observed by X-ray facilities show a similar deficiency of X-ray flux. We discuss the origin of the X-ray weakness and the absorption properties of Hot DOGs. Hot DOGs at $z\lesssim1$ could be excellent laboratories to probe the characteristics of the accretion flow and of the X-ray emitting plasma at extreme values of the Eddington ratio.

We present an analysis of the joint $XMM$-$Newton$ and $NuSTAR$ observations of the radio-quiet quasar QSO B2202$-$209. Using an optical observation from the Hale Telescope at the Palomar Observatory, we revise the redshift of the source from the previously reported $z=1.77$ to $z=0.532$, and we estimate the mass of the central black hole, $\log (M_{\rm BH}/M_{\odot}) = 9.08 \pm 0.18$. The X-ray spectrum of this source can be well described by a power-law of photon index $\Gamma = 1.82 \pm 0.05$ with $E_{\rm cut} = 152_{-54}^{+103}\,{\rm keV}$, in the rest frame of the source. Assuming a Comptonisation model, we estimate the coronal temperature to be $kT_{\rm e}=42\pm 3 \,{\rm keV}$ and $kT_{\rm e}= 56 \pm 3\,{\rm keV}$ for a spherical and a slab geometry, respectively. The coronal properties are comparable to the ones derived for local AGN, despite a difference of around one order of magnitude in black hole mass and X-ray luminosity ($L_{2-10} = 1.93\times 10^{45}\,{\rm erg\,s^{-1}}$). The quasar is X-ray loud, with an unusually flat observed optical-to-X-ray spectral slope $\alpha_{\rm OX} = 1.00 \pm 0.02$, and has an exceptionally strong optical [O III] line. Assuming that both the X-ray emission and the [O III] line are isotropic, these two extreme properties can be explained by a nearly edge-on disk, leading to a reduction in the observed UV continuum light.

Luminous quasars at z>5.6 can be studied in detail with the current generation of telescopes and provide us with unique information on the first gigayear of the universe. Thus far these studies have been statistically limited by the number of quasars known at these redshifts. Such quasars are rare and therefore wide-field surveys are required to identify them and multiwavelength data are needed to separate them efficiently from their main contaminants, the far more numerous cool dwarfs. In this paper, we update and extend the selection for z~6 quasars presented in Banados et al. (2014) using the Pan-STARRS1 (PS1) survey. We present the PS1 distant quasar sample, which currently consists of 124 quasars in the redshift range 5.6<z<6.7 that satisfy our selection criteria. Seventy-seven of these quasars have been discovered with PS1, and 63 of them are newly identified in this paper. We present composite spectra of the PS1 distant quasar sample. This sample spans a factor of ~20 in luminosity and shows a variety of emission line properties. The number of quasars at z>5.6 presented in this work almost double the quasars previously known at these redshifts, marking a transition phase from studies of individual sources to statistical studies of the high-redshift quasar population, which was impossible with earlier, smaller samples.

We investigate the process of reionisation in a model in which the dark matter is a warm elementary particle such as a sterile neutrino. We focus on models that are consistent with the dark matter decay interpretation of the recently detected line at 3.5 keV in the X-ray spectra of galaxies and clusters. In warm dark matter models the primordial spectrum of density perturbations has a cut-off on the scale of dwarf galaxies. Structure formation therefore begins later than in the standard cold dark matter (CDM) model and very few objects form below the cut-off mass scale. To calculate the number of ionising photons, we use the Durham semi-analytic model of galaxy formation, GALFORM. We find that even the most extreme 7 keV sterile neutrino we consider is able to reionise the Universe early enough to be compatible with the bounds on the epoch of reionisation from Planck. This, perhaps surprising, result arises from the rapid build-up of high redshift galaxies in the sterile neutrino models which is also reflected in a faster evolution of their far-UV luminosity function between $10>z>7$ than in CDM. The dominant sources of ionising photons are systematically more massive in the sterile neutrino models than in CDM. As a consistency check on the models, we calculate the present-day luminosity function of satellites of Milky Way-like galaxies. When the satellites recently discovered in the DES survey are taken into account, strong constraints are placed on viable sterile neutrino models.

The $\Lambda$CDM cosmological model successfully reproduces many aspects of the galaxy and structure formation of the universe. However, the growth of large-scale structures (LSSs) in the early universe is not well tested yet with observational data. Here, we have utilized wide and deep optical--near-infrared data in order to search for distant galaxy clusters and superclusters ($0.8<z<1.2$). From the spectroscopic observation with the Inamori Magellan Areal Camera and Spectrograph (IMACS) on the Magellan telescope, three massive clusters at $z\sim$0.91 are confirmed in the SSA22 field. Interestingly, all of them have similar redshifts within $\Delta z\sim$0.01 with velocity dispersions ranging from 470 to 1300 km s$^{-1}$. Moreover, as the maximum separation is $\sim$15 Mpc, they compose a supercluster at $z\sim$0.91, meaning that this is one of the most massive superclusters at this redshift to date. The galaxy density map implies that the confirmed clusters are embedded in a larger structure stretching over $\sim$100 Mpc. $\Lambda$CDM models predict about one supercluster like this in our surveyed volume, consistent with our finding so far. However, there are more supercluster candidates in this field, suggesting that additional studies are required to determine if the $\Lambda$CDM cosmological model can successfully reproduce the LSSs at high redshift.

Recent studies suggest that faint active galactic nuclei may be responsible for the reionization of the universe. Confirmation of this scenario requires spectroscopic identification of faint quasars ($M_{1450}>-24$ mag) at $z \gtrsim6$, but only a very small number of such quasars have been spectroscopically identified so far. Here, we report the discovery of a faint quasar IMS J220417.92+011144.8 at z~6 in a 12.5 deg$^{2}$ region of the SA22 field of the Infrared Medium-deep Survey (IMS). The spectrum of the quasar shows a sharp break at $\sim8443~\rm{\AA}$, with emission lines redshifted to $z=5.944 \pm 0.002$ and rest-frame ultraviolet continuum magnitude $M_{1450}=-23.59\pm0.10$ AB mag. The discovery of IMS J220417.92+011144.8 is consistent with the expected number of quasars at z~6 estimated from quasar luminosity functions based on previous observations of spectroscopically identified low-luminosity quasars . This suggests that the number of $M_{1450}\sim-23$ mag quasars at z~6 may not be high enough to fully account for the reionization of the universe. In addition, our study demonstrates that faint quasars in the early universe can be identified effectively with a moderately wide and deep near-infrared survey such as the IMS.

The optical light curve of the quasar PG 1302-102 at $z = 0.278$ shows a strong, smooth 5.2 yr periodic signal, detectable over a period of $\sim 20$ yr. Although the interpretation of this phenomenon is still uncertain, the most plausible mechanisms involve a binary system of two supermassive black holes with a subparsec separation. At this close separation, the nuclear black holes in PG 1302-102 will likely merge within $\sim 10^{5}$ yr due to gravitational wave emission alone. Here we report the rest-frame near-infrared time lags for PG 1302-102. Compiling data from {\it WISE} and {\it Akari}, we confirm that the periodic behavior reported in the optical light curve from Graham et al. (2015) is reproduced at infrared wavelengths, with best-fit observed-frame 3.4 and $4.6 \mu$m time lags of $(2219 \pm 153, 2408 \pm 148)$ days for a near face-on orientation of the torus, or $(4103\pm 153, 4292 \pm 148)$ days for an inclined system with relativistic Doppler boosting in effect. The periodicity in the infrared light curves and the light-travel time of the accretion disk photons to reach the dust glowing regions support that a source within the accretion disk is responsible for the optical variability of PG 1302-102, echoed at the further out dusty regions. The implied distance of this dusty, assumed toroidal region is $\sim$ 1.5 pc for a near face-on geometry, or $\sim$1.1 pc for the relativistic Doppler boosted case.

We present the rest-frame optical spectral properties of 155 luminous quasars at 3.3<z<6.4 taken with the AKARI space telescope, including the first detection of H$\alpha$ emission line as far out as z~6. We extend the scaling relation between the rest-frame optical continuum and line luminosity of active galactic nuclei (AGNs) to the high luminosity, high redshift regime that has rarely been probed before. Remarkably, we find that a single log-linear relation can be applied to the 5100${\rm \AA}$ and H$\alpha$ AGN luminosities over a wide range of luminosity (10$^{42}$<$L_{5100}$<10$^{47}$ergs/s) or redshift (0<z<6), suggesting that the physical mechanism governing this relation is unchanged from z=0 to 6, over five decades in luminosity. Similar scaling relations are found between the optical and the UV continuum luminosities or line widths. Applying the scaling relations to the H$\beta$ black hole mass ($M_{\rm BH}$) estimator of local AGNs, we derive the $M_{\rm BH}$ estimators based on H$\alpha$, MgII, and CIV lines, finding that the UV-line based masses are overall consistent with the Balmer-line based, but with a large intrinsic scatter of 0.40dex for the CIV estimates. Our 43 $M_{\rm BH}$ estimates from H$\alpha$ confirm the existence of BHs as massive as ~10$^{10}M_{\odot}$ out to z~5, and provide a secure footing for previous MgII-line based studies that a rapid $M_{\rm BH}$ growth has occurred in the early universe.

We present 2.5-5.0 $\mu$m spectra of 83 nearby ($0.002\,<\,z\,<\,0.48$) and bright ($K<14$mag) type-1 active galactic nuclei (AGNs) taken with the Infrared Camera (IRC) on board $\it{AKARI}$. The 2.5-5.0 $\mu$m spectral region contains emission lines such as Br$\beta$ (2.63 $\mu$m), Br$\alpha$ (4.05 $\mu$m), and polycyclic aromatic hydrocarbons (PAH; 3.3 $\mu$m), which can be used for studying the black hole (BH) masses and star formation activities in the host galaxies of AGNs. The spectral region also suffers less dust extinction than in the ultra violet (UV) or optical wavelengths, which may provide an unobscured view of dusty AGNs. Our sample is selected from bright quasar surveys of Palomar-Green (PG) and SNUQSO, and AGNs with reverberation-mapped BH masses from Peterson et al. (2004). Using 11 AGNs with reliable detection of Brackett lines, we derive the Brackett-line-based BH mass estimators. We also find that the observed Brackett line ratios can be explained with the commonly adopted physical conditions of the broad line region (BLR). Moreover, we fit the hot and warm dust components of the dust torus by adding photometric data of SDSS, 2MASS, $\it{WISE}$, and $\it{ISO}$ to the $\it{AKARI}$ spectra, finding hot and warm dust temperatures of $\sim1100\,\rm{K}$ and $\sim220\,\rm{K}$, respectively, rather than the commonly cited hot dust temperature of 1500 K.

Activity at the centers of galaxies, during which the central supermassive black hole is accreting material, is nowadays accepted to be rather ubiquitous and most probably a phase of every galaxy's evolution. It has been suggested that galactic mergers and interactions may be the culprits behind the triggering of nuclear activity. We use near-infrared data from the new Infrared Medium-Deep Survey (IMS) and the Deep eXtragalactic Survey (DXS) of the VIMOS-SA22 field and radio data at 1.4 GHz from the FIRST survey and a deep VLA survey to study the environments of radio-AGN over an area of ~25 sq. degrees and down to a radio flux limit of 0.1 mJy and a J-band magnitude of 23 mag AB. Radio-AGN are predominantly found in environments similar to those of control galaxies at similar redshift, J-band magnitude, and U-R rest-frame absolute color. However, a sub-population of radio-AGN is found in environments up to 100 times denser than their control sources. We thus preclude merging as the dominant triggering mechanism of radio-AGN. Through the fitting of the broadband spectral energy distribution of radio-AGN in the least and most dense environments, we find that those in the least dense environments show higher radio-loudness, higher star formation efficiencies, and higher accretion rates, typical of the so-called high-excitation radio-AGN. These differences tend to disappear at z>1. We interpret our results in terms of a different triggering mechanism for these sources that is driven by mass-loss through winds of young stars created during the observed ongoing star formation.

We identify and characterize a population of luminous dust poor quasars at 0<z<5, similar in photometric properties to the objects found at z>6 previously. This class of active galactic nuclei has been known to show little IR emission from a dusty structure, but is yet poorly understood in terms of number evolution or of dependence on physical quantities. In order to better understand the luminous dust poor quasar properties, we compiled a rest-frame UV to IR library of 41,000 optically selected type-1 quasars with $L_{bol}>10^{45.7} erg s^{-1}$. After fitting the broad-band spectral energy distributions (SEDs) with accretion disk and dust components, we find 0.6% of our sample to be hot dust poor with a rest-frame 2.3 micron to 0.51 micron flux density ratio of -0.5 dex or less. The dust poor SEDs are blue in the UV-optical and weak in the MIR, such that their accretion disks are less obscured, and that hot dust emission traces that of warm dust down to the dust poor regime. At a given bolometric luminosity, dust poor quasars are lower in black hole mass and higher in Eddington ratio than general luminous quasars, suggesting that they are in a rapidly growing evolutionary state in which the dust poor phase appears as a short or rare phenomenon. The dust poor fraction increases with redshift, and possible implications for the evolution of the dust poor fraction are discussed.

Polycyclic Aromatic Hydrocarbon (PAH) emission features dominate the mid-infrared spectra of star-forming galaxies and can be useful to calibrate star formation rates and diagnose ionized states of grains. However, the PAH 3.3 micron feature has not been studied as much as other PAH features since it is weaker than others and resides outside of Spitzer capability. In order to detect and calibrate the 3.3 micron PAH emission and investigate its potential as a star formation rate indicator, we carried out an AKARI mission program, AKARI mJy Unbiased Survey of Extragalactic Survey (AMUSES) and compare its sample with various literature samples. We obtained 2 ~5 micron low resolution spectra of 20 flux-limited galaxies with mixed SED classes, which yields the detection of the 3.3 micron PAH emission from three out of 20 galaxies. For the combined sample of AMUSES and literature samples, the 3.3 micron PAH luminosities correlate with the infrared luminosities of star-forming galaxies, albeit with a large scatter (1.5 dex). The correlation appears to break down at the domain of ultra-luminous infrared galaxies (ULIRGs), and the power of the 3.3 micron PAH luminosity as a proxy for the infrared luminosity is hampered at log[L(PAH3.3)/(erg/sec)] > -42.0. Possible origins for this deviation in the correlation are discussed, including contribution from AGN and strongly obscured YSOs, and the destruction of PAH molecules in ULIRGs.

Three observables of early-type galaxies - size ($r_{e}$), surface brightness ($I_{e}$), and velocity dispersion ($\sigma_{0}$) - form a tight planar correlation known as the fundamental plane (FP), which has provided great insights into the galaxy formation and the evolution processes. However, the FP has been found to be tilted against the simple virial expectation, prompting debates on its origin. In order to investigate the contribution of systematic stellar population variation to the FP tilt, we study here the FP relations of early-type galaxies in mid-infrared (MIR) which may represent the stellar mass well. We examined the wavelength dependence of the FP coefficients, $a$ and $b$ in $\log r_{e}= a\log\sigma_{0} + b\log< I >_{e} + c$, using a sample of 56 early-type galaxies for which visible (V-band), near-infrared (K-band), and MIR (Spitzer IRAC, 3.6--8.0$\mu$m) data are available. We find that the coefficient $a$ increases as a function of wavelength as $da/d\lambda=0.11\pm0.04\mu m^{-1}$, while the coefficient $b$ reaches the closest to -1 at 3.6--5.8$\mu$m. When applied to the visible FP coefficients derived from a larger sample of nearby early-type galaxies, we get the FP relation with $(a,b) \simeq $(1.6--1.8,-0.9) at 3.6$\mu$m. Our result suggests that the stellar population effect can explain more than half of the FP tilt, closing the gap between the virial expectation and the optical FP. The reduction in the FP tilt is reflected in the dynamical mass-to-light ratio, $M_{dyn}/L$, dependence on $L$ which decreases toward 3.6--5.8$\mu$m, suggesting that the MIR light better represents mass than the shorter wavelengths.