Source author record

Yue Shen

Yue Shen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

112works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

112 published item(s)

preprint2026arXiv

The Advanced X-ray Imaging Satellite (AXIS) Community Science Book

The AXIS Community Science Book represents the collective effort of 592 scientists worldwide to define the transformative science enabled by the Advanced X-ray Imaging Satellite (AXIS), a next-generation X-ray mission selected by NASA's Astrophysics Probe Program for Phase A study. AXIS will advance the legacy of high-angular-resolution X-ray astronomy with ~1.5'' imaging over a wide 24' field of view and an order of magnitude greater collecting area than Chandra in the 0.3-12 keV band. Combining sharp imaging, high throughput, and rapid response capabilities, AXIS will open new windows on virtually every aspect of modern astrophysics, exploring the birth and growth of supermassive black holes, the feedback processes that shape galaxies, the life cycles of stars and exoplanet environments, and the nature of compact stellar remnants, supernova remnants, and explosive transients. This book compiles 138 community-contributed science cases developed by five Science Working Groups focused on AGN and supermassive black holes, galaxy evolution and feedback, compact objects and supernova remnants, stellar physics and exoplanets, and time-domain and multi-messenger astrophysics. Together, these studies establish the scientific foundation for next-generation X-ray exploration in the 2030s and highlight strong synergies with facilities of the 2030s, such as JWST, Roman, Rubin/LSST, SKA, ALMA, ngVLA, and next-generation gravitational-wave and neutrino networks.

preprint2024arXiv

From Beginner to Expert: Modeling Medical Knowledge into General LLMs

Recently, large language model (LLM) based artificial intelligence (AI) systems have demonstrated remarkable capabilities in natural language understanding and generation. However, these models face a significant challenge when it comes to sensitive applications, such as reasoning over medical knowledge and answering medical questions in a physician-like manner. Prior studies attempted to overcome this challenge by increasing the model size (>100B) to learn more general medical knowledge, while there is still room for improvement in LLMs with smaller-scale model sizes (<100B). In this work, we start from a pre-trained general LLM model (AntGLM-10B) and fine-tune it from a medical beginner towards a medical expert (called AntGLM-Med-10B), which leverages a 3-stage optimization procedure, i.e., general medical knowledge injection, medical domain instruction tuning, and specific medical task adaptation. Our contributions are threefold: (1) We specifically investigate how to adapt a pre-trained general LLM in medical domain, especially for a specific medical task. (2) We collect and construct large-scale medical datasets for each stage of the optimization process. These datasets encompass various data types and tasks, such as question-answering, medical reasoning, multi-choice questions, and medical conversations. (3) Specifically for multi-choice questions in the medical domain, we propose a novel Verification-of-Choice approach for prompting engineering, which significantly enhances the reasoning ability of LLMs. Remarkably, by combining the above approaches, our AntGLM-Med-10B model can outperform the most of LLMs on PubMedQA, including both general and medical LLMs, even when these LLMs have larger model size.

preprint2024arXiv

RJUA-QA: A Comprehensive QA Dataset for Urology

We introduce RJUA-QA, a novel medical dataset for question answering (QA) and reasoning with clinical evidence, contributing to bridge the gap between general large language models (LLMs) and medical-specific LLM applications. RJUA-QA is derived from realistic clinical scenarios and aims to facilitate LLMs in generating reliable diagnostic and advice. The dataset contains 2,132 curated Question-Context-Answer pairs, corresponding about 25,000 diagnostic records and clinical cases. The dataset covers 67 common urological disease categories, where the disease coverage exceeds 97.6\% of the population seeking medical services in urology. Each data instance in RJUA-QA comprises: (1) a question mirroring real patient to inquiry about clinical symptoms and medical conditions, (2) a context including comprehensive expert knowledge, serving as a reference for medical examination and diagnosis, (3) a doctor response offering the diagnostic conclusion and suggested examination guidance, (4) a diagnosed clinical disease as the recommended diagnostic outcome, and (5) clinical advice providing recommendations for medical examination. RJUA-QA is the first medical QA dataset for clinical reasoning over the patient inquiries, where expert-level knowledge and experience are required for yielding diagnostic conclusions and medical examination advice. A comprehensive evaluation is conducted to evaluate the performance of both medical-specific and general LLMs on the RJUA-QA dataset. Our data is are publicly available at \url{https://github.com/alipay/RJU_Ant_QA}.

preprint2023arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: The Black Hole Mass$-$Stellar Mass Relations at $0.2\lesssim z\lesssim 0.8$

We measure the correlation between black-hole mass $M_{\rm BH}$ and host stellar mass $M_*$ for a sample of 38 broad-line quasars at $0.2\lesssim z\lesssim 0.8$ (median redshift $z_{\rm med}=0.5$). The black-hole masses are derived from a dedicated reverberation mapping program for distant quasars, and the stellar masses are estimated from two-band optical+IR HST imaging. Most of these quasars are well centered within $\lesssim 1$kpc from the host galaxy centroid, with only a few cases in merging/disturbed systems showing larger spatial offsets. Our sample spans two orders of magnitude in stellar mass ($\sim 10^9-10^{11}\,M_\odot$) and black-hole mass ($\sim 10^7-10^9\,M_\odot$), and reveals a significant correlation between the two quantities. We find a best-fit intrinsic (i.e., selection effects corrected) $M_{\rm BH}-M_{\rm *,host}$ relation of $\log (M_{\rm BH}/M_{\rm \odot})=7.01_{-0.33}^{+0.23} + 1.74_{-0.64}^{+0.64}\log (M_{\rm *,host}/10^{10}M_{\rm \odot})$, with an intrinsic scatter of $0.47_{-0.17}^{+0.24}$dex. Decomposing our quasar hosts into bulges and disks, there is a similar $M_{\rm BH}-M_{\rm *,bulge}$ relation with a slightly larger scatter, likely caused by systematic uncertainties in the bulge-disk decomposition. The $M_{\rm BH}-M_{\rm *,host}$ relation at $z_{\rm med}=0.5$ is similar to that in local quiescent galaxies, with negligible evolution over the redshift range probed by our sample. With direct black-hole masses from reverberation mapping and a large dynamical range of the sample, selection biases do not appear to affect our conclusions significantly. Our results, along with other samples in the literature, suggest that the locally-measured black-hole mass$-$host stellar mass relation is already in place at $z\sim 1$.

preprint2022arXiv

Binary Neural Networks as a general-propose compute paradigm for on-device computer vision

For binary neural networks (BNNs) to become the mainstream on-device computer vision algorithm, they must achieve a superior speed-vs-accuracy tradeoff than 8-bit quantization and establish a similar degree of general applicability in vision tasks. To this end, we propose a BNN framework comprising 1) a minimalistic inference scheme for hardware-friendliness, 2) an over-parameterized training scheme for high accuracy, and 3) a simple procedure to adapt to different vision tasks. The resultant framework overtakes 8-bit quantization in the speed-vs-accuracy tradeoff for classification, detection, segmentation, super-resolution and matching: our BNNs not only retain the accuracy levels of their 8-bit baselines but also showcase 1.3-2.4$\times$ faster FPS on mobile CPUs. Similar conclusions can be drawn for prototypical systolic-array-based AI accelerators, where our BNNs promise 2.8-7$\times$ fewer execution cycles than 8-bit and 2.1-2.7$\times$ fewer cycles than alternative BNN designs. These results suggest that the time for large-scale BNN adoption could be upon us.

preprint2022arXiv

Dwarf AGNs from Optical Variability for the Origins of Seeds (DAVOS): Insights from the Dark Energy Survey Deep Fields

We present a sample of 706, $z < 1.5$ active galactic nuclei (AGNs) selected from optical photometric variability in three of the Dark Energy Survey (DES) deep fields (E2, C3, and X3) over an area of 4.64 deg$^2$. We construct light curves using difference imaging aperture photometry for resolved sources and non-difference imaging PSF photometry for unresolved sources, respectively, and characterize the variability significance. Our DES light curves have a mean cadence of 7 days, a 6 year baseline, and a single-epoch imaging depth of up to $g \sim 24.5$. Using spectral energy distribution (SED) fitting, we find 26 out of total 706 variable galaxies are consistent with dwarf galaxies with a reliable stellar mass estimate ($M_{\ast}<10^{9.5}\ M_\odot$; median photometric redshift of 0.9). We were able to constrain rapid characteristic variability timescales ($\sim$ weeks) using the DES light curves in 15 dwarf AGN candidates (a subset of our variable AGN candidates) at a median photometric redshift of 0.4. This rapid variability is consistent with their low black hole masses. We confirm the low-mass AGN nature of one source with a high S/N optical spectrum. We publish our catalog, optical light curves, and supplementary data, such as X-ray properties and optical spectra, when available. We measure a variable AGN fraction versus stellar mass and compare to results from a forward model. This work demonstrates the feasibility of optical variability to identify AGNs with lower black hole masses in deep fields, which may be more "pristine" analogs of supermassive black hole seeds.

preprint2022arXiv

Optical Variability of Quasars with 20-Year Photometric Light Curves

We study the optical $gri$ photometric variability of a sample of 190 quasars within the SDSS Stripe 82 region that have long-term photometric coverage during $\sim 1998-2020$ with SDSS, PanSTARRS-1, the Dark Energy Survey, and dedicated follow-up monitoring with Blanco 4m/DECam. With on average $\sim 200$ nightly epochs per quasar per filter band, we improve the parameter constraints from a Damped Random Walk (DRW) model fit to the light curves over previous studies with 10-15 yr baselines and $\lesssim 100$ epochs. We find that the average damping timescale $τ_{\rm DRW}$ continues to rise with increased baseline, reaching a median value of $\sim 750$ days ($g$ band) in the rest-frame of these quasars using the 20-yr light curves. Some quasars may have gradual, long-term trends in their light curves, suggesting that either the DRW fit requires very long baselines to converge, or that the underlying variability is more complex than a single DRW process for these quasars. Using a subset of quasars with better-constrained $τ_{\rm DRW}$ (less than 20\% of the baseline), we confirm a weak wavelength dependence of $τ_{\rm DRW}\propto λ^{0.51\pm0.20}$. We further quantify optical variability of these quasars over days to decades timescales using structure function (SF) and power spectrum density (PSD) analyses. The SF and PSD measurements qualitatively confirm the measured (hundreds of days) damping timescales from the DRW fits. However, the ensemble PSD is steeper than that of a DRW on timescales less than $\sim$ a month for these luminous quasars, and this second break point correlates with the longer DRW damping timescale.

preprint2022arXiv

The eROSITA Final Equatorial-Depth Survey (eFEDS): Galaxy Clusters and Groups in Disguise

The eROSITA Final Equatorial-Depth Survey (eFEDS), executed during the performance verification phase of the Spectrum-Roentgen-Gamma (SRG)/eROSITA telescope, was completed in Nov. 2019. One of the science goals of this survey is to demonstrate the ability of eROSITA to detect samples of clusters and groups at the final depth of the eROSITA all-sky survey. Because of the sizeable point-spread function of eROSITA, high-redshift clusters of galaxies or compact nearby groups hosting bright active galactic nuclei (AGN) can be misclassified as point sources by the source detection algorithms. A total of 346 galaxy clusters and groups in the redshift range of 0.1<z<1.3 were identified based on their red sequence in the point source catalog. We examine the multiwavelength properties of these clusters and groups to understand the potential biases in our selection process and the completeness of the extent-selected sample. The majority of the clusters and groups in the point source sample are indeed underluminous and compact compared to the extent-selected sample. Their faint X-ray emission, well below the flux limit of the extent-selected eFEDS clusters, and their compact X-ray emission are likely to be the main reason for this misclassification. In the sample, we confirm that 10% of the sources host AGN in their brightest cluster galaxies (BCGs) through optical spectroscopy and visual inspection. By studying their X-ray, optical, infrared, and radio properties, we establish a method for identifying clusters and groups that host AGN in their BCGs. We successfully test this method on the current point source catalog through the Sloan Digital Sky Survey optical spectroscopy and find eight clusters and groups with active radio-loud AGN that are particularly bright in the infrared. They include eFEDSJ091437.8+024558, eFEDSJ083520.1+012516, and eFEDSJ092227.1+043339 at redshifts 0.3-0.4. [ABRIDGED]

preprint2022arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: UV-Optical Accretion Disk Measurements with Hubble Space Telescope

We present accretion-disk structure measurements from UV-optical reverberation mapping observations of a sample of eight quasars at 0.24<z<0.85. Ultraviolet photometry comes from two cycles of Hubble Space Telescope monitoring, accompanied by multi-band optical monitoring by the Las Cumbres Observatory network and Liverpool Telescopes. The targets were selected from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project sample with reliable black-hole mass measurements from Hbeta reverberation mapping results. We measure significant lags between the UV and various optical griz bands using JAVELIN and CREAM methods. We use the significant lag results from both methods to fit the accretion-disk structure using a Markov chain Monte Carlo approach. We study the accretion disk as a function of disk normalization, temperature scaling, and efficiency. We find direct evidence for diffuse nebular emission from Balmer and FeII lines over discrete wavelength ranges. We also find that our best-fit disk color profile is broadly consistent with the Shakura \& Sunyaev disk model. We compare our UV-optical lags to the disk sizes inferred from optical-optical lags of the same quasars and find that our results are consistent with these quasars being drawn from a limited high-lag subset of the broader population. Our results are therefore broadly consistent with models that suggest longer disk lags in a subset of quasars, for example, due to a nonzero size of the ionizing corona and/or magnetic heating contributing to the disk response.

preprint2021arXiv

Metallicity in Quasar Broad Line Regions at Redshift $\sim$ 6

Broad line regions (BLRs) in high-redshift quasars provide crucial information of chemical enrichment in the early universe. Here we present a study of BLR metallicities in 33 quasars at redshift $5.7<z<6.4$. Using the near-IR spectra of the quasars obtained from the Gemini telescope, we measure their rest-frame UV emission line flux and calculate flux ratios. We then estimate BLR metallicities with empirical calibrations based on photoionization models. The inferred median metallicity of our sample is a few times the solar value, indicating that the BLR gas had been highly metal-enriched at $z\sim6$. We compare our sample with a low-redshift quasar sample with similar luminosities and find no evidence of redshift evolution in quasar BLR metallicities. This is consistent with previous studies. The Fe II$/$Mg II flux ratio, a proxy for the Fe$/α$ element abundance ratio, shows no redshift evolution as well, further supporting rapid nuclear star formation at $z\sim6$. We also find that the black hole mass-BLR metallicity relation at $z\sim6$ is consistent with the relation measured at $2<z<5$, suggesting that our results are not biased by a selection effect due to this relation.

preprint2021arXiv

The eROSITA Final Equatorial-Depth Survey (eFEDS): The AGN Catalogue and its X-ray Spectral Properties

Context: After the successful launch of the Spectrum-Roentgen-Gamma (SRG) mission in July 2019, eROSITA, the soft X-ray instrument aboard SRG, performed scanning observations of a large contiguous field, namely the eROSITA Final Equatorial Depth Survey (eFEDS), ahead of the planned four-year all-sky survey. eFEDS yielded a large sample of X-ray sources with very rich multi-band photometric and spectroscopic coverage. Aims: We present here the eFEDS Active Galactic Nuclei (AGN) catalog and the eROSITA X-ray spectral properties of the eFEDS sources. Methods: Using a Bayesian method, we perform a systematic X-ray spectral analysis for all eFEDS sources. The appropriate model is chosen based on the source classification and the spectral quality, and, in the case of AGN, including the possibility of intrinsic (rest-frame) absorption and/or soft excess emission. Hierarchical Bayesian modeling (HBM) is used to estimate the spectral parameter distribution of the sample. Results: X-ray spectral properties are presented for all eFEDS X-ray sources. There are 21952 candidate AGN, which comprise 79% of the eFEDS sample. Despite a large number of faint sources with low photon counts, our spectral fitting provides meaningful measurements of fluxes, luminosities, and spectral shapes for a majority of the sources. This AGN catalog is dominated by X-ray unobscured sources, with an obscured (logNH>21.5) fraction of 10% derived by HBM. The power-law slope of the catalog can be described by a Gaussian distribution of 1.94+-0.22. Above a photon counts threshold of 500, nine out of 50 AGN have soft excess detected. For the sources with blue UV to optical color (type-I AGN), the X-ray emission is well correlated with the UV emission with the usual anti-correlation between the X-ray to UV spectral slope α_{OX} and the UV luminosity.

preprint2021arXiv

Varstrometry for Off-nucleus and Dual sub-Kpc AGN (VODKA): Hubble Space Telescope Discovers Double Quasars

Dual supermassive black holes (SMBHs) at $\sim$kpc scales are the progenitor population of SMBH mergers and play an important role in understanding the pairing and dynamical evolution of massive black holes in galaxy mergers. Because of the stringent resolution requirement and the apparent rareness of these small-separation pairs, there are scarce observational constraints on this population, with few confirmed dual SMBHs at $<10$kpc separations at $z>1$. Here we present results from a pilot search for kpc-scale dual quasars selected with Gaia Data release 2 (DR2) astrometry and followed up with Hubble Space Telescope (HST) Wide Field Camera 3 dual-band (F475W and F814W) snapshot imaging. Our targets are quasars primarily selected with the varstrometry technique, i.e., light centroid jitter caused by asynchronous variability from both members in an unresolved quasar pair, supplemented by sub-arcsec pairs already resolved by Gaia DR2. We find an overall high fraction of HST-resolved pairs among the varstrometry-selected quasars (unresolved in Gaia DR2), $\sim 30-50\%$, increasing toward high redshift ($\sim 60-80\%$ at $z>1.5$). We discuss the nature of the 43 resolved sub-arcsec pairs based on HST and supplementary data. A substantial fraction ($\sim 40\%$) of these pairs are likely physical quasar pairs or gravitationally lensed quasars. We also discover a triple quasar candidate and a quadruply lensed quasar, which is among the smallest-separation quadruple lenses. These results provide important guidelines to improve varstrometry selection and follow-up confirmation of $\sim$kpc-scale dual SMBHs at high redshift.

preprint2020arXiv

Characterization of Optical Light Curves of Extreme Variability Quasars Over a ~16-yr Baseline

We study the optical light curves - primarily probing the variable emission from the accretion disk - of ~ 900 extreme variability quasars (EVQs, with maximum flux variations more than 1 mag) over an observed-frame baseline of ~ 16 years using public data from the SDSS Stripe 82, PanSTARRS-1 and the Dark Energy Survey. We classify the multi-year long-term light curves of EVQs into three categories roughly in the order of decreasing smoothness: monotonic decreasing or increasing (3.7%), single broad peak and dip (56.8%), and more complex patterns (39.5%). The rareness of monotonic cases suggests that the major mechanisms driving the extreme optical variability do not operate over timescales much longer than a few years. Simulated light curves with a damped random walk model generally under-predict the first two categories with smoother long-term trends. Despite the different long-term behaviors of these EVQs, there is little dependence of the long-term trend on the physical properties of quasars, such as their luminosity, BH mass, and Eddington ratio. The large dynamic range of optical flux variability over multi-year timescales of these EVQs allows us to explore the ensemble correlation between the short-term (< 6 months) variability and the seasonal-average flux across the decade-long baseline (the rms-mean flux relation). We find that unlike the results for X-ray variability studies, the linear short-term flux variations do not scale with the seasonal-average flux, indicating different mechanisms that drive the short-term flickering and long-term extreme variability of accretion disk emission. Finally, we present a sample of 16 EVQs, where the approximately bell-shaped large amplitude variation in the light curve can be reasonably well fit by a simple microlensing model.

preprint2020arXiv

Dark Energy Survey Identification of A Low-Mass Active Galactic Nucleus at Redshift 0.823 from Optical Variability

We report the identification of a low-mass AGN, DES J0218$-$0430, in a redshift $z = 0.823$ galaxy in the Dark Energy Survey (DES) Supernova field. We select DES J0218$-$0430 as an AGN candidate by characterizing its long-term optical variability alone based on DES optical broad-band light curves spanning over 6 years. An archival optical spectrum from the fourth phase of the Sloan Digital Sky Survey shows both broad Mg II and broad H$β$ lines, confirming its nature as a broad-line AGN. Archival XMM-Newton X-ray observations suggest an intrinsic hard X-ray luminosity of $L_{\rm 2-12\,keV}\sim7.6\pm0.4\times10^{43}$ erg s$^{-1}$, which exceeds those of the most X-ray luminous starburst galaxies, in support of an AGN driving the optical variability. Based on the broad H$β$ from SDSS spectrum, we estimate a virial BH mass of $M_{\bullet}\approx10^{6.43}$-$10^{6.72}M_{\odot}$ (with the error denoting 1$σ$ statistical uncertainties only), consistent with the estimation from OzDES, making it the lowest mass AGN with redshift $>$ 0.4 detected in optical. We estimate the host galaxy stellar mass to be $M_{\ast}\sim10^{10.5\pm0.3}M_{\odot}$ based on modeling the multi-wavelength spectral energy distribution. DES J0218$-$0430 extends the $M_{\bullet}$-$M_{\ast}$ relation observed in luminous AGNs at $z\sim1$ to masses lower than being probed by previous work. Our work demonstrates the feasibility of using optical variability to identify low-mass AGNs at higher redshift in deeper synoptic surveys with direct implications for the upcoming Legacy Survey of Space and Time at Vera C. Rubin Observatory.

preprint2020arXiv

Dust Reverberation Mapping in Distant Quasars from Optical and Mid-Infrared Imaging Surveys

The size of the dust torus in Active Galactic Nuclei (AGN) and their high-luminosity counterparts, quasars, can be inferred from the time delay between UV/optical accretion disk continuum variability and the response in the mid-infrared (MIR) torus emission. This dust reverberation mapping (RM) technique has been successfully applied to $\sim 70$ $z\lesssim 0.3$ AGN and quasars. Here we present first results of our dust RM program for distant quasars covered in the SDSS Stripe 82 region combining $\sim 20$-yr ground-based optical light curves with 10-yr MIR light curves from the WISE satellite. We measure a high-fidelity lag between W1-band (3.4 $μ$m) and $g$ band for 587 quasars over $0.3\lesssim z\lesssim 2$ ($\left<z\right>\sim 0.8$) and two orders of magnitude in quasar luminosity. They tightly follow (intrinsic scatter $\sim 0.17$ dex in lag) the IR lag-luminosity relation observed for $z<0.3$ AGN, revealing a remarkable size-luminosity relation for the dust torus over more than four decades in AGN luminosity, with little dependence on additional quasar properties such as Eddington ratio and variability amplitude. This study motivates further investigations in the utility of dust RM for cosmology, and strongly endorses a compelling science case for the combined 10-yr Vera C. Rubin Observatory Legacy Survey of Space and Time (optical) and 5-yr Nancy Grace Roman Space Telescope 2$μ$m light curves in a deep survey for low-redshift AGN dust RM with much lower luminosities and shorter, measurable IR lags. The compiled optical and MIR light curves for 7,384 quasars in our parent sample are made public with this work.

preprint2020arXiv

Graph Representation Learning for Merchant Incentive Optimization in Mobile Payment Marketing

Mobile payment such as Alipay has been widely used in our daily lives. To further promote the mobile payment activities, it is important to run marketing campaigns under a limited budget by providing incentives such as coupons, commissions to merchants. As a result, incentive optimization is the key to maximizing the commercial objective of the marketing campaign. With the analyses of online experiments, we found that the transaction network can subtly describe the similarity of merchants' responses to different incentives, which is of great use in the incentive optimization problem. In this paper, we present a graph representation learning method atop of transaction networks for merchant incentive optimization in mobile payment marketing. With limited samples collected from online experiments, our end-to-end method first learns merchant representations based on an attributed transaction networks, then effectively models the correlations between the commercial objectives each merchant may achieve and the incentives under varying treatments. Thus we are able to model the sensitivity to incentive for each merchant, and spend the most budgets on those merchants that show strong sensitivities in the marketing campaign. Extensive offline and online experimental results at Alipay demonstrate the effectiveness of our proposed approach.

preprint2020arXiv

High-Order Paired-ASPP Networks for Semantic Segmenation

Current semantic segmentation models only exploit first-order statistics, while rarely exploring high-order statistics. However, common first-order statistics are insufficient to support a solid unanimous representation. In this paper, we propose High-Order Paired-ASPP Network to exploit high-order statistics from various feature levels. The network first introduces a High-Order Representation module to extract the contextual high-order information from all stages of the backbone. They can provide more semantic clues and discriminative information than the first-order ones. Besides, a Paired-ASPP module is proposed to embed high-order statistics of the early stages into the last stage. It can further preserve the boundary-related and spatial context in the low-level features for final prediction. Our experiments show that the high-order statistics significantly boost the performance on confusing objects. Our method achieves competitive performance without bells and whistles on three benchmarks, i.e, Cityscapes, ADE20K and Pascal-Context with the mIoU of 81.6%, 45.3% and 52.9%.

preprint2020arXiv

Optical Variability of the Dwarf AGN NGC 4395 from the Transiting Exoplanet Survey Satellite

We present optical light curves from the Transiting Exoplanet Survey Satellite (TESS) for the archetypical dwarf active galactic nucleus (AGN) in the nearby galaxy NGC 4395 hosting a $\sim 10^5\,M_\odot$ supermassive black hole (SMBH). Significant variability is detected on timescales from weeks to hours before reaching the background noise level. The $\sim$month-long, 30 minute-cadence, high-precision TESS light curve can be well fit by a simple damped random walk (DRW) model, with the damping timescale $τ_{\rm DRW}$ constrained to be $2.3_{-0.7}^{+1.8}$~days ($1σ$). NGC 4395 lies almost exactly on the extrapolation of the $τ_{\rm DRW}-M_{\rm BH}$ relation measured for AGNs with BH masses that are more than three orders of magnitude larger. The optical variability periodogram can be well fit by a broken power law with the high-frequency slope ($-1.88\pm0.15$) and the characteristic timescale ($τ_{\rm br}\equiv 1/(2πf_{\rm br})=1.4_{-0.5}^{+1.9}\,$days) consistent with the DRW model within 1$σ$. This work demonstrates the power of TESS light curves in identifying low-mass accreting SMBHs with optical variability, and a potential global $τ_{\rm DRW}-M_{\rm BH}$ relation that can be used to estimate SMBH masses with optical variability measurements.

preprint2020arXiv

Placing High-Redshift Quasars in Perspective: a Catalog of Spectroscopic Properties from the Gemini Near Infrared Spectrograph -- Distant Quasar Survey

We present spectroscopic measurements for 226 sources from the Gemini Near Infrared Spectrograph - Distant Quasar Survey (GNIRS-DQS). Being the largest uniform, homogeneous survey of its kind, it represents a flux-limited sample ($m_{i}$ ${\lesssim}$ 19.0 mag, $H$ ${\lesssim}$ 16.5 mag) of Sloan Digital Sky Survey (SDSS) quasars at 1.5 ${\lesssim}$ $z$ ${\lesssim}$ 3.5 with a monochromatic luminosity ($λL_λ$) at 5100 ${\unicode{xC5}}$ in the range of $10^{44} - 10^{46}$ erg $\rm{s}^{-1}$. A combination of the GNIRS and SDSS spectra covers principal quasar diagnostic features, chiefly the C IV $λ$1549, Mg II $λλ$2798, 2803, H$β$ $λ$4861, and [O III] $λλ$4959, 5007 emission lines, in each source. The spectral inventory will be utilized primarily to develop prescriptions for obtaining more accurate and precise redshifts, black hole masses, and accretion rates for all quasars. Additionally, the measurements will facilitate an understanding of the dependence of rest-frame ultraviolet-optical spectral properties of quasars on redshift, luminosity, and Eddington ratio, and test whether the physical properties of the quasar central engine evolve over cosmic time.

preprint2020arXiv

Strong Mg II and Fe II Absorbers at 2.2 < z < 6.0

We present a study of strong intervening absorption systems in the near-IR spectra of 31 luminous quasars at $z>5.7$. The quasar spectra were obtained with {\it Gemini} GNIRS that provide continuous wavelength coverage from $\sim$0.9 to $\sim$2.5 $μ$m. We detect 32 strong Mg II doublet absorbers with rest-frame equivalent width $W_r$ ($\lambda2796$) $>1.0$ Åat $2.2 < z < 6.0$. Each Mg II absorber is confirmed by at least two associated Fe II absorption lines in the rest-frame wavelength range of $\sim 1600-2600$ Å. We find that the comoving line density ($dN/dX$) of the strong Fe II-bearing Mg II absorbers decreases towards higher redshift at $z>3$, consistent with previous studies. Compared with strong Mg II absorbers detected in damped Ly$α$ systems at 2 $<z<$ 4, our absorbers are potentially less saturated and show much larger rest-frame velocity widths. This suggests that the gas traced by our absorbers are potentially affected by galactic superwinds. We analyze the {\it Hubble Space Telescope} near-IR images of the quasars and identify possible associated galaxies for our strong absorbers. There are a maximum of two galaxy candidates found within 5" radius of each absorber. The median F105W-band magnitude of these galaxy candidates is 24.8 mag, which is fainter than the $L^*$ galaxy luminosity at $z\sim$ 4. By using our observed $dN/dX$ of strong Mg II absorbers and galaxy candidates median luminosity, we suggest that at high redshift, strong Mg II absorbers tend to have a more disturbed environment but smaller halo size than that at $z <$ 1.

preprint2020arXiv

The Curious Case of PHL 293B: A Long-Lived Transient in a Metal-Poor Blue Compact Dwarf Galaxy

We report on small-amplitude optical variability and recent dissipation of the unusually persistent broad emission lines in the blue compact dwarf galaxy PHL 293B. The galaxy's unusual spectral features (P Cygni-like profiles with $\sim$800 km s$^{-1}$ blueshifted absorption lines) have resulted in conflicting interpretations of the nature of this source in the literature. However, analysis of new Gemini spectroscopy reveals the broad emission has begun to fade after being persistent for over a decade prior. Precise difference imaging light curves constructed with the Sloan Digital Sky Survey and the Dark Energy Survey reveal small-amplitude optical variability of $\sim$0.1 mag in the g band offset by $100\pm21$ pc from the brightest pixel of the host. The light curve is well-described by an active galactic nuclei (AGN)-like damped random walk process. However, we conclude that the origin of the optical variability and spectral features of PHL 293B is due to a long-lived stellar transient, likely a Type IIn supernova or non-terminal outburst, mimicking long-term AGN-like variability. This work highlights the challenges of discriminating between scenarios in such extreme environments, relevant to searches for AGNs in dwarf galaxies. This is the second long-lived transient discovered in a blue compact dwarf, after SDSS1133. Our result implies such long-lived stellar transients may be more common in metal-deficient galaxies. Systematic searches for low-level variability in dwarf galaxies will be possible with the upcoming Legacy Survey of Space and Time at Vera C. Rubin Observatory.

preprint2020arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Photometric g and i Light Curves

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.

preprint2020arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: The H$β$ Radius-Luminosity Relation

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.

preprint2020arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: The M_BH-Host Relations at 0.2<z<0.6 from Reverberation Mapping and Hubble Space Telescope Imaging

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.

preprint2020arXiv

TNT: Target-driveN Trajectory Prediction

Predicting the future behavior of moving agents is essential for real world applications. It is challenging as the intent of the agent and the corresponding behavior is unknown and intrinsically multimodal. Our key insight is that for prediction within a moderate time horizon, the future modes can be effectively captured by a set of target states. This leads to our target-driven trajectory prediction (TNT) framework. TNT has three stages which are trained end-to-end. It first predicts an agent's potential target states $T$ steps into the future, by encoding its interactions with the environment and the other agents. TNT then generates trajectory state sequences conditioned on targets. A final stage estimates trajectory likelihoods and a final compact set of trajectory predictions is selected. This is in contrast to previous work which models agent intents as latent variables, and relies on test-time sampling to generate diverse trajectories. We benchmark TNT on trajectory prediction of vehicles and pedestrians, where we outperform state-of-the-art on Argoverse Forecasting, INTERACTION, Stanford Drone and an in-house Pedestrian-at-Intersection dataset.

preprint2019arXiv

An Unusual Mid-Infrared Flare in a Type 2 AGN: An Obscured Turning-on AGN or Tidal Disruption Event?

We report the discovery of an exceptional MIR flare in a Type 2 AGN, SDSS J165726.81+234528.1, at $z=0.059$. This object brightened by 3 mag in the Wide-field Infrared Survey Explorer (WISE) $W1$ and $W2$ bands between 2015 and 2017 (and is fading since 2018), without significant changes ($\lesssim$ 0.2 mag) in the optical over the same period of time. Based on the WISE light curves and near-IR imaging, the flare is more significant at longer wavelengths, suggesting an origin of hot dust emission. The estimated black hole mass ($\sim 10^{6.5}\,M_\odot$) from different methods places its peak bolometric luminosity around the Eddington limit. The high luminosity of the MIR flare and its multi-year timescale suggest that it most likely originated from reprocessed dust radiation in an extended torus surrounding the AGN, instead of from stellar explosions. The MIR color variability is consistent with known changing-look AGN and tidal disruption events (TDEs), but inconsistent with normal supernovae. We suggest that it is a turning-on Type 2 AGN or TDE, where the optical variability is obscured by the dust torus during the transition. This MIR flare event reveals a population of dramatic nuclear transients that are missed in the optical.

preprint2019arXiv

Spectral Variability of a Sample of Extreme Variability Quasars and Implications for the MgII Broad-line Region

We present new Gemini/GMOS optical spectroscopy of 16 extreme variability quasars (EVQs) that dimmed by more than 1.5 mag in the $g$ band between the Sloan Digital Sky Survey (SDSS) and the Dark Energy Survey (DES) epochs (separated by a few years in the quasar rest frame). The quasar sample covers a redshift range of $0.5 < z < 2.1$. Nearly half of these EVQs brightened significantly (by more than 0.5 mag in the $g$ band) in a few years after reaching their previous faintest state, and some EVQs showed rapid (non-blazar) variations of greater than 1-2 mag on timescales of only months. Leveraging on the large dynamic range in continuum variability between the earlier SDSS and the new GMOS spectra, we explore the associated variations in the broad Mg II,$\lambda2798$ line, whose variability properties have not been well studied before. The broad Mg II flux varies in the same direction as the continuum flux, albeit with a smaller amplitude, which indicates at least some portion of Mg II is reverberating to continuum changes. However, the width (FWHM) of Mg II does not vary accordingly as continuum changes for most objects in the sample, in contrast to the case of the broad Balmer lines. Using the width of broad Mg II to estimate the black hole mass therefore introduces a luminosity-dependent bias.

preprint2019arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Accretion and Broad Emission Line Physics from a Hypervariable Quasar

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.

preprint2019arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Initial CIV Lag Results from Four Years of Data

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.

preprint2019arXiv

Understanding Broad Mg II Variability in Quasars with Photoionization: Implications for Reverberation Mapping and Changing-Look Quasars

The broad Mg II line in quasars has distinct variability properties compared with broad Balmer lines: it is less variable, and usually does not display a "breathing" mode, the increase in the average cloud distance when luminosity increases. We demonstrate that these variability properties of Mg II can be reasonably well explained by simple Locally Optimally Emitting Cloud (LOC) photoionization models, confirming earlier photoionization results. In the fiducial LOC model, the Mg II-emitting gas is on average more distant from the ionizing source than the H$α$/H$β$ gas, and responds with a lower amplitude to continuum variations. If the broad-line region (BLR) is truncated at a physical radius of $\sim 0.3$ pc (for a $10^{8.5}M_{\odot}$ BH accreting at Eddington ratio of 0.1), most of the Mg II flux will always be emitted near this outer boundary and hence will not display breathing. These results indicate that reverberation mapping results on broad Mg II, while generally more difficult to obtain due to the lower line responsivity, can still be used to infer the Mg II BLR size and hence black hole mass. But it is possible that Mg II does not have a well defined intrinsic BLR size-luminosity relation for individual quasars, even though a global one for the general population may still exist. The dramatic changes in broad H$α$/H$β$ emission in the observationally-rare changing-look quasars are fully consistent with photoionization responses to extreme continuum variability, and the LOC model provides natural explanations for the persistence of broad Mg II in changing-look quasars defined on H$α$/H$β$, and the rare population of broad Mg II emitters in the spectra of massive inactive galaxies.

preprint2019arXiv

Varstrometry for Off-nucleus and Dual sub-Kpc AGN (VODKA): Methodology and Initial Results with Gaia DR2

Gaia's milli-arcsec (mas) astrometric precision allows systematic identification of optically-selected sub-kpc dual active galactic nuclei (AGN), off-nucleus AGN, and small-scale lensed quasars by `varstrometry' -- where variability-induced astrometric jitter, i.e., temporal displacements of photocenter in unresolved sources, can be reasonably well detected or constrained. This approach extends systematic searches for small-scale ($\gtrsim$ mas) dual and off-nucleus AGN to poorly explored regime between $\sim 10$ pc and $\sim 1$ kpc, with Gaia's full sky coverage and depth to $G\sim 21$. We outline the general principles of this method and calculate the expected astrometric signals from the full time series of photocenter measurements and light curves. We demonstrate the feasibility of varstrometry by using Gaia DR2 data on a sample of variable pre-main sequence stars with known close companions. We find that extended host galaxies have a significant impact on the accuracy of astrometric and photometric variability in Gaia DR2, a situation to be improved in future Gaia releases. Using spectroscopically confirmed SDSS quasars, we present several examples of candidate sub-kpc off-nucleus or dual AGN selected from Gaia DR2. We discuss the merits and limitations of this method and follow-up strategy for promising candidates. We highlight Gaia's potential of systematically discovering and characterizing the sub-kpc off-nucleus and dual AGN population in the entire optical sky.

preprint2016arXiv

Chandra X-ray and Hubble Space Telescope Imaging of Optically Selected Kiloparsec-Scale Binary Active Galactic Nuclei II: Host Galaxy Morphology and AGN Activity

Binary active galactic nuclei (AGNs) provide clues to how gas-rich mergers trigger and fuel AGNs and how supermassive black hole (SMBH) pairs evolve in a gas-rich environment. While significant effort has been invested in their identification, the detailed properties of binary AGNs and their host galaxies are still poorly constrained. In a companion paper, we examined the nature of ionizing sources in the double nuclei of four kpc-scale binary AGNs with redshifts between 0.1~0.2. Here, we present their host galaxy morphology based on F336W (U-band) and F105W (Y-band) images taken by the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope. Our targets have double-peaked narrow emission lines and were confirmed to host binary AGNs with follow up observations. We find that kpc-scale binary AGNs occur in galaxy mergers with diverse morphological types. There are three major mergers with intermediate morphologies and a minor merger with a dominant disk component. We estimate the masses of the SMBHs from their host bulge stellar masses and obtain Eddington ratios for each AGN. Compared with a representative control sample drawn at the same redshift and stellar mass, the AGN luminosities and Eddington ratios of our binary AGNs are similar to those of single AGNs. The U-Y color maps indicate that clumpy star forming regions could significantly affect the X-ray detection of binary AGNs, e.g., the hardness ratio. Considering the weak X-ray emission in AGNs triggered in merger systems, we suggest that samples of X-ray selected AGNs may be biased against gas-rich mergers.

preprint2016arXiv

SDSS J0159+0105: A Radio-Quiet Quasar with a Centi-Parsec Supermassive Black Hole Binary Candidate

We report a candidate centi-parsec supermassive black hole binary (SMBHB) in the radio-quiet quasar SDSS J0159+0105 at z=0.217. With a modified lomb-scargle code GLSdeDRW and the auto-correlation analysis ACF, we detect two significant (at P>99%) periodic signals at ~741 day and ~1500 day from the 8.1-year Catalina V-band light curve of this quasar. The period ratio, which is close to 1:2, is typical of a black-hole binary system with a mass ratio of 0.05<q<0.8 according to recent numerical simulations. SDSS J0159+0105 has two SDSS spectroscopic observations separated by ~10 years. There is a significant change in the broad H-beta profile between the two epochs, which can be explained by a single broad-line region (BLR) around the binary system illuminated by the aforementioned mini-disks, or a stream of gas flowing from the circumbinary disk to one of the SMBHs. From the single BLR assumption and the orbital period t_orb ~1500 day, we estimate the total virial masses of M_SMBHB ~ 1.3x10^8 M_sun, the average distances of BLR of ~0.04pc (~50 light-day, with +/-0.3 dex uncertainty), and a SMBHB separation of d= (0.01pc)M_{8,tot}^{1/3} (T_rest/3.3yr)^{2/3} ~ 0.013 pc (15 light-day). Based on analytical work, the postulated circumbinary disk has an inner radius of 2d = 0.026 pc (30 light-day). SDSS J0159+0105 also displays unusual spectral energy distribution. The unique properties of SDSS J0159+0105 are consistent with it being a centi-parsec SMBHB. The GLSdeDRW code link: http://butler.lab.asu.edu/qso_period/

preprint2016arXiv

The Detailed Science Case for the Maunakea Spectroscopic Explorer: the Composition and Dynamics of the Faint Universe

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

preprint2016arXiv

The Final SDSS High-Redshift Quasar Sample of 52 Quasars at z>5.7

We present the discovery of nine quasars at $z\sim6$ identified in the Sloan Digital Sky Survey (SDSS) imaging data. This completes our survey of $z\sim6$ quasars in the SDSS footprint. Our final sample consists of 52 quasars at $5.7<z\le6.4$, including 29 quasars with $z_{\rm AB}\le20$ mag selected from 11,240 deg$^2$ of the SDSS single-epoch imaging survey (the main survey), 10 quasars with $20\le z_{\rm AB}\le20.5$ selected from 4223 deg$^2$ of the SDSS overlap regions (regions with two or more imaging scans), and 13 quasars down to $z_{\rm AB}\approx22$ mag from the 277 deg$^2$ in Stripe 82. They span a wide luminosity range of $-29.0\le M_{1450}\le-24.5$. This well-defined sample is used to derive the quasar luminosity function (QLF) at $z\sim6$. After combining our SDSS sample with two faint ($M_{1450}\ge-23$ mag) quasars from the literature, we obtain the parameters for a double power-law fit to the QLF. The bright-end slope $β$ of the QLF is well constrained to be $β=-2.8\pm0.2$. Due to the small number of low-luminosity quasars, the faint-end slope $α$ and the characteristic magnitude $M_{1450}^{\ast}$ are less well constrained, with $α=-1.90_{-0.44}^{+0.58}$ and $M^{\ast}=-25.2_{-3.8}^{+1.2}$ mag. The spatial density of luminous quasars, parametrized as $ρ(M_{1450}<-26,z)=ρ(z=6)\,10^{k(z-6)}$, drops rapidly from $z\sim5$ to 6, with $k=-0.72\pm0.11$. Based on our fitted QLF and assuming an IGM clumping factor of $C=3$, we find that the observed quasar population cannot provide enough photons to ionize the $z\sim6$ IGM at $\sim90$\% confidence. Quasars may still provide a significant fraction of the required photons, although much larger samples of faint quasars are needed for more stringent constraints on the quasar contribution to reionization.

preprint2016arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: An Investigation of Biases in CIV Emission-Line Properties

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.

preprint2016arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Velocity Shifts of Quasar Emission Lines

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.

preprint2016arXiv

The Time-Domain Spectroscopic Survey: Understanding the Optically Variable Sky with SEQUELS in SDSS-III

The Time-Domain Spectroscopic Survey (TDSS) is an SDSS-IV eBOSS subproject primarily aimed at obtaining identification spectra of ~220,000 optically-variable objects systematically selected from SDSS/Pan-STARRS1 multi-epoch imaging. We present a preview of the science enabled by TDSS, based on TDSS spectra taken over ~320 deg^2 of sky as part of the SEQUELS survey in SDSS-III, which is in part a pilot survey for eBOSS in SDSS-IV. Using the 15,746 TDSS-selected single-epoch spectra of photometrically variable objects in SEQUELS, we determine the demographics of our variability-selected sample, and investigate the unique spectral characteristics inherent in samples selected by variability. We show that variability-based selection of quasars complements color-based selection by selecting additional redder quasars, and mitigates redshift biases to produce a smooth quasar redshift distribution over a wide range of redshifts. The resulting quasar sample contains systematically higher fractions of blazars and broad absorption line quasars than from color-selected samples. Similarly, we show that M-dwarfs in the TDSS-selected stellar sample have systematically higher chromospheric active fractions than the underlying M-dwarf population, based on their H-alpha emission. TDSS also contains a large number of RR Lyrae and eclipsing binary stars with main-sequence colors, including a few composite-spectrum binaries. Finally, our visual inspection of TDSS spectra uncovers a significant number of peculiar spectra, and we highlight a few cases of these interesting objects. With a factor of ~15 more spectra, the main TDSS survey in SDSS-IV will leverage the lessons learned from these early results for a variety of time-domain science applications.

preprint2016arXiv

The UV-bright Quasar Survey (UVQS): DR1

We present the first data release (DR1) from our UV-bright Quasar Survey (UVQS) for new $z \sim 1$ active galactic nuclei (AGN) across the sky. Using simple GALEX UV and WISE near-IR color selection criteria, we generated a list of 1450 primary candidates with $FUV < 18.5$ mag. We obtained discovery spectra, primarily on 3m-class telescopes, for 1040 of these candidates and confirmed 86% as AGN with redshifts generally at $z>0.5$. Including a small set of observed secondary candidates, we report the discovery of 217 AGN with $FUV < 18$ mag that had no previously reported spectroscopic redshift. These are excellent potential targets for UV spectroscopy before the end of the {\it Hubble Space Telescope} mission. The main data products are publicly released through the Mikulski Archive for Space Telescopes.

preprint2016arXiv

Towards an Understanding of Changing-Look Quasars: An Archival Spectroscopic Search in SDSS

The uncertain origin of the recently-discovered `changing-looking' quasar phenomenon -- in which a luminous quasar dims significantly to a quiescent state in repeat spectroscopy over ~10 year timescales -- may present unexpected challenges to our understanding of quasar accretion. To better understand this phenomenon, we take a first step to building a sample of changing-look quasars with a systematic but simple archival search for these objects in the Sloan Digital Sky Survey Data Release 12. By leveraging the >10 year baselines for objects with repeat spectroscopy, we uncover two new changing-look quasars, and a third discovered previously. Decomposition of the multi-epoch spectra and analysis of the broad emission lines suggest that the quasar accretion disk emission dims due to rapidly decreasing accretion rates (by factors of >2.5), while disfavoring changes in intrinsic dust extinction for the two objects where these analyses are possible. Broad emission line energetics also support intrinsic dimming of quasar emission as the origin for this phenomenon rather than transient tidal disruption events or supernovae. Although our search criteria included quasars at all redshifts and transitions from either quasar-like to galaxy-like states or the reverse, all of the clear cases of changing-look quasars discovered were at relatively low-redshift (z ~ 0.2 - 0.3) and only exhibit quasar-like to galaxy-like transitions.

preprint2016arXiv

X-ray spectral properties of the AGN sample in the northern XMM-XXL field

In this paper we describe and publicly release a catalogue consisting of 8445 point-like X-ray sources detected in the XMM-XXL north survey. For the 2512 AGN which have reliable spectroscopy, we present the X-ray spectral fitting which has been computed with a Bayesian approach. We have also applied an X-ray spectral stacking method to different sub-samples. We confirm the well-known Iwasawa-Taniguchi effect in our luminosity-redshift sub-samples, and argue that such an effect is due to a decrease in the covering factor of a distant obscuring torus with increasing X-ray luminosity. By comparing the distribution of the reflection fraction, the ratio of the normalization of the reflected component to the direct radiation, we find that the low-luminosity, low-redshift sub-sample had systematically higher reflection fraction values than the high-redshift, high-luminosity one. On the other hand, no significant difference is found between samples having similar luminosity but different redshift, suggesting that the structure of the torus does not evolve strongly with redshift. Contrary to previous works, we do not find evidence for an increasing photon index at high Eddington ratio. This may be an indication that the structure of the accretion disc changes as the Eddington ratio approaches unity. Comparing our X-ray spectral analysis results with the optical spectral classification, we find that $\sim20$ per cent of optical type-1 AGN show an X-ray absorbing column density higher than $10^{21.5}\,\rm{cm^{-2}}$, and about 50 per cent of type-2 AGN have an X-ray absorbing column density less than $10^{21.5}\,\rm{cm^{-2}}$. We suggest that the excess X-ray absorption shown in the high-luminosity optical type-1 AGN can be due to small-scale dust-free gas within (or close to) the broad line region, while in the low-luminosity ones it can be due to a clumpy torus with a large covering factor.

preprint2015arXiv

Clustering of intermediate redshift quasars using the final SDSS III-BOSS sample

We measure the two-point clustering of spectroscopically confirmed quasars from the final sample of the Baryon Oscillation Spectroscopic Survey (BOSS) on comoving scales of 4 < s < 22 Mpc/h. The sample covers 6950 deg^2 (~ 19 (Gpc/h)^3) and, over the redshift range 2.2 < z < 2.8, contains 55,826 homogeneously selected quasars, which is twice as many as in any similar work. We deduce b_Q = 3.54 +/- 0.10 ; the most precise measurement of quasar bias to date at these redshifts. This corresponds to a host halo mass of ~ 2 x 10^12 ~ M_sun/h with an implied quasar duty cycle of ~1 percent. The real-space projected correlation function is well-fit by a power law of index -2 and correlation length r0 = (8.12 +/- 0.22), Mpc/h over scales of 4 < rp < 25 ~ Mpc/h. To better study the evolution of quasar clustering at moderate redshift, we extend the redshift range of our study to z ~ 3.4 and measure the bias and correlation length of three subsamples over 2.2 < z < 3.4. We find no significant evolution of r0 or bias over this range, implying that the host halo mass of quasars decreases somewhat with increasing redshift. We find quasar clustering remains similar over a decade in luminosity, contradicting a scenario in which quasar luminosity is monotonically related to halo mass at z ~ 2.5. Our results are broadly consistent with previous BOSS measurements, but they yield more precise constraints based upon a larger and more uniform data set.

preprint2015arXiv

Dissecting the quasar main sequence: insight from host galaxy properties

The diverse properties of broad-line quasars appear to follow a well-defined main sequence along which the optical FeII strength increases. It has been suggested that this sequence is mainly driven by the Eddington ratio (L/L_Edd) of the black hole (BH) accretion. Shen & Ho demonstrated with quasar clustering analysis that the average BH mass decreases with increasing FeII strength when quasar luminosity is fixed, consistent with this suggestion. Here we perform an independent test by measuring the stellar velocity dispersion sigma* (hence the BH mass via the M-sigma* relation) from decomposed host spectra in low-redshift Sloan Digital Sky Survey quasars. We found that at fixed quasar luminosity, sigma* systematically decreases with increasing FeII strength, confirming that Eddington ratio increases with FeII strength. We also found that at fixed luminosity and FeII strength, there is little dependence of sigma* on the broad Hbeta FWHM. These new results reinforce the framework put forward by Shen & Ho that Eddington ratio and orientation govern most of the diversity seen in broad-line quasar properties.

preprint2015arXiv

Now You See It, Now You Don't: The Disappearing Central Engine of the Quasar J1011+5442

We report the discovery of a new "changing-look" quasar, SDSS J101152.98+544206.4, through repeat spectroscopy from the Time Domain Spectroscopic Survey. This is an addition to a small but growing set of quasars whose blue continua and broad optical emission lines have been observed to decline by a large factor on a time scale of approximately a decade. The 5100 Angstrom monochromatic continuum luminosity of this quasar drops by a factor of > 9.8 in a rest-frame time interval of < 9.7 years, while the broad H-alpha luminosity drops by a factor of 55 in the same amount of time. The width of the broad H-alpha line increases in the dim state such that the black hole mass derived from the appropriate single-epoch scaling relation agrees between the two epochs within a factor of 3. The fluxes of the narrow emission lines do not appear to change between epochs. The light curve obtained by the Catalina Sky Survey suggests that the transition occurs within a rest-frame time interval of approximately 500 days. We examine three possible mechanisms for this transition suggested in the recent literature. An abrupt change in the reddening towards the central engine is disfavored by the substantial difference between the timescale to obscure the central engine and the observed timescale of the transition. A decaying tidal disruption flare is consistent with the decay rate of the light curve but not with the prolonged bright state preceding the decay, nor can this scenario provide the power required by the luminosities of the emission lines. An abrupt drop in the accretion rate onto the supermassive black hole appears to be the most plausible explanation for the rapid dimming.

preprint2015arXiv

Rest-frame Optical Properties of Luminous 1.5<z<3.5 Quasars: the Hbeta-[OIII] Region

We study the rest-frame optical properties of 74 luminous (L_bol=10^46.2-48.2 erg/s), 1.5<z<3.5 broad-line quasars with near-IR (JHK) slit spectroscopy. Systemic redshifts based on the peak of the [OIII]5007 line reveal that redshift estimates from the rest-frame UV broad emission lines (mostly MgII) are intrinsically uncertain by ~ 200 km/s (measurement errors accounted for). The overall full-width-at-half-maximum of the narrow [OIII] line is ~ 1000 km/s on average. A significant fraction of the total [OIII] flux (~ 40%) is in a blueshifted wing component with a median velocity offset of ~ 700 km/s, indicative of ionized outflows within a few kpc from the nucleus; we do not find evidence of significant [OIII] flux beyond ~ 10 kpc in our slit spectroscopy. The [OIII] line is noticeably more asymmetric and weaker than that in typical less luminous low-z quasars. However, when matched in quasar continuum luminosity, low-z quasars have similar [OIII] profiles and strengths as these high-z systems. Therefore the exceptionally large width and blueshifted wing, and the relatively weak strength of [OIII] in high-z luminous quasars are mostly a luminosity effect rather than redshift evolution. The Hbeta-[OIII] region of these high-z quasars displays a similar spectral diversity and Eigenvector 1 correlations with anti-correlated [OIII] and optical FeII strengths, as seen in low-z quasars; but the average broad Hbeta width is larger by 25% than typical low-z quasars, indicating more massive black holes in these high-z systems. These results highlight the importance of understanding [OIII] emission in the general context of quasar parameter space in order to understand quasar feedback in the form of [OIII] outflows. The calibrated one-dimensional near-IR spectra are made publicly available, along with a composite spectrum.

preprint2015arXiv

Reverberation Mapping with Intermediate-Band Photometry: Detection of Broad-Line H$α$ Time Lags for Quasars at $0.2<z<0.4$

We present a reverberation mapping (RM) experiment that combines broad- and intermediate-band photometry; it is the first such attempt targeting a sample of 13 quasars at $0.2<z<0.9$. The quasars were selected to have strong H$α$ or H$β$ emission lines that are located in one of three intermediate bands (with FWHM around 200 Å) centered at 8045, 8505, and 9171 Å. The imaging observations were carried out in the intermediate bands and the broad $i$ and $z$ bands using the prime-focus imager 90Prime on the 2.3m Bok telescope. Because of the large ($\sim$1 deg$^2$) field-of-view (FoV) of 90Prime, we were able to include the 13 quasars within only five telescope pointings or fields. The five fields were repeatedly observed over 20--30 epochs that were unevenly distributed over a duration of 5--6 months. The combination of the broad- and intermediate-band photometry allows us to derive accurate light curves for both optical continuum (from the accretion disk) and line (from the broad-line region, or BLR) emission. We detect H$α$ time lags between the continuum and line emission in 6 quasars. These quasars are at a relatively low redshift range $0.2<z<0.4$. The measured lags are consistent with the current BLR size-luminosity relation for H$β$ at $z<0.3$. While this experiment appears successful in detecting lags of the bright H$α$ line, further investigation is required to see if it can also be applied to the fainter H$β$ line for quasars at higher redshifts. Finally we demonstrate that by using a small telescope with a large FoV, intermediate-band photometric RM can be efficiently executed for a large sample of quasars at $z>0.2$.

preprint2015arXiv

The Eleventh and Twelfth Data Releases of the Sloan Digital Sky Survey: Final Data from SDSS-III

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

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Ensemble Spectroscopic Variability of Quasar Broad Emission Lines

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 ($Δ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 $δ$-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.

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: First Broad-line Hbeta and MgII Lags at z>~0.3 from six-Month Spectroscopy

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.

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: No Evidence for Evolution in the M-sigma Relation to z~1

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.

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Post-Starburst Signatures in Quasar Host Galaxies at z < 1

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.

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Rapid CIV Broad Absorption Line Variability

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.

preprint2015arXiv

The Time Domain Spectroscopic Survey: Variable Object Selection and Anticipated Results

We present the selection algorithm and anticipated results for the Time Domain Spectroscopic Survey (TDSS). TDSS is an SDSS-IV eBOSS subproject that will provide initial identification spectra of approximately 220,000 luminosity-variable objects (variable stars and AGN) across 7,500 square degrees selected from a combination of SDSS and multi-epoch Pan-STARRS1 photometry. TDSS will be the largest spectroscopic survey to explicitly target variable objects, avoiding pre-selection on the basis of colors or detailed modeling of specific variability characteristics. Kernel Density Estimate (KDE) analysis of our target population performed on SDSS Stripe 82 data suggests our target sample will be 95% pure (meaning 95% of objects we select have genuine luminosity variability of a few magnitudes or more). Our final spectroscopic sample will contain roughly 135,000 quasars and 85,000 stellar variables, approximately 4,000 of which will be RR Lyrae stars which may be used as outer Milky Way probes. The variability-selected quasar population has a smoother redshift distribution than a color-selected sample, and variability measurements similar to those we develop here may be used to make more uniform quasar samples in large surveys. The stellar variable targets are distributed fairly uniformly across color space, indicating that TDSS will obtain spectra for a wide variety of stellar variables including pulsating variables, stars with significant chromospheric activity, cataclysmic variables and eclipsing binaries. TDSS will serve as a pathfinder mission to identify and characterize the multitude of variable objects that will be detected photometrically in even larger variability surveys such as LSST.

preprint2015arXiv

X-ray Insights into the Nature of PHL 1811 Analogs and Weak Emission-Line Quasars: Unification with a Geometrically Thick Accretion Disk?

We present an X-ray and multiwavelength study of 33 weak emission-line quasars (WLQs) and 18 quasars that are analogs of the extreme WLQ, PHL 1811, at z ~ 0.5-2.9. New Chandra 1.5-9.5 ks exploratory observations were obtained for 32 objects while the others have archival X-ray observations. Significant fractions of these luminous type 1 quasars are distinctly X-ray weak compared to typical quasars, including 16 (48%) of the WLQs and 17 (94%) of the PHL 1811 analogs with average X-ray weakness factors of 17 and 39, respectively. We measure a relatively hard ($Γ=1.16_{-0.32}^{+0.37}$) effective power-law photon index for a stack of the X-ray weak subsample, suggesting X-ray absorption, and spectral analysis of one PHL 1811 analog, J1521+5202, also indicates significant intrinsic X-ray absorption. We compare composite SDSS spectra for the X-ray weak and X-ray normal populations and find several optical-UV tracers of X-ray weakness; e.g., Fe II rest-frame equivalent width and relative color. We describe how orientation effects under our previously proposed "shielding-gas" scenario can likely unify the X-ray weak and X-ray normal populations. We suggest that the shielding gas may naturally be understood as a geometrically thick inner accretion disk that shields the broad line region from the ionizing continuum. If WLQs and PHL 1811 analogs have very high Eddington ratios, the inner disk could be significantly puffed up (e.g., a slim disk). Shielding of the broad emission-line region by a geometrically thick disk may have a significant role in setting the broad distributions of C IV rest-frame equivalent width and blueshift for quasars more generally.

preprint2014arXiv

A hunt for dual radio active galactic nuclei in the VLASS

This whitepaper describes how the VLASS could be designed in a manner to allow the identification of candidate dual active galactic nuclei (AGN) at separations <7 kpc. Dual AGN represent a clear marker of two supermassive black holes within an ongoing merger. A dual AGN survey will provide a wealth of studies in structure growth and gravitational-wave science. Radio wavelengths are ideal for identifying close pairs, as disturbed stellar and gaseous material can obscure their presence in optical and shorter wavelengths. With sufficiently high resolution and sensitivity, a large-scale radio imaging survey like the VLASS will uncover many of these systems and provide the means to broadly study the radio properties of candidate dual systems revealed at other wavelengths. We determine that the ideal survey for our purposes will be at as high a resolution as possible, with significantly more science return in A array at L-band or higher, or B array at C-band or higher. We describe a range of potential survey parameters within this document. Based on the analysis outlined in this whitepaper, our ideal survey would create a catalogue of $\gtrsim$100 dual AGN in either: 1) a medium-sensitivity (~1 mJy detection threshold), wide-field (few thousand square degree) survey, or 2) a high-sensitivity (~10 uJy threshold) survey of several hundred square degrees.

preprint2014arXiv

A real-time fast radio burst: polarization detection and multiwavelength follow-up

Fast radio bursts (FRBs) are one of the most tantalizing mysteries of the radio sky; their progenitors and origins remain unknown and until now no rapid multiwavelength follow-up of an FRB has been possible. New instrumentation has decreased the time between observation and discovery from years to seconds, and enables polarimetry to be performed on FRBs for the first time. We have discovered an FRB (FRB 140514) in real-time on 14 May, 2014 at 17:14:11.06 UTC at the Parkes radio telescope and triggered follow-up at other wavelengths within hours of the event. FRB 140514 was found with a dispersion measure (DM) of 562.7(6) cm$^{-3}$ pc, giving an upper limit on source redshift of $z \lesssim 0.5$. FRB 140514 was found to be 21$\pm$7% (3-$σ$) circularly polarized on the leading edge with a 1-$σ$ upper limit on linear polarization $<10%$. We conclude that this polarization is intrinsic to the FRB. If there was any intrinsic linear polarization, as might be expected from coherent emission, then it may have been depolarized by Faraday rotation caused by passing through strong magnetic fields and/or high density environments. FRB 140514 was discovered during a campaign to re-observe known FRB fields, and lies close to a previous discovery, FRB 110220; based on the difference in DMs of these bursts and time-on-sky arguments, we attribute the proximity to sampling bias and conclude that they are distinct objects. Follow-up conducted by 12 telescopes observing from X-ray to radio wavelengths was unable to identify a variable multiwavelength counterpart, allowing us to rule out models in which FRBs originate from nearby ($z < 0.3$) supernovae and long duration gamma-ray bursts.

preprint2014arXiv

Constraining Sub-Parsec Binary Supermassive Black Holes in Quasars with Multi-Epoch Spectroscopy. II. The Population with Kinematically Offset Broad Balmer Emission Lines

A small fraction of quasars have long been known to show bulk velocity offsets in the broad Balmer lines with respect to the systemic redshift of the host galaxy. Models to explain these offsets usually invoke broad-line region gas kinematics/asymmetry around single black holes (BHs), orbital motion of massive (~sub-pc) binary black holes (BBHs), or recoil BHs, but single-epoch spectra are unable to distinguish between these scenarios. The line-of-sight (LOS) radial velocity (RV) shifts from long-term spectroscopic monitoring can be used to test the BBH hypothesis. We have selected a sample of 399 quasars with offset broad H-beta lines from the SDSS DR7 quasar catalog, and have conducted second-epoch optical spectroscopy for 50 of them. Combined with the existing SDSS spectra, the new observations enable us to constrain the LOS RV shifts of broad H-beta lines with a rest-frame baseline of a few years to nearly a decade. Using cross-correlation analysis, we detect significant (99% confidence) radial accelerations in the broad H-beta lines in 24 of the 50 objects. We suggest that 9 of the 24 detections are sub-pc BBH candidates, which show consistent velocity shifts independently measured from a second broad line (either H-alpha or Mg II) without significant changes in the broad-line profiles. Combining the results on the general quasar population studied in Paper I, we find a tentative anti-correlation between the velocity offset in the first-epoch spectrum and the average acceleration between two epochs, which could be explained by orbital phase modulation when the time separation between two epochs is a non-negligible fraction of the orbital period of the motion causing the line displacement. We discuss the implications of our results for the identification of sub-pc BBH candidates in offset-line quasars and for the constraints on their frequency and orbital parameters. [Abridged]

preprint2014arXiv

Molecular-Scale Hydrophilicity Induced by Solute: Molecular-thick Charged Pancakes of Aqueous Salt Solution on Hydrophobic Carbon-based Surfaces

We directly observed molecular-thick aqueous salt-solution pancakes on a hydrophobic graphite surface under ambient conditions employing atomic force microscopy. This observation indicates the unexpected molecular-scale hydrophilicity of the salt solution on graphite surfaces, which is different from the macroscopic wetting property of a droplet standing on the graphite surface. Interestingly, the pancakes spontaneously displayed strong positively charged behavior. Theoretical studies showed that the formation of such positively charged pancakes is attributed to cation-π interactions between Na+ ions in the aqueous solution and aromatic rings on the graphite surface, promoting the adsorption of water molecules together with cations onto the graphite surface; i.e., Na+ ions as a medium adsorbed to the graphite surface through cation-π interactions on one side while at the same time bonding to water molecules through hydration interaction on the other side at a molecular scale. These findings suggest that actual interactions regarding carbon-based graphitic surfaces including those of graphene, carbon nanotubes, and biochar may be significantly different from existing theory and they provide new insight into the control of surface wettability, interactions and related physical, chemical and biological processes.

preprint2014arXiv

Relativistic redshifts in quasar broad lines

The broad emission lines commonly seen in quasar spectra have velocity widths of a few per cent of the speed of light, so special- and general-relativistic effects have a significant influence on the line profile. We have determined the redshift of the broad H-beta line in the quasar rest frame (determined from the core component of the [OIII] line) for over 20,000 quasars from the Sloan Digital Sky Survey DR7 quasar catalog. The mean redshift as a function of line width is approximately consistent with the relativistic redshift that is expected if the line originates in a randomly oriented Keplerian disk that is obscured when the inclination of the disk to the line of sight exceeds ~30-45 degrees, consistent with simple AGN unification schemes. This result also implies that the net line-of-sight inflow/outflow velocities in the broad-line region are much less than the Keplerian velocity when averaged over a large sample of quasars with a given line width.

preprint2014arXiv

Robust OS-ELM with a novel selective ensemble based on particle swarm optimization

In this paper, a robust online sequential extreme learning machine (ROS-ELM) is proposed. It is based on the original OS-ELM with an adaptive selective ensemble framework. Two novel insights are proposed in this paper. First, a novel selective ensemble algorithm referred to as particle swarm optimization selective ensemble (PSOSEN) is proposed. Noting that PSOSEN is a general selective ensemble method which is applicable to any learning algorithms, including batch learning and online learning. Second, an adaptive selective ensemble framework for online learning is designed to balance the robustness and complexity of the algorithm. Experiments for both regression and classification problems with UCI data sets are carried out. Comparisons between OS-ELM, simple ensemble OS-ELM (EOS-ELM) and the proposed ROS-ELM empirically show that ROS-ELM significantly improves the robustness and stability.

preprint2014arXiv

Single-Epoch Black Hole Mass Estimators For Broad-Line Active Galactic Nuclei: Recalibrating Hbeta with A New Approach

Based on an updated Hbeta reverberation mapping (RM) sample of 44 nearby active galactic nuclei (AGN), we propose a novel approach for black hole (BH) mass estimation using two filtered luminosities computed from single-epoch (SE) AGN spectra around the Hbeta region. We found that the two optimal-filter luminosities extract virial information (size and virial velocity of the broad line region, BLR) from the spectra, justifying their usage in this empirical BH mass estimator. The major advantages of this new recipe over traditional SE BH mass estimators utilizing continuum luminosity and broad line width are: 1) it has a smaller intrinsic scatter of 0.28 dex calibrated against RM masses; 2) it is extremely simple to use in practice, without any need to decompose the spectrum; 3) it produces unambiguous and highly repeatable results even with low signal-to-noise spectra. The combination of the two luminosities can also cancel out, to some extent, systematic luminosity errors potentially introduced by uncertainties in distance or flux calibration. In addition, we recalibrated the traditional SE mass estimators using broad Hbeta full width at half maximum (FWHM) and monochromatic continuum luminosity at 5100A (L_5100). We found that using the best-fit slopes on FWHM and L_5100 (derived from fitting the BLR radius-luminosity relation and the correlation between rms line dispersion and SE FWHM, respectively) rather than simple assumptions (e.g., 0.5 for L_5100 and 2 for FWHM) leads to more precise SE mass estimates, improving the intrinsic scatter from 0.41 dex to 0.36 dex with respect to the RM masses. We compared different estimators and discussed their applications to the Sloan Digital Sky Survey (SDSS) quasar sample. Due to the limitations of the current RM sample, application of any SE recipe calibrated against RM masses to distant quasars should be treated with caution.

preprint2014arXiv

The Dependence of C IV Broad Absorption Line Properties on Accompanying Si IV and Al III Absorption: Relating Quasar-Wind Ionization Levels, Kinematics, and Column Densities

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.

preprint2014arXiv

The diversity of quasars unified by accretion and orientation

Quasars are rapidly accreting supermassive black holes at the center of massive galaxies. They display a broad range of properties across all wavelengths, reflecting the diversity in the physical conditions of the regions close to the central engine. These properties, however, are not random, but form well-defined trends. The dominant trend is known as Eigenvector 1, where many properties correlate with the strength of optical iron and [OIII] emission. The main physical driver of Eigenvector 1 has long been suspected to be the quasar luminosity normalized by the mass of the hole (the Eddington ratio), an important quantity of the black hole accretion process. But a definitive proof has been missing. Here we report an analysis of archival data that reveals that Eddington ratio indeed drives Eigenvector 1. We also find that orientation plays a significant role in determining the observed kinematics of the gas, implying a flattened, disklike geometry for the fast-moving clouds close to the hole. Our results show that most of the diversity of quasar phenomenology can be unified with two simple quantities, Eddington ratio and orientation.

preprint2014arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Technical Overview

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.

preprint2014arXiv

The Tenth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Apache Point Observatory Galactic Evolution Experiment

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

preprint2013arXiv

Broad Absorption Line Quasars with Redshifted Troughs: High-Velocity Infall or Rotationally Dominated Outflows?

We report the discovery in the Sloan Digital Sky Survey and the SDSS-III Baryon Oscillation Spectroscopic Survey of seventeen broad absorption line (BAL) quasars with high-ionization troughs that include absorption redshifted relative to the quasar rest frame. The redshifted troughs extend to velocities up to v=12,000 km/s and the trough widths exceed 3000 km/s in all but one case. Approximately 1 in 1000 BAL quasars with blueshifted C IV absorption also has redshifted C IV absorption; objects with C IV absorption present only at redshifted velocities are roughly four times rarer. In more than half of our objects, redshifted absorption is seen in C II or Al III as well as C IV, making low-ionization absorption at least ten times more common among BAL quasars with redshifted troughs than among standard BAL quasars. However, the C IV absorption equivalent widths in our objects are on average smaller than those of standard BAL quasars with low-ionization absorption. We consider several possible ways of generating redshifted absorption. The two most likely possibilities may be at work simultaneously, in the same objects or in different ones. Rotationally dominated outflows seen against a quasar's extended continuum source can produce redshifted and blueshifted absorption, but variability consistent with this scenario is seen in only one of the four objects with multiple spectra. The infall of relatively dense and low-ionization gas to radii as small as 400 Schwarzschild radii can in principle explain the observed range of trough profiles, but current models do not easily explain the origin and survival of such gas. Whatever the origin(s) of the absorbing gas in these objects, it must be located at small radii to explain its large redshifted velocities, and thus offers a novel probe of the inner regions of quasars.

preprint2013arXiv

Broad Absorption Line Variability on Multi-Year Timescales in a Large Quasar Sample

We present a detailed investigation of the variability of 428 C IV and 235 Si IV Broad Absorption Line (BAL) troughs identified in multi-epoch observations of 291 quasars by the Sloan Digital Sky Survey-I/II/III. These observations primarily sample rest-frame timescales of 1-3 yr over which significant rearrangement of the BAL wind is expected. We derive a number of observational results on, e.g., the frequency of BAL variability, the velocity range over which BAL variability occurs, the primary observed form of BAL-trough variability, the dependence of BAL variability upon timescale, the frequency of BAL strengthening vs. weakening, correlations between BAL variability and BAL-trough profiles, relations between C IV and Si IV BAL variability, coordinated multi-trough variability, and BAL variations as a function of quasar properties. We assess implications of these observational results for quasar winds. Our results support models where most BAL absorption is formed within an order-of-magnitude of the wind-launching radius, although a significant minority of BAL troughs may arise on larger scales. We estimate an average lifetime for a BAL trough along our line-of-sight of a few thousand years. BAL disappearance and emergence events appear to be extremes of general BAL variability, rather than being qualitatively distinct phenomena. We derive the parameters of a random-walk model for BAL EW variability, finding that this model can acceptably describe some key aspects of EW variability. The coordinated trough variability of BAL quasars with multiple troughs suggests that changes in "shielding gas" may play a significant role in driving general BAL variability.

preprint2013arXiv

Constraining sub-Parsec Binary Supermassive Black Holes in Quasars with Multi-Epoch Spectroscopy. I. The General Quasar Population

We perform a systematic search for sub-parsec binary supermassive black holes (BHs) in normal broad line quasars at z<0.8, using multi-epoch SDSS spectroscopy of the broad Hbeta line. Our working model is that: only one of the two BHs in the binary is active, and dynamically dominates its own broad line region (BLR); the inactive companion BH is orbiting at a distance of a few R_BLR, where R_BLR~0.01-0.1 pc is the BLR size. We search for the expected line-of-sight acceleration of the broad line velocity from binary orbital motion by cross-correlating SDSS spectra from two epochs separated by up to several years in the quasar restframe. Out of ~700 pairs of spectra for which we have good measurements of the velocity shift between two epochs (1-sigma error~40 km/s), we detect 28 systems with significant velocity shifts in broad Hbeta, among which seven are the best candidates for the hypothesized binaries. We use the distribution of the observed accelerations (mostly non-detections) to place constraints on the abundance of such binary systems among the general quasar population. Excess variance in the velocity shift is inferred for observations separated by longer than 0.4 yr (quasar restframe). Attributing all the excess to binary motion would imply that most of the quasars in this sample must be in binaries; that the inactive BH must be on average more massive than the active one; and that the binary separation is at most a few times the size of the BLR. However, if this excess variance is partly or largely due to long-term broad line variability, the requirement of a large population of close binaries is much weakened or even disfavored for massive companions. [Abridged]

preprint2013arXiv

The Distribution of AGN Covering Factors

We review our knowledge of the most basic properties of the AGN obscuring region - its location, scale, symmetry, and mean covering factor - and discuss new evidence on the distribution of covering factors in a sample of ~9000 quasars with WISE, UKIDSS, and SDSS photometry. The obscuring regions of AGN may be in some ways more complex than we thought - multi-scale, not symmetric, chaotic - and in some ways simpler - with no dependence on luminosity, and a covering factor distribution that may be determined by the simplest of considerations - e.g. random misalignments.

preprint2013arXiv

The Mass of Quasars

I review the current status of quasar black hole (BH) mass estimations. Spectroscopic methods have been developed to estimate BH mass in broad line quasars to an accuracy of ~0.5 dex. Despite their popularity, significant issues and confusion remain regarding these mass estimators. I provide an in-depth discussion on the merits and caveats of the single-epoch (SE) virial BH mass estimators, and a detailed derivation of the statistical biases of these SE mass estimates resulting from their errors. I show that error-induced sample biases on the order of a factor of several are likely present in the SE mass estimates for flux-limited, statistical quasar samples, and the distribution of SE masses in finite luminosity bins can be narrower than the nominal uncertainty of these mass estimates. I then discuss the latest applications of SE virial masses in quasar studies, including the early growth of supermassive black holes, quasar demography in the mass-luminosity plane, and the evolution of the BH-host scaling relations, with specific emphases on selection effects and sample biases in the SE masses. I conclude that there is a pressing need to understand and deal with the errors in these BH mass estimates, and to improve these BH weighing methods with substantially more and better reverberation mapping data.

preprint2013arXiv

The z=5 Quasar Luminosity Function from SDSS Stripe 82

We present a measurement of the Type I quasar luminosity function at z=5 using a large sample of spectroscopically confirmed quasars selected from optical imaging data. We measure the bright end (M_1450<-26) with Sloan Digital Sky Survey (SDSS) data covering ~6000 deg^2, then extend to lower luminosities (M_1450<-24) with newly discovered, faint z~5 quasars selected from 235 deg^2 of deep, coadded imaging in the SDSS Stripe 82 region (the celestial equator in the Southern Galactic Cap). The faint sample includes 14 quasars with spectra obtained as ancillary science targets in the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), and 59 quasars observed at the MMT and Magellan telescopes. We construct a well-defined sample of 4.7<z<5.1 quasars that is highly complete, with 73 spectroscopic identifications out of 92 candidates. Our color selection method is also highly efficient: of the 73 spectra obtained, 71 are high redshift quasars. These observations reach below the break in the luminosity function (M_1450* ~ -27). The bright end slope is steep (beta <~ -4), with a constraint of beta < -3.1 at 95% confidence. The break luminosity appears to evolve strongly at high redshift, providing an explanation for the flattening of the bright end slope reported previously. We find a factor of ~2 greater decrease in the number density of luminous quasars (M_1450<-26) from z=5 to z=6 than from z=4 to z=5, suggesting a more rapid decline in quasar activity at high redshift than found in previous surveys. Our model for the quasar luminosity function predicts that quasars generate ~30% of the ionizing photons required to keep the universe ionized at z=5.

preprint2012arXiv

Astrometric Reverberation Mapping

Spatially extended emission regions of Active Galactic Nuclei (AGN) respond to continuum variations, if such emission regions are powered by energy reprocessing of the continuum. The response from different parts of the reverberating region arrives at different times lagging behind the continuum variation. The lags can be used to map the geometry and kinematics of the emission region (i.e., reverberation mapping, RM). If the extended emission region is not spherically symmetric in configuration and velocity space, reverberation may produce astrometric offsets in the emission region photocenter as a function of time delay and velocity, detectable with future micro-arcsec to tens of micro-arcsec astrometry. Such astrometric responses provide independent constraints on the geometric and kinematic structure of the extended emission region, complementary to traditional reverberation mapping. In addition, astrometric RM is more sensitive to infer the inclination of a flattened geometry and the rotation angle of the extended emission region.

preprint2012arXiv

Broad Absorption Line Disappearance on Multi-Year Timescales in a Large Quasar Sample

We present 21 examples of C IV Broad Absorption Line (BAL) trough disappearance in 19 quasars selected from systematic multi-epoch observations of 582 bright BAL quasars (1.9 < z < 4.5) by the Sloan Digital Sky Survey-I/II (SDSS-I/II) and SDSS-III. The observations span 1.1-3.9 yr rest-frame timescales, longer than have been sampled in many previous BAL variability studies. On these timescales, ~2.3% of C IV BAL troughs disappear and ~3.3% of BAL quasars show a disappearing trough. These observed frequencies suggest that many C IV BAL absorbers spend on average at most a century along our line of sight to their quasar. Ten of the 19 BAL quasars showing C IV BAL disappearance have apparently transformed from BAL to non-BAL quasars; these are the first reported examples of such transformations. The BAL troughs that disappear tend to be those with small-to-moderate equivalent widths, relatively shallow depths, and high outflow velocities. Other non-disappearing C IV BALs in those nine objects having multiple troughs tend to weaken when one of them disappears, indicating a connection between the disappearing and non-disappearing troughs, even for velocity separations as large as 10000-15000 km/s. We discuss possible origins of this connection including disk-wind rotation and changes in shielding gas.

preprint2012arXiv

Chandra X-ray and Hubble Space Telescope Imaging of Optically Selected kiloparsec-Scale Binary Active Galactic Nuclei I. Nature of the Nuclear Ionizing Sources

Kiloparsec-scale binary active galactic nuclei (AGNs) signal active supermassive black hole (SMBH) pairs in merging galaxies. Despite their significance, unambiguously confirmed cases remain scarce and most have been discovered serendipitously. In a previous systematic search, we optically identified four kpc-scale binary AGNs from candidates selected with double-peaked narrow emission lines at redshifts of 0.1--0.2. Here we present Chandra and Hubble Space Telescope Wide Field Camera 3 (WFC3) imaging of these four systems. We critically examine and confirm the binary-AGN scenario for two of the four targets, by combining high angular resolution X-ray imaging spectroscopy with Chandra ACIS-S, better nuclear position constraints from WFC3 F105W imaging, and direct starburst estimates from WFC3 F336W imaging; for the other two targets, the existing data are still consistent with the binary-AGN scenario, but we cannot rule out the possibility of only one AGN ionizing gas in both merging galaxies. We find tentative evidence for a systematically smaller X-ray-to-[O III] luminosity ratio and/or higher Compton-thick fraction in optically selected kpc-scale binary AGNs than in single AGNs, possibly caused by a higher nuclear gas column due to mergers and/or a viewing angle bias related to the double-peak narrow line selection. While our result lends some further support to the general approach of optically identifying kpc-scale binary AGNs, it also highlights the challenge and ambiguity of X-ray confirmation.

preprint2012arXiv

Comparing Single-Epoch Virial Black Hole Mass Estimators for Luminous Quasars

Single-epoch virial black hole (BH) mass estimators utilizing broad emission lines have been routinely applied to high-redshift quasars to estimate their BH masses. Depending on the redshift, different line estimators (Halpha, Hbeta, MgII, CIV) are often used with optical/near-infrared spectroscopy. Here we use a homogeneous sample of 60 intermediate-redshift (z~1.5-2.2) SDSS quasars with optical and near-infrared spectra covering CIV through Halpha to investigate the consistency between different line estimators. We critically compare restframe UV line estimators (CIV, CIII], and MgII) with optical estimators (Hbeta and Halpha) in terms of correlations between line widths and between continuum/line luminosities, for the high-luminosity regime (L_5100>10^45.4 erg/s) probed by our sample. The continuum luminosities of L_1350 and L_3000, and the broad line luminosities are well correlated with L_5100. We found that the MgII FWHM correlates well with the FWHMs of the Balmer lines, and that the MgII line estimator can be calibrated to yield consistent virial mass estimates with those based on the Hbeta/Halpha estimators, thus extending earlier results on less luminous objects. The CIV FWHM is poorly correlated with the Balmer line FWHMs, and the scatter between the CIV and Hbeta FWHMs consists of an irreducible part (~0.12 dex), and a part that correlates with the blueshift of the CIV centroid relative to that of Hbeta. The CIII] FWHM is found to correlate with the CIV FWHM, and hence is also poorly correlated with the Hbeta FWHM. While the CIV and CIII] lines can be calibrated to yield consistent virial mass estimates as Hbeta on average, the scatter is substantially larger than MgII, and the usage of CIV/CIII] FWHM in the mass estimators does not improve the agreement with the Hbeta estimator. (Abridged)

preprint2012arXiv

Cross-Correlation of SDSS DR7 Quasars and DR10 BOSS Galaxies: The Weak Luminosity Dependence of Quasar Clustering at z~0.5

We present the measurement of the two-point cross-correlation function (CCF) of 8,198 Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) quasars and 349,608 DR10 CMASS galaxies from the Baryonic Oscillation Spectroscopic Survey (BOSS) at redshift <z>~0.5 (0.3<z<0.9). The cross-correlation function can be reasonably well fit by a power-law model xi_QG(r)=(r/r_0)^(-gamma) on projected scales of r_p=2-25 Mpc/h with r_0=6.61+-0.25 Mpc/h and gamma=1.69+-0.07. We estimate a quasar linear bias of b_Q=1.38+-0.10 at <z>=0.53 from the CCF measurements. This linear bias corresponds to a characteristic host halo mass of ~4x10^12 M_sun/h, compared to ~10^13 M_sun/h characteristic host halo mass for CMASS galaxies. We divide the quasar sample in luminosity and constrain the luminosity dependence of quasar bias to be db_Q/dlogL=0.20+-0.34 or 0.11+-0.32 (depending on different luminosity divisions) for quasar luminosities -23.5>M_i(z=2)>-25.5, implying a weak luminosity dependence of quasar clustering for the bright end of the quasar population at <z>~0.5. We compare our measurements with theoretical predictions, Halo Occupation Distribution (HOD) models and mock catalogs. These comparisons suggest quasars reside in a broad range of host halos, and the host halo mass distributions significantly overlap with each other for quasars at different luminosities, implying a poor correlation between halo mass and instantaneous quasar luminosity. We also find that the quasar HOD parameterization is largely degenerate such that different HODs can reproduce the CCF equally well, but with different outcomes such as the satellite fraction and host halo mass distribution. These results highlight the limitations and ambiguities in modeling the distribution of quasars with the standard HOD approach and the need for additional information in populating quasars in dark matter halos with HOD. [Abridged]

preprint2012arXiv

IR-derived covering factors for a large sample of quasars from WISE-UKIDSS-SDSS

We investigate the range of covering factors (determined from the ratio of IR to UV/optical luminosity) seen in luminous type 1 quasars using a combination of data from the WISE, UKIDSS and SDSS surveys. Accretion disk (UV/optical) and obscuring dust (IR) luminosities are measured via the use of a simple three component SED model. We use these estimates to investigate the distribution of covering factors and its relationship to both accretion luminosity and IR SED shape. The distribution of covering factors (f_C) is observed to be log-normal, with a bias-corrected mean of <log10 f_C>=-0.41 and standard deviation of 0.2. The fraction of IR luminosity emitted in the near-IR (1--5 micron) is found to be high (~40 per cent), and strongly dependant on covering factor.

preprint2012arXiv

On the link between associated MgII absorbers and star formation in quasar hosts

A few percent of quasars show strong associated MgII absorption, with velocities (v_off) lying within a few thousand km/s from the quasar systemic redshift. These associated absorption line systems (AALs) are usually interpreted as absorbers that are either intrinsic to the quasar and its host, or arising from external galaxies clustering around the quasar. Using composite spectra of ~1,800 MgII AAL quasars selected from SDSS DR7 at 0.4<~ z<~2, we show that quasars with AALs with v_off<1500 km/s have a prominent excess in [OII]3727 emission (detected at >7sigma) at rest to the quasar host, compared to unabsorbed quasars. We interpret this [OII] excess as due to enhanced star formation in the quasar host. Our results suggest that a significant fraction of AALs with v_off<1500 km/s are physically associated with the quasar and its host, and are most likely large-scale flows with typical velocities of a few hundred km/s. AAL quasars also have dust reddening lying between normal quasars and the so-called dust-reddened quasars. We suggest that the unique properties of AAL quasars can be explained if they are the transitional population from heavily dust-reddened quasar to normal quasars in the formation process of quasars and their hosts. This scenario predicts a larger fraction of young bulges, disturbed morphologies and interactions of AAL quasar hosts compared to normal quasars. The intrinsic link between associated absorbers and quasar hosts opens a new window to probe massive galaxy formation and galactic-scale feedback processes, and provides a crucial test of the evolutionary picture of quasars.

preprint2012arXiv

The Baryon Oscillation Spectroscopic Survey of SDSS-III

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

preprint2012arXiv

The clustering of intermediate redshift quasars as measured by the Baryon Oscillation Spectroscopic Survey

We measure the quasar two-point correlation function over the redshift range 2.2<z<2.8 using data from the Baryon Oscillation Spectroscopic Survey. We use a homogeneous subset of the data consisting of 27,129 quasars with spectroscopic redshifts---by far the largest such sample used for clustering measurements at these redshifts to date. The sample covers 3,600 square degrees, corresponding to a comoving volume of 9.7(Gpc/h)^3 assuming a fiducial LambdaCDM cosmology, and it has a median absolute i-band magnitude of -26, k-corrected to z=2. After accounting for redshift errors we find that the redshift space correlation function is fit well by a power-law of slope -2 and amplitude s_0=(9.7\pm 0.5)Mpc/h over the range 3<s<25Mpc/h. The projected correlation function, which integrates out the effects of peculiar velocities and redshift errors, is fit well by a power-law of slope -1 and r_0=(8.4\pm 0.6)Mpc/h over the range 4<R<16Mpc/h. There is no evidence for strong luminosity or redshift dependence to the clustering amplitude, in part because of the limited dynamic range in our sample. Our results are consistent with, but more precise than, previous measurements at similar redshifts. Our measurement of the quasar clustering amplitude implies a bias factor of b~3.5 for our quasar sample. We compare the data to models to constrain the manner in which quasars occupy dark matter halos at z~2.4 and infer that such quasars inhabit halos with a characteristic mass of <M>~10^{12}Msun/h with a duty cycle for the quasar activity of 1 per cent.

preprint2012arXiv

The Demographics of Broad Line Quasars in the Mass-Luminosity Plane II. Black Hole Mass and Eddington Ratio Functions

We employ a flexible Bayesian technique to estimate the black hole mass and Eddington ratio functions for Type 1 (i.e., broad line) quasars from a uniformly-selected data set of ~58,000 quasars from the SDSS DR7. We find that the SDSS becomes significantly incomplete at M_{BH} < 3 x 10^8 M_{Sun} or L / L_{Edd} < 0.07, and that the number densities of Type 1 quasars continue to increase down to these limits. Both the mass and Eddington ratio functions show evidence of downsizing, with the most massive and highest Eddington ratio black holes experiencing Type 1 quasar phases first, although the Eddington ratio number densities are flat at z < 2. We estimate the maximum Eddington ratio of Type 1 quasars in the observable Universe to be L / L_{Edd} ~ 3. Consistent with our results in Paper I, we do not find statistical evidence for a so-called "sub-Eddington boundary" in the mass-luminosity plane of broad line quasars, and demonstrate that such an apparent boundary in the observed distribution can be caused by selection effect and errors in virial BH mass estimates. Based on the typical Eddington ratio in a given mass bin, we estimate typical growth times for the black holes in Type 1 quasars and find that they are typically comparable to or longer than the age of the universe, implying an earlier phase of accelerated (i.e., with higher Eddington ratios) and possibly obscured growth. The large masses probed by our sample imply that most of our black holes reside in what are locally early type galaxies, and we interpret our results within the context of models of self-regulated black hole growth.

preprint2012arXiv

The Demographics of Broad-Line Quasars in the Mass-Luminosity Plane. I. Testing FWHM-Based Virial Black Hole Masses

We jointly constrain the luminosity function (LF) and black hole mass function (BHMF) of broad-line quasars with forward Bayesian modeling in the quasar mass-luminosity plane, based on a homogeneous sample of ~ 58,000 SDSS DR7 quasars at z ~ 0.3-5. We take into account the selection effect of the sample flux limit; more importantly, we deal with the statistical scatter between true BH masses and FWHM-based single-epoch virial mass estimates, as well as potential luminosity-dependent biases of these mass estimates. The LF is tightly constrained in the regime sampled by SDSS, and makes reasonable predictions when extrapolated to ~ 3 magnitudes fainter. Downsizing is seen in the model LF. On the other hand, we find it difficult to constrain the BHMF to within a factor of a few at z>~0.7 (with MgII and CIV-based virial BH masses). This is mainly driven by the unknown luminosity-dependent bias of these mass estimators and its degeneracy with other model parameters, and secondly driven by the fact that SDSS quasars only sample the tip of the active BH population at high redshift. Nevertheless, the most likely models favor a positive luminosity-dependent bias for MgII and possibly for CIV, such that at fixed true BH mass, objects with higher-than-average luminosities have over-estimated FWHM-based virial masses. There is tentative evidence that downsizing also manifests itself in the active BHMF, and the BH mass density in broad-line quasars contributes an insignificant amount to the total BH mass density at all times. Within our model uncertainties, we do not find a strong BH mass dependence of the mean Eddington ratio; but there is evidence that the mean Eddington ratio (at fixed BH mass) increases with redshift.

preprint2012arXiv

The Halo Occupation Distribution of SDSS Quasars

We present an estimate of the projected two-point correlation function (2PCF) of quasars in the Sloan Digital Sky Survey (SDSS) over the full range of one- and two-halo scales, 0.02-120 Mpc/h. This was achieved by combining data from SDSS DR7 on large scales and Hennawi et al. (2006; with appropriate statistical corrections) on small scales. Our combined clustering sample is the largest spectroscopic quasar clustering sample to date, containing ~48,000 quasars in the redshift range 0.4<z<2.5 with median redshift 1.4. We interpret these precise 2PCF measurements within the halo occupation distribution (HOD) framework and constrain the occupation functions of central and satellite quasars in dark matter halos. In order to explain the small-scale clustering, the HOD modeling requires that a small fraction of z~1.4 quasars, fsat=(7.4+/-1.4) 10^(-4), be satellites in dark matter halos. At z~1.4, the median masses of the host halos of central and satellite quasars are constrained to be Mcen=(4.1+0.3/-0.4) 10^12 Msun/h and Msat=(3.6+0.8/-1.0) 10^14 Msun/h, respectively. To investigate the redshift evolution of the quasar-halo relationship, we also perform HOD modeling of the projected 2PCF measured by Shen et al. (2007) for SDSS quasars with median redshift 3.2. We find tentative evidence for an increase in the mass scale of quasar host halos---the inferred median mass of halos hosting central quasars at z~3.2 is Mcen=(14.1+5.8/-6.9) 10^12 Msun/h. The cutoff profiles of the mean occupation functions of central quasars reveal that quasar luminosity is more tightly correlated with halo mass at higher redshifts. The average quasar duty cycle around the median host halo mass is inferred to be fq=(7.3+0.6/-1.5) 10^(-4) at z~1.4 and fq=(8.6+20.4/-7.2) 10^(-2) at z~3.2. We discuss the implications of our results for quasar evolution and quasar-galaxy co-evolution.

preprint2012arXiv

The Ninth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Baryon Oscillation Spectroscopic Survey

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

preprint2012arXiv

The SDSS-III Baryon Oscillation Spectroscopic Survey: The Quasar Luminosity Function from Data Release Nine

We present a new measurement of the optical Quasar Luminosity Function (QLF), using data from the Sloan Digital Sky Survey-III: Baryon Oscillation Spectroscopic Survey (SDSS-III: BOSS). From the SDSS-III Data Release Nine (DR9), we select a uniform sample of 22,301 i<=21.8 quasars over an area of 2236 sq. deg with confirmed spectroscopic redshifts between 2.2<z<3.5, filling in a key part of the luminosity-redshift plane for optical quasar studies. We derive the completeness of the survey through simulated quasar photometry, and check this completeness estimate using a sample of quasars selected by their photometric variability within the BOSS footprint. We investigate the level of systematics associated with our quasar sample using the simulations, in the process generating color-redshift relations and a new quasar k-correction. We probe the faint end of the QLF to M_i(z=2.2) = -24.5 and see a clear break in the QLF at all redshifts up to z=3.5. We find that a log-linear relation (in log[Phi*] - M*) for a luminosity and density evolution (LEDE) model adequately describes our data within the range 2.2<z<3.5; across this interval the break luminosity increases by a factor of ~2.3 while Phi* declines by a factor of ~6. At z<2.2 our data is reasonably well fit by a pure luminosity evolution (PLE) model. We see only a weak signature of "AGN downsizing", in line with recent studies of the hard X-ray luminosity function. We compare our measured QLF to a number of theoretical models and find that models making a variety of assumptions about quasar triggering and halo occupation can fit our data over a wide range of redshifts and luminosities.

preprint2012arXiv

The Sloan Digital Sky Survey quasar catalog: ninth data release

We present the Data Release 9 Quasar (DR9Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the survey, are spectrocopically confirmed as quasars via visual inspection, have luminosities Mi[z=2]<-20.5 (in a $Λ$CDM cosmology with H0 = 70 km/s/Mpc, $Ω_{\rm M}$ = 0.3, and $Ω_Λ$ = 0.7) and either display at least one emission line with full width at half maximum (FWHM) larger than 500 km/s or, if not, have interesting/complex absorption features. It includes as well, known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. This catalog contains 87,822 quasars (78,086 are new discoveries) detected over 3,275 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra newly derived from a training set of 8,632 spectra from SDSS-DR7. The number of quasars with $z>2.15$ (61,931) is ~2.8 times larger than the number of z>2.15 quasars previously known. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 7,533 broad absorption line quasars and gives their characteristics. For each object the catalog presents five-band (u,g,r,i,z) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

preprint2011arXiv

A Catalog of Quasar Properties from SDSS DR7

We present a compilation of properties of the 105,783 quasars in the SDSS Data Release 7 (DR7) quasar catalog. In this value-added product, we compile continuum and emission line measurements around the Halpha, Hbeta, MgII and CIV regions, as well as other quantities such as radio properties, broad absorption line quasar (BALQSO) flags, and disk emitters. We also compile virial black hole mass estimates based on various calibrations. For the fiducial virial mass estimates we use the Vestergaard & Peterson (VP06) calibrations for Hbeta and CIV, and our own calibration for MgII which matches the VP06 Hbeta masses on average. We describe the construction of this catalog, and discuss its limitations. The catalog and its future updates will be made publicly available online.

preprint2011arXiv

Active Galactic Nucleus Pairs from the Sloan Digital Sky Survey. I. The Frequency on ~ 5-100 kpc Scales

Galaxy-galaxy mergers and close interactions have long been regarded as a viable mechanism for channeling gas toward the central supermassive black holes (SMBHs) of galaxies which are triggered as active galactic nuclei (AGNs). AGN pairs, in which the central SMBHs of a galaxy merger are both active, are expected to be common from such events. We conduct a systematic study of 1286 AGN pairs at z ~ 0.1 with line-of-sight velocity offsets dv < 600 km/s and projected separations r_p < 100 kpc, selected from the Seventh Data Release of the Sloan Digital Sky Survey (SDSS). This AGN pair sample was drawn from 138,070 AGNs optically identified based on diagnostic emission line ratios and/or line widths. The fraction of AGN pairs with 5 kpc ~< r_p < 100 kpc among all spectroscopically selected AGNs at 0.02 < z < 0.16 is 3.6% after correcting for SDSS spectroscopic incompleteness; ~ 30% of these pairs show morphological tidal features in their SDSS images, and the fraction becomes >~ 80% for pairs with the brightest nuclei. Our sample increases the number of known AGN pairs on these scales by more than an order of magnitude. We present their AGN and host-galaxy star formation properties in a companion paper.

preprint2011arXiv

Active Galactic Nucleus Pairs from the Sloan Digital Sky Survey. II. Evidence for Tidally Enhanced Star Formation and Black Hole Accretion

Active galactic nuclei (AGNs) are occasionally seen in pairs, suggesting that tidal encounters are responsible for the accretion of material by both central supermassive black holes (BHs). In Paper I of this series, we selected a sample of AGN pairs with projected separations r_p < 100 kpc and velocity offsets < 600 km/s from the SDSS DR7 and quantified their frequency. In this paper, we address the BH accretion and recent star formation properties in their host galaxies. AGN pairs experience stronger BH accretion, as measured by their [O III]5007 luminosities (corrected for contribution from star formation) and Eddington ratios, than do control samples of single AGNs matched in redshift and host stellar mass. Their host galaxies have stronger post-starburst activity and younger mean stellar ages, as indicated by stronger H-delta absorption and smaller 4000 A break in their spectra. The BH accretion and recent star formation in the host galaxies both increase with decreasing projected separation in AGN pairs, for r_p ~< 10-30 kpc. The intensity of BH accretion, the post-starburst strength, and the mean stellar ages are correlated between the two AGNs in a pair. The luminosities and Eddington ratios of AGN pairs are correlated with recent star formation in their host galaxies, with a scaling relation consistent with that observed in single AGNs. Our results suggest that galaxy tidal interactions enhance both BH accretion and star formation in AGN pairs, even though the majority of low redshift AGNs is not coincident with on-going interactions.

preprint2011arXiv

Cosmic Train Wreck by Massive Black Holes: Discovery of a kpc-Scale Triple Active Galactic Nucleus

Hierarchical galaxy mergers will lead to the formation of binary and, in the case of a subsequent merger before a binary coalesce, triple supermassive black holes (SMBHs), given that most massive galaxies harbor SMBHs. A triple of SMBHs becomes visible as a triple active galactic nucleus (AGN) when the BHs accrete large amounts of gas at the same time. Here we report the discovery of a kpc-scale triple AGN, SDSSJ1027+1749 at z = 0.066, from our systematic search for hierarchical mergers of AGNs. The galaxy contains three emission-line nuclei, two of which are offset by ~ 450 and 110 km/s in velocity and by 2.4 and 3.0 kpc in projected separation from the central nucleus. All three nuclei are classified as obscured AGNs based on optical diagnostic emission line ratios, with black hole mass estimates M_BH ~> 10^8 M_sun from stellar velocity dispersions measured in the associated stellar components. Based on dynamical friction timescale estimates, the three stellar components in SDSSJ1027+1749 will merge in ~ 40 Myr, and their associated SMBHs may evolve into a gravitationally interacting triple system in ~< 200 Myr. Our result sets a lower limit of ~ 5 x 10^(-5) for the fraction of kpc-scale triples in optically selected AGNs at z ~ 0.1.

preprint2011arXiv

Type 2 Active Galactic Nuclei with Double-Peaked [OIII] Lines. II. Single AGNs with Complex Narrow-Line Region Kinematics are More Common than Binary AGNs

(Abridged) Approximately 1% of low redshift (z<0.3) optically-selected type 2 AGNs show a double-peaked [OIII] narrow emission line profile in their spatially-integrated spectra. Such features are usually interpreted as due either to kinematics, such as biconical outflows and/or disk rotation of the narrow line region (NLR) around single black holes, or to the relative motion of two distinct NLRs in a merging pair of AGNs. Here we report follow-up near infrared (NIR) imaging and optical slit spectroscopy of 31 double-peaked [OIII] type 2 AGNs drawn from the SDSS parent sample presented in Liu et al (2010). These data reveal a mixture of origins for the double-peaked feature. Roughly 10% of our objects are best explained by binary AGNs at (projected) kpc-scale separations, where two stellar components with spatially coincident NLRs are seen. ~ 50% of our objects have [OIII] emission offset by a few kpc, corresponding to the two velocity components seen in the SDSS spectra, but there are no corresponding double stellar components seen in the NIR imaging. For those objects with sufficiently high quality slit spectra, we see velocity and/or velocity dispersion gradients in [OIII] emission, suggestive of the kinematic signatures of a single NLR. The remaining ~40% of our objects are ambiguous, and will need higher spatial resolution observations to distinguish between the two scenarios. Our observations therefore favor the kinematics scenario with a single AGN for the majority of these double-peaked [OIII] type 2 AGNs. We emphasize the importance of combining imaging and slit spectroscopy in identifying kpc binary AGNs, i.e., in no cases does one of these alone allow an unambiguous identification. We estimate that ~ 0.5-2.5% of the z<0.3 type 2 AGNs are kpc-scale binary AGNs of comparable luminosities, with a relative orbital velocity >~150 km/s.

preprint2011arXiv

Unification of Luminous Type 1 Quasars through CIV Emission

Using a sample of 30,000 quasars from SDSS-DR7, we explore the range of properties exhibited by high-ionization, broad emission lines, such as CIV 1549. Specifically we investigate the anti-correlation between L_UV and emission line EQW (the Baldwin Effect) and the "blueshifting" of high-ionization emission lines. The blueshift of the CIV emission line is nearly ubiquitous, with a mean shift of 810 km/s for radio-quiet (RQ) quasars and 360 km/s for radio-loud (RL) quasars, and the Baldwin Effect is present in both RQ and RL samples. Composite spectra are constructed as a function of CIV emission line properties in attempt to reveal empirical relationships between different line species and the SED. Within a two-component disk+wind model of the broad emission line region (BELR), where the wind filters the continuum seen by the disk component, we find that RL quasars are consistent with being dominated by the disk component, while BALQSOs are consistent with being dominated by the wind component. Some RQ objects have emission line features similar to RL quasars; they may simply have insufficient black hole (BH) spin to form radio jets. Our results suggest that there could be significant systematic errors in the determination of L_bol and BH mass that make it difficult to place these findings in a more physical context. However, it is possible to classify quasars in a paradigm where the diversity of BELR parameters are due to differences in an accretion disk wind between quasars (and over time); these differences are underlain primarily by the SED, which ultimately must be tied to BH mass and accretion rate.

preprint2010arXiv

Constraints on black hole duty cycles and the black hole-halo relation from SDSS quasar clustering

We use Shen et al.'s (2009) measurements of luminosity-dependent clustering in the SDSS Data Release 5 Quasar Catalog, at redshifts 0.4 < z < 2.5, to constrain the relation between quasar luminosity and host halo mass and to infer the duty cycle f_opt, the fraction of black holes that shine as optically luminous quasars at a given time. We assume a monotonic mean relation between quasar luminosity and host halo mass, with log-normal scatter Σ. For specified f_opt and Σ, matching the observed quasar space density determines the normalization of the luminosity-halo mass relation, from which we predict the clustering bias. The data show no change of bias between the faint and bright halves of the quasar sample but a modest increase in bias for the brightest 10%. At the mean redshift z=1.45 of the sample, the data can be well described either by models with small intrinsic scatter (Σ=0.1 dex) and a duty cycle f_opt=6*10^(-4) or by models with much larger duty cycles and larger values of the scatter. "Continuity equation" models of the black hole mass population imply f_opt > 2*10^(-3) in this range of masses and redshifts, and the combination of this constraint with the clustering measurements implies scatter Σ> 0.4 dex. These findings contrast with those inferred from the much stronger clustering of high-luminosity quasars at z~4, which require minimal scatter between luminosity and halo mass and duty cycles close to one.

preprint2010arXiv

Discovery of Four kpc-Scale Binary AGNs

We report the discovery of four kpc-scale binary AGNs. These objects were originally selected from the Sloan Digital Sky Survey based on double-peaked [O III] 4959,5007 emission lines in their fiber spectra. The double peaks could result from pairing active supermassive black holes (SMBHs) in a galaxy merger, or could be due to bulk motions of narrow-line region gas around a single SMBH. Deep near-infrared (NIR) images and optical slit spectra obtained from the Magellan 6.5 m and the APO 3.5 m telescopes strongly support the binary SMBH scenario for the four objects. In each system, the NIR images reveal tidal features and double stellar components with a projected separation of several kpc, while optical slit spectra show two Seyfert 2 nuclei spatially coincident with the stellar components, with line-of-sight velocity offsets of a few hundred km/s. These objects were drawn from a sample of only 43 objects, demonstrating the efficiency of this technique to find kpc-scale binary AGNs.

preprint2010arXiv

Dust-free quasars in the early Universe

The most distant quasars known, at redshifts z=6, generally have properties indistinguishable from those of lower-redshift quasars in the rest-frame ultraviolet/optical and X-ray bands. This puzzling result suggests that these distant quasars are evolved objects even though the Universe was only seven per cent of its current age at these redshifts. Recently one z=6 quasar was shown not to have any detectable emission from hot dust, but it was unclear whether that indicated different hot-dust properties at high redshift or if it is simply an outlier. Here we report the discovery of a second quasar without hot-dust emission in a sample of 21 z=6 quasars. Such apparently hot-dust-free quasars have no counterparts at low redshift. Moreover, we demonstrate that the hot-dust abundance in the 21 quasars builds up in tandem with the growth of the central black hole, whereas at low redshift it is almost independent of the black hole mass. Thus z=6 quasars are indeed at an early evolutionary stage, with rapid mass accretion and dust formation. The two hot-dust-free quasars are likely to be first-generation quasars born in dust-free environments and are too young to have formed a detectable amount of hot dust around them.

preprint2010arXiv

Identifying Supermassive Black Hole Binaries with Broad Emission Line Diagnosis

Double-peaked broad emission lines in Active Galactic Nuclei (AGNs) may indicate the existence of a bound supermassive black hole (SMBH) binary where two distinct broad line regions (BLRs) contribute together to the line profile. An alternative interpretation is a disk emitter origin for the double-peaked line profile. Using simple BLR models, we calculate the expected broad line profile for a SMBH binary at different separations. Under reasonable assumptions that both BLRs are illuminated by the two active SMBHs and that the ionizing flux at the BLR location is roughly constant, we confirm the emergence of double-peaked features and radial velocity drifts of the two peaks due to the binary orbital motion. However, such a clear double-peaked feature only arises in a particular stage of the binary evolution when the two BHs are close enough such that the line-of-sight orbital velocity difference is larger than the FWHM of the individual broad components, while the two BLRs are still mostly distinct. Prior to this stage, the velocity splitting due to the orbit motion of the binary is too small to separate the emission from the two BLRs, leading to asymmetric broad line profiles in general. When the two BHs are even closer such that the two BLRs can no longer be distinct, the line profile becomes more complex and the splitting of the peaks does not correspond to the orbital motion of the binary. In this regime there are no coherent radial velocity drifts in the peaks with time. Asymmetric line profiles are probably a far more common signature of binary SMBHs than are double-peaked profiles. We discuss the temporal variations of the broad line profile for binary SMBHs and highlight the different behaviors of reverberation mapping in the binary and disk emitter cases, which may serve as a feasible tool to disentangle these two scenarios.

preprint2010arXiv

On the Populations of Radio Galaxies with Extended Morphology at z<0.3

Extended extragalactic radio sources have traditionally been classified into FR I and II types, based on the ratio r of the separation S between the brightest regions on either sides of the host galaxy and the total size T of the radio source. Here we examine the distribution of various physical properties as a function of r of 1040 luminous radio galaxies (RGs) at z<0.3 from the SDSS, NVSS, and FIRST. About 2/3 of the RGs are lobe-dominated (LD), and 1/3 have prominent jets. If we follow the original definition of the FR types (a division based solely on r), FR I and FR II RGs overlap in their host galaxy properties. However, the rare, LD sources with r>0.8 AND OIII5007 line luminosity >10^6 Lsun are markedly different on average from the rest of the RGs, in the sense that they are hosted in lower mass galaxies, live in relatively sparse environments, and have higher accretion rates onto the central SMBH. Thus these objects and the rest of RGs form a well-defined dichotomy. Motivated by the stark differences in the nuclear emission line properties of the RG subsamples, we suggest that the accretion rate onto the SMBH may play the primary role in creating the different morphologies. At relatively high accretion rates, the accretion system may produce powerful jets that create the "classical double" morphology (roughly corresponding to the LD sources with r>0.8 and emission lines); at lower accretion rates the jets from a radiatively inefficient accretion flow generate radio lobes without apparent "hot spots" at the edge (corresponding to the majority of LD sources). At slightly lower accretion rates AND in galaxies with dense galactic structure, sources with prominent jets result. It is possible that while the high accretion rate systems could affect sub-Mpc scale environments, the jets from lower accretion rate systems may efficiently suppress activity within the host galaxies.

preprint2010arXiv

Space density of optically-selected type 2 quasars

Type 2 quasars are luminous active galactic nuclei (AGN) whose central regions are obscured by large amounts of gas and dust. In this paper, we present a catalog of type 2 quasars from the Sloan Digital Sky Survey (SDSS), selected based on their optical emission lines. The catalog contains 887 objects with redshifts z < 0.83; this is six times larger than the previous version and is by far the largest sample of type 2 quasars in the literature. We derive the [OIII]5008 luminosity function for 10^8.3 Lsun < L[OIII] < 10^10 Lsun (corresponding to intrinsic luminosities up to M[2400A]-28 mag or bolometric luminosities up to 4x10^47 erg/sec). This luminosity function provides strong lower limits to the actual space density of obscured quasars, due to our selection criteria, the details of the spectroscopic target selection, as well as other effects. We derive the equivalent luminosity function for the complete sample of type 1 (unobscured) quasars; then, we determine the ratio of type 2/type 1 quasar number densities. Our best data constrain this ratio to be at least 1.5:1 for 10^8.3 Lsun < L[OIII] < 10^9.5 Lsun at z < 0.3, and at least 1.2:1 for L[OIII]=10^10 Lsun at 0.3 < z < 0.83. Type 2 quasars are at least as abundant as type 1 quasars in the relatively nearby Universe (z < 0.8) for the highest luminosities.

preprint2010arXiv

The Impact of the Uncertainty in Single-Epoch Virial Black Hole Mass Estimates on the Observed Evolution of the Black Hole - Bulge Scaling Relations

Recent observations of the black hole (BH) - bulge scaling relations usually report positive redshift evolution, with higher redshift galaxies harboring more massive BHs than expected from the local relations. All of these studies focus on broad line quasars with BH mass estimated from virial estimators based on single-epoch spectra. Since the sample selection is largely based on quasar luminosity, the cosmic scatter in the BH-bulge relation introduces a statistical bias leading to on average more massive BHs given galaxy properties at high redshift (Lauer et al. 2007). We here emphasize a previously under-appreciated statistical bias resulting from the uncertainty of single-epoch virial BH mass estimators and the shape of the underlying (true) BH mass function, which leads to on average overestimation of the true BH masses at the high-mass end (Shen et al. 2008). We demonstrate that the latter virial mass bias can contribute a substantial amount to the observed excess in BH mass at fixed bulge properties, comparable to the Lauer et al. bias. The virial mass bias is independent of the Lauer et al. bias, hence if both biases are at work, they can largely (or even fully) account for the observed BH mass excess at high redshift.

preprint2010arXiv

The Sloan Digital Sky Survey Quasar Catalog V. Seventh Data Release

We present the fifth edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog, which is based upon the SDSS Seventh Data Release. The catalog, which contains 105,783 spectroscopically confirmed quasars, represents the conclusion of the SDSS-I and SDSS-II quasar survey. The catalog consists of the SDSS objects that have luminosities larger than M_i = -22.0 (in a cosmology with H_0 = 70 km/s/Mpc Omega_M = 0.3, and Omega_Lambda = 0.7) have at least one emission line with FWHM larger than 1000 km/s or have interesting/complex absorption features, are fainter than i > 15.0 and have highly reliable redshifts. The catalog covers an area of 9380 deg^2. The quasar redshifts range from 0.065 to 5.46, with a median value of 1.49; the catalog includes 1248 quasars at redshifts greater than four, of which 56 are at redshifts greater than five. The catalog contains 9210 quasars with i < 18; slightly over half of the entries have i< 19. For each object the catalog presents positions accurate to better than 0.1" rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3800-9200 Ang. at a spectral resolution R = 2000 the spectra can be retrieved from the SDSS public database using the information provided in the catalog. Over 96% of the objects in the catalog were discovered by the SDSS. We also include a supplemental list of an additional 207 quasars with SDSS spectra whose archive photometric information is incomplete.

preprint2009arXiv

Binary Quasars at High Redshift I: 24 New Quasar Pairs at z ~ 3-4

The clustering of quasars on small scales yields fundamental constraints on models of quasar evolution and the buildup of supermassive black holes. This paper describes the first systematic survey to discover high redshift binary quasars. Using color-selection and photometric redshift techniques, we searched 8142 deg^2 of SDSS imaging data for binary quasar candidates, and confirmed them with follow-up spectroscopy. Our sample of 27 high redshift binaries (24 of them new discoveries) at redshifts 2.9 < z < 4.3 with proper transverse separations 10 kpc < R_{\perp} < 650 kpc increases the number of such objects known by an order of magnitude. Eight members of this sample are very close pairs with R_{\perp} < 100 kpc, and of these close systems four are at z > 3.5. The completeness and efficiency of our well-defined selection algorithm are quantified using simulated photometry and we find that our sample is ~ 50% complete. Our companion paper uses this knowledge to make the first measurement of the small scale clustering (R < 1 Mpc/h comoving) of high-redshift quasars. High redshift binaries constitute exponentially rare coincidences of two extreme (M >~ 10^9 Msun) supermassive black holes. At z ~ 4 there is about one close binary per 10 Gpc^3, thus these could be the highest sigma peaks, the analogs of superclusters, in the early Universe.

preprint2009arXiv

Binary Quasars at High Redshift II: Sub-Mpc Clustering at z ~ 3-4

We present measurements of the small-scale (0.1<~ r <~ 1 Mpc/h) quasar two-point correlation function at z>2.9, for a flux-limited (i<21) sample of 15 binary quasars compiled by Hennawi et al. (2009). The amplitude of the small-scale clustering increases from z ~ 3 to z ~ 4. The small-scale clustering amplitude is comparable to or lower than power-law extrapolations (with slope gamma=2) from the large-scale correlation function of the i<20.2 quasar sample from the Sloan Digital Sky Survey. Using simple prescriptions relating quasars to dark matter halos, we model the observed small-scale clustering with halo occupation models. Reproducing the large-scale clustering amplitude requires that the active fraction of the black holes in the central galaxies of halos is near unity, but the level of small-scale clustering favors an active fraction of black holes in satellite galaxies 0.1 <~ f_s <~ 0.5 at z >~ 3.

preprint2009arXiv

Clustering of Low-Redshift (z <= 2.2) Quasars from the Sloan Digital Sky Survey

We present measurements of the quasar two-point correlation function, ξ_{Q}, over the redshift range z=0.3-2.2 based upon data from the SDSS. Using a homogeneous sample of 30,239 quasars with spectroscopic redshifts from the DR5 Quasar Catalogue, our study represents the largest sample used for this type of investigation to date. With this redshift range and an areal coverage of approx 4,000 deg^2, we sample over 25 h^-3 Gpc^3 (comoving) assuming the current LCDM cosmology. Over this redshift range, we find that the redshift-space correlation function, xi(s), is adequately fit by a single power-law, with s_{0}=5.95+/-0.45 h^-1 Mpc and γ_{s}=1.16+0.11-0.16 when fit over s=1-25 h^-1 Mpc. Using the projected correlation function we calculate the real-space correlation length, r_{0}=5.45+0.35-0.45 h^-1 Mpc and γ=1.90+0.04-0.03, over scales of rp=1-130 h^-1 Mpc. Dividing the sample into redshift slices, we find very little, if any, evidence for the evolution of quasar clustering, with the redshift-space correlation length staying roughly constant at s_{0} ~ 6-7 h^-1 Mpc at z<2.2 (and only increasing at redshifts greater than this). Comparing our clustering measurements to those reported for X-ray selected AGN at z=0.5-1, we find reasonable agreement in some cases but significantly lower correlation lengths in others. We find that the linear bias evolves from b~1.4 at z=0.5 to b~3 at z=2.2, with b(z=1.27)=2.06+/-0.03 for the full sample. We compare our data to analytical models and infer that quasars inhabit dark matter haloes of constant mass M ~2 x 10^12 h^-1 M_Sol from redshifts z~2.5 (the peak of quasar activity) to z~0. [ABRIDGED]

preprint2009arXiv

High-Redshift SDSS Quasars with Weak Emission Lines

We identify a sample of 74 high-redshift quasars (z>3) with weak emission lines from the Fifth Data Release of the Sloan Digital Sky Survey and present infrared, optical, and radio observations of a subsample of four objects at z>4. These weak emission-line quasars (WLQs) constitute a prominent tail of the Lya+NV equivalent width distribution, and we compare them to quasars with more typical emission-line properties and to low-redshift active galactic nuclei with weak/absent emission lines, namely BL Lac objects. We find that WLQs exhibit hot (T~1000 K) thermal dust emission and have rest-frame 0.1-5 micron spectral energy distributions that are quite similar to those of normal quasars. The variability, polarization, and radio properties of WLQs are also different from those of BL Lacs, making continuum boosting by a relativistic jet an unlikely physical interpretation. The most probable scenario for WLQs involves broad-line region properties that are physically distinct from those of normal quasars.

preprint2009arXiv

Modeling Porous Dust Grains with Ballistic Aggregates. II. Light Scattering Properties

We study the light scattering properties of random ballistic aggregates constructed in Shen et al. (Paper I). Using the discrete-dipole-approximation, we compute the scattering phase function and linear polarization for random aggregates with various sizes and porosities, and with two different compositions: 100% silicate and 50% silicate-50% graphite. We investigate the dependence of light scattering properties on wavelength, cluster size and porosity using these aggregate models. We find that while the shape of the phase function depends mainly on the size parameter of the aggregates, the linear polarization depends on both the size parameter and the porosity of the aggregates, with increasing degree of polarization as the porosity increases. Contrary to previous studies, we argue that monomer size has negligible effects on the light scattering properties of ballistic aggregates, as long as the constituent monomer is smaller than the incident wavelength up to 2*pi*a_0/lambda\sim 1.6 where a_0 is the monomer radius. Previous claims for such monomer size effects are in fact the combined effects of size parameter and porosity. Finally, we present aggregate models that can reproduce the phase function and polarization of scattered light from the AU Mic debris disk and from cometary dust, including the negative polarization observed for comets at scattering angles 160<theta<180 deg. These aggregates have moderate porosities, P\sim 0.6, and are of sub-micron-size for the debris disk case, or micron-size for the comet case.

preprint2009arXiv

Type 2 AGNs with Double-Peaked [O III] Lines: Narrow Line Region Kinematics or Merging Supermassive Black Hole Pairs?

We present a sample of 167 type 2 AGNs with double-peaked [O III] 4959,5007 narrow emission lines, selected from the Seventh Data Release of the Sloan Digital Sky Survey. The double-peaked profiles can be well modeled by two velocity components, blueshifted and redshifted from the systemic velocity. Half of these objects have a more prominent redshifted component. In cases where the H-beta emission line is strong, it also shows two velocity components whose line-of-sight (LOS) velocity offsets are consistent with those of [O III]. The relative LOS velocity offset between the two components is typically a few hundred km/s, larger by a factor of ~ 1.5 than the line full width at half maximum of each component. The offset correlates with the host stellar velocity dispersion sigma*. The host galaxies of this sample show systematically larger sigma*, stellar masses, and concentrations, and older luminosity-weighted mean stellar ages than a regular type 2 AGN sample matched in redshift, [O III] 5007 equivalent width and luminosity; they show no significant difference in radio properties. These double-peaked features could be due to narrow-line region kinematics, or binary black holes. The statistical properties do not show strong preference for or against either scenario, and spatially resolved optical imaging, spectroscopy, radio or X-ray followup are needed to draw firm conclusions.

preprint2008arXiv

Quasar Clustering from SDSS DR5: Dependences on Physical Properties

Using a homogenous sample of 38,208 quasars with a sky coverage of $4000 {\rm deg^2}$ drawn from the SDSS Data Release Five quasar catalog, we study the dependence of quasar clustering on luminosity, virial black hole mass, quasar color, and radio loudness. At $z<2.5$, quasar clustering depends weakly on luminosity and virial black hole mass, with typical uncertainty levels $\sim 10%$ for the measured correlation lengths. These weak dependences are consistent with models in which substantial scatter between quasar luminosity, virial black hole mass and the host dark matter halo mass has diluted any clustering difference, where halo mass is assumed to be the relevant quantity that best correlates with clustering strength. However, the most luminous and most massive quasars are more strongly clustered (at the $\sim 2σ$ level) than the remainder of the sample, which we attribute to the rapid increase of the bias factor at the high-mass end of host halos. We do not observe a strong dependence of clustering strength on quasar colors within our sample. On the other hand, radio-loud quasars are more strongly clustered than are radio-quiet quasars matched in redshift and optical luminosity (or virial black hole mass), consistent with local observations of radio galaxies and radio-loud type 2 AGN. Thus radio-loud quasars reside in more massive and denser environments in the biased halo clustering picture. Using the Sheth et al.(2001) formula for the linear halo bias, the estimated host halo mass for radio-loud quasars is $\sim 10^{13} h^{-1}M_\odot$, compared to $\sim 2\times 10^{12} h^{-1}M_\odot$ for radio-quiet quasar hosts at $z\sim 1.5$.

preprint2006arXiv

Differences in the AGN Populations of Groups and Clusters: Clues to AGN Evolution

We combine optical and X-ray data for eight low redshift ($z\sim 0.06$) poor groups of galaxies from the {\it XI}({\it XMM/IMACS}) Groups Project to study the AGN population in the group environment. Among $\sim 140$ group members, we identify five AGN based on their optical emission lines. None of these optically-selected AGN are detected by {\it XMM-Newton}. One additional AGN is discovered in the {\it XMM-Newton} observations. This X-ray detected AGN, which has no obvious AGN emission line signatures in its optical spectrum, is a member of the only X-ray luminous group in our sample. The lack of a significant population of X-ray bright, but optically dull AGN among less dynamically evolved groups is in stark contrast to the large fraction of such objects in rich clusters of galaxies (Martini et al. 2006). We suggest this result can be explained by a physical scenario for AGN accretion evolution: AGN activity is initially triggered by galaxy merging, leading to a high accretion rate and an optically dominant phase (via thin disk accretion). As the accretion rate drops in time, the AGN gradually enters an X-ray dominant low-accretion phase (via a radiative inefficient accretion flow). In this picture, optical- and X-ray-selected AGN are the same population of supermassive black holes observed at different epochs. Within the context of this scenario, the majority of AGN in poor groups are in the high-accretion optically dominant phase, while the AGN population in rich clusters is mostly in the low-accretion X-ray dominant phase.

preprint2006arXiv

Three Dimensional Compressible Hydrodynamic Simulations of Vortices in Disks

We carry out three-dimensional, high resolution (up to $1024^2\times 256$) hydrodynamic simulations of the evolution of vortices in vertically unstratified Keplerian disks using the shearing sheet approximation. The transient amplification of incompressible, linear amplitude leading waves (which has been proposed as a possible route to nonlinear hydrodynamical turbulence in disks) is used as one test of our algorithms; our methods accurately capture the predicted amplification, converges at second-order, and is free from aliasing. Waves expected to reach nonlinear amplitude at peak amplification become unstable to Kelvin-Helmholtz modes when $\mid W_{\rm max}\mid\gtrsim Ω$ (where $W_{\rm max}$ is the local maximum of vorticity and $Ω$ the angular velocity). We study the evolution of a power-law distribution of vorticity consistent with Kolmogorov turbulence; in two-dimensions long-lived vortices emerge and decay slowly, similar to previous studies. In three-dimensions, however, vortices are unstable to bending modes, leading to rapid decay. Only vortices with a length to width ratio smaller than one survive; in three-dimensions the residual kinetic energy and shear stress is at least one order of magnitude smaller than in two-dimensions. No evidence for sustained hydrodynamical turbulence and transport is observed in three-dimensions. Instead, at late times the residual transport is determined by the amplitude of slowly decaying, large-scale vortices (with horizontal extent comparable to the scale height of the disk), with additional contributions from nearly incompressible inertial waves possible. Evaluating the role that large-scale vortices play in astrophysical accretion disks will require understanding the mechanisms that generate and destroy them.

preprint2004arXiv

Global axisymmetric stability analysis for a composite system of two gravitationally coupled scale-free discs

In a composite system of gravitationally coupled stellar and gaseous discs, we perform linear stability analysis for axisymmetric coplanar perturbations using the two-fluid formalism. The background stellar and gaseous discs are taken to be scale-free with all physical variables varying as powers of cylindrical radius $r$ with compatible exponents. The unstable modes set in as neutral modes or stationary perturbation configurations with angular frequency $ω=0$.

preprint2004arXiv

Outflows and inflows in astrophysical systems

We seek for self-similar solutions describing the time-dependent evolution of self-gravity systems with either spherical symmetry or axisymmetric disk geometry. By assuming self-similar variable $x\equiv r/at$ where $a$ is isothermal sound speed we find self-similar solutions extending from the initial instant $t=0$ to the final stage $t\to \infty$ using standard semi-analytical methods. Different types of solutions are constructed, which describe overall expansion or collapse, envelope expansion with core collapse (EECC), the formation of central rotationally supported quasi-equilibrium disk as well as shocks. Though infinitely many, these self-similarity solutions have similar asymptotic behaviors which may impose diagnosis on the velocity and density structures in astrophysical systems.