Search

Quasar outflows may play a crucial role in regulating the host galaxy, although the spatial scale of quasar outflows remain a major enigma, with their acceleration mechanism poorly understood. The kinematic information of outflow is the key to understanding its origin and acceleration mechanism. Here, we report the galactocentric distances of different outflow components for both a sample and an individual quasar. We find that the outflow distance increases with velocity, with a typical value from several parsecs to more than one hundred parsecs, providing direct evidence for an acceleration happening at a scale of the order of 10 parsecs. These outflows carry ~1% of the total quasar energy, while their kinematics are consistent with a dust driven model with a launching radius comparable to the scale of a dusty torus, indicating that the coupling between dust and quasar radiation may produce powerful feedback that is crucial to galaxy evolution., 10 pages, 4 figures in main text. Published in Science Advances on 11 Feb 2022

The nearby dwarf spiral galaxy NGC 4395 contains a broad-lined active galactic nucleus (AGN) of exceptionally low luminosity powered by accretion onto a central black hole of very low mass ($\sim10^4-10^5$ M$_\odot$). In order to constrain the size of the optical continuum emission region through reverberation mapping, we carried out high-cadence photometric monitoring of NGC 4395 in the $griz$ filter bands on two consecutive nights in 2022 April using the four-channel MuSCAT3 camera on the Faulkes Telescope North at Haleakal\={a} Observatory. Correlated variability across the $griz$ bands is clearly detected, and the $r$, $i$, and $z$ band light curves show lags of $8.4^{+1.0}_{-1.1}$, $14.2^{+1.2}_{-1.4}$, and $20.4^{+2.0}_{-2.1}$ minutes with respect to the $g$ band when measured using the full-duration light curves. When lags are measured for each night separately, the Night 2 data exhibit lower cross-correlation amplitudes and shorter lags than the Night 1 light curves. Using the full-duration lags, we find that the lag-wavelength relationship is consistent with the $\tau\propto\lambda^{4/3}$ dependence found for more luminous AGN. Combining our results with continuum lags measured for other objects, the lag between $g$ and $z$ band scales with optical continuum luminosity as $\tau_{gz} \propto L^{0.56\pm0.05}$, similar to the scaling of broad-line region size with luminosity, reinforcing recent evidence that diffuse continuum emission from the broad-line region may contribute substantially to optical continuum variability and reverberation lags., Comment: 8 pages, 3 figures, 2 tables, submitted to ApJL. Comments welcome

Extremely variable quasars (EVQs) are a population of sources showing large optical photometric variability revealed by time-domain surveys. The physical origin of such extreme variability is yet unclear. In this first paper of a series, we construct the largest-ever sample of 14,012 EVQs using photometric data spanning over $>$ 15 years from SDSS and Pan-STARRS1. We divide them into five sub-samples according to the relative brightness of an EVQ during SDSS spectroscopic observation compared to the mean brightness from photometric observations. Corresponding control samples of normal quasars are built with matched redshift, bolometric luminosity and supermassive black hole mass. We obtain the composite SDSS spectra of EVQs in various states and their corresponding control samples. We find EVQs exhibit clearly bluer (redder) SDSS spectra during bright (dim) states, consistent with the "bluer-when-brighter" trend widely seen in normal quasars. We further find the line EWs of broad MgII, CIV and OIII (but not broad Hb which is yet puzzling) gradually decrease from dim state to bright state, similar to the so-called intrinsic Baldwin effect commonly seen in normal AGNs. Meanwhile, EVQs have systematically larger line EWs compared with the control samples. We also see that EVQs exhibit subtle excess in the very broad line component compared with control samples. Possible explanations for the discoveries are discussed. Our findings support the hypothesis that EVQs are in the tail of a broad distribution of quasar properties, but not a distinct population., 25 pages and 21 figures; Accepted for publication in ApJ

Photoionization modeling of active galactic nuclei (AGN) predicts that diffuse continuum (DC) emission from the broad-line region makes a substantial contribution to the total continuum emission from ultraviolet through near-infrared wavelengths. Evidence for this DC component is present in the strong Balmer jump feature in AGN spectra, and possibly from reverberation measurements that find longer lags than expected from disk emission alone. However, the Balmer jump region contains numerous blended emission features, making it difficult to isolate the DC emission strength. In contrast, the Paschen jump region near 8200 \r{A} is relatively uncontaminated by other strong emission features. Here, we examine whether the Paschen jump can aid in constraining the DC contribution, using Hubble Space Telescope STIS spectra of six nearby Seyfert 1 nuclei. The spectra appear smooth across the Paschen edge, and we find no evidence of a Paschen spectral break or jump in total flux. We fit multi-component spectral models over the range $6800-9700$ \r{A} and find that the spectra can still be compatible with a significant DC contribution if the DC Paschen jump is offset by an opposite spectral break resulting from blended high-order Paschen emission lines. The fits imply DC contributions ranging from $\sim$10% to 50% at 8000 \r{A}, but the fitting results are highly dependent on assumptions made about other model components. These degeneracies can potentially be alleviated by carrying out fits over a broader wavelength range, provided that models can accurately represent the disk continuum shape, Fe II emission, high-order Balmer line emission, and other components., Comment: 17 pages, 6 figures, and 3 tables; ApJ in press

We perform a systematic survey of active galactic nuclei (AGNs) continuum lags using ∼3 days cadence gri-band light curves from the Zwicky Transient Facility. We select a sample of 94 type 1 AGNs at z < 0.8 with significant and consistent inter-band lags based on the interpolated cross-correlation function method and the Bayesian method JAVELIN. Within the framework of the “lamp-post” reprocessing model, our findings are: (1) The continuum emission (CE) sizes inferred from the data are larger than the disk sizes predicted by the standard thin-disk model. (2) For a subset of the sample, the CE size exceeds the theoretical limit of the self-gravity radius (12 lt-days) for geometrically thin disks. (3) The CE size scales with continuum luminosity as R CE ∝ L 0.48±0.04 with a scatter of 0.2 dex, analogous to the well-known radius–luminosity relation of broad Hβ. These findings suggest a significant contribution of diffuse continuum emission from the broad-line region (BLR) to AGN continuum lags. We find that the R CE–L relation can be explained by a photoionization model that assumes ∼23% of the total flux comes from the diffuse BLR emission. In addition, the ratio of the CE size and model-predicted disk size anticorrelates with the continuum luminosity, which is indicative of a potential nondisk BLR lag contribution evolving with the luminosity. Finally, a robust positive correlation between the CE size and black hole mass is detected.

Periodic quasars have been suggested as candidates for hosting binary supermassive black holes (SMBHs), although alternative scenarios remain possible to explain the optical light curve periodicity. To test the alternative hypothesis of precessing radio jet, we present deep 6 GHz radio imaging conducted with NSF's Karl G. Jansky Very Large Array (VLA) in its C configuration for the three candidate periodic quasars, DES J024703.24$-$010032.0, DES J024944.66$-$000036.8, and DES J025214.67$-$002813.7. Our targets were selected based on their optical variability using 20-yr long multi-color light curves from the Dark Energy Survey (DES) and the Sloan Digital Sky Survey (SDSS). The new VLA observations show that all three periodic quasars are radio-quiet with the radio loudness parameters measured to be $R\equiv f_{6\,{\rm cm}}/f_{{\rm 2500}}$ of $\lesssim$1.0$-$1.5 and the $k$-corrected luminosities $\nu L_\nu$[6 GHz] of $\lesssim$5$-$21 $\times$ 10$^{39}$ erg s$^{-1}$. They are in stark contrast to previously known periodic quasars proposed as binary SMBH candidates such as the blazar OJ 287 and PG1302$-$102. Our results rule out optical emission contributed from precessing radio jets as the origin of the optical periodicity in the three DES$-$SDSS-selected candidate periodic quasars. Future continued optical monitoring and complementary multi-wavelength observations are still needed to further test the binary SMBH hypothesis as well as other competing scenarios to explain the optical periodicity., Comment: 8 pages, 5 figures, 2 tables, accepted for publication in MNRAS

We report the results of near-infrared spectroscopic observations of 37 quasars in the redshift range $6.3< z\le7.64$, including 32 quasars at $z>6.5$, forming the largest quasar near-infrared spectral sample at this redshift. The spectra, taken with Keck, Gemini, VLT, and Magellan, allow investigations of central black hole mass and quasar rest-frame ultraviolet spectral properties. The black hole masses derived from the MgII emission lines are in the range $(0.3-3.6)\times10^{9}\,M_{\odot}$, which requires massive seed black holes with masses $\gtrsim10^{3-4}\,M_{\odot}$, assuming Eddington accretion since $z=30$. The Eddington ratio distribution peaks at $\lambda_{\rm Edd}\sim0.8$ and has a mean of 1.08, suggesting high accretion rates for these quasars. The CIV - MgII emission line velocity differences in our sample show an increase of CIV blueshift towards higher redshift, but the evolutionary trend observed from this sample is weaker than the previous results from smaller samples at similar redshift. The FeII/MgII flux ratios derived for these quasars up to $z=7.6$, compared with previous measurements at different redshifts, do not show any evidence of strong redshift evolution, suggesting metal-enriched environments in these quasars. Using this quasar sample, we create a quasar composite spectrum for $z>6.5$ quasars and find no significant redshift evolution of quasar broad emission lines and continuum slope, except for a blueshift of the CIV line. Our sample yields a strong broad absorption line quasar fraction of $\sim$24%, higher than the fractions in lower redshift quasar samples, although this could be affected by small sample statistics and selection effects., Comment: Accepted for publication in ApJ. 26 pages, 9 figures

We present an analysis of the variability of broad absorption lines (BALs) in a quasar SDSS J141955.26+522741.1 at $z=2.145$ with 72 observations from the Sloan Digital Sky Survey Data Release 16 (SDSS DR16). The strong correlation between the equivalent widths of BAL and the continuum luminosity, reveals that the variation of BAL trough is dominated by the photoionization. The photoionization model predicts that when the time interval $\Delta T$ between two observations is longer than the recombination timescale $t_{\rm rec}$, the BAL variations can be detected. This can be characterized as a "sharp rise" in the detection rate of BAL variation at $\Delta T=t_{\rm rec}$. For the first time, we detect such a "sharp rise" signature in the detection rate of BAL variations. As a result, we propose that the $t_{\rm rec}$ can be obtained from the "sharp rise" of the detection rate of BAL variation. It is worth mentioning that the BAL variations are detected at the time-intervals less than the $t_{\rm rec}$ for half an order of magnitude in two individual troughs. This result indicates that there may be multiple components with different $t_{\rm rec}$ but the same velocity in an individual trough., Comment: 7 pages, 5 figures, Accepted for Publication in ApJL

We report on continued, $\sim$15 year-long, broad Balmer emission lines in three metal-poor dwarf emission-line galaxies selected from Sloan Digital Sky Survey spectroscopy. The persistent luminosity of the broad Balmer emission indicates the galaxies are active galactic nuclei (AGNs) with virial black hole masses of $\sim 10^{6.7}{-10^{7.0}\ M_{\odot}}$. The lack of observed hard X-ray emission and the possibility that the Balmer emission could be due to a long-lived stellar transient motivated additional follow-up spectroscopy. We also identify a previously-unreported blueshifted narrow absorption line in the broad H$\alpha$ feature in one of the AGNs, indicating an AGN-driven outflow with hydrogen column densities of order $10^{17}$ cm$^{-2}$. We also extract light curves from the Catalina Real-Time Transient Survey and the Zwicky Transient Facility. We detect probable AGN-like variability in three galaxies, further supporting the AGN scenario. This also suggests the AGNs are not strongly obscured. This sample of galaxies are among the most metal-poor which host an AGN ($Z=0.05 - 0.16\ Z_\odot$). We speculate they may be analogues to seed black holes which formed in unevolved galaxies at high redshift. Given the rarity of metal-poor AGNs and small sample size available, we investigate prospects for their identification in future spectroscopic and photometric surveys., Comment: 8 pages, 4 figures, 1 table, 1 appendix; accepted for publication in MNRAS

Binary supermassive black holes (BSBHs) are expected to be a generic byproduct from hierarchical galaxy formation. The final coalescence of BSBHs is thought to be the loudest gravitational wave (GW) siren, yet no confirmed BSBH is known in the GW-dominated regime. While periodic quasars have been proposed as BSBH candidates, the physical origin of the periodicity has been largely uncertain. Here, we report discovery of a periodicity (p = 1607 ± 7 d) at 99.95 per cent significance (with a global p value of ∼10-3 accounting for the look elsewhere effect) in the optical light curves of a redshift 1.53 quasar, SDSS J025214.67-002813.7. Combining archival Sloan Digital Sky Survey data with new, sensitive imaging from the Dark Energy Survey, the total ∼20-yr time baseline spans ∼4.6 cycles of the observed 4.4-yr (rest frame 1.7-yr) periodicity. The light curves are best fit by a bursty model predicted by hydrodynamic simulations of circumbinary accretion discs. The periodicity is likely caused by accretion rate modulation by a milli-parsec BSBH emitting GWs, dynamically coupled to the circumbinary accretion disc. A bursty hydrodynamic variability model is statistically preferred over a smooth, sinusoidal model expected from relativistic Doppler boost, a kinematic effect proposed for PG1302-102. Furthermore, the frequency dependence of the variability amplitudes disfavours Doppler boost, lending independent support to the circumbinary accretion variability hypothesis. Given our detection rate of one BSBH candidate from circumbinary accretion variability out of 625 quasars, it suggests that future large, sensitive synoptic surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time may be able to detect hundreds to thousands of candidate BSBHs from circumbinary accretion with direct implications for Laser Interferometer Space Antenna., Liao, Wei-Ting, et al. DES Collaboration, The DES data management system is supported by the National Science Foundation under grant numbers AST-1138766 and AST-1536171. The DES participants from Spanish institutions are partially supported by MINECO under grants AYA2015-71825, ESP2015-66861, FPA2015-68048, SEV-2016-0588, SEV-2016-0597, and MDM-2015-0509, some of which include ERDF funds from the European Union. IFAE is partially funded by the CERCA program of the Generalitat de Catalunya. Research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Program (FP7/2007-2013) including ERC grant agreements 240672, 291329, and 306478. We acknowledge support from the Australian Research Council Centre of Excellence for All-sky Astrophysics (CAASTRO), through project number CE110001020, and the Brazilian Instituto Nacional de Ciênciae Tecnologia (INCT) e-Universe (CNPq grant 465376/2014-2).

We perform a systematic search for high-redshift ($z >$ 1.5) extreme variability quasars (EVQs) using repeat spectra from the Sixteenth Data Release of Sloan Digital Sky Survey, which provides a baseline spanning up to $\sim$18 yrs in the observed frame. We compile a sample of 348 EVQs with a maximum continuum variability at rest frame 1450 Angstrom of more than 100% (i.e., $\delta$V $\equiv$ (Max$-$Min)/Mean $>$1). The EVQs show a range of emission line variability, including 23 where at least one line in our redshift range disappears below detectability, which can then be seen as analogous to low-redshift changing-look quasars (CLQs)". Importantly, spurious CLQs caused by SDSS problematic spectral flux calibration, e.g., fiber drop issue, have been rejected. The similar properties (e.g., continuum/line, difference-composite spectra and Eddington ratio) of normal EVQs and CLQs, implies that they are basically the same physical population with analogous intrinsic variability mechanisms, as a tail of a continuous distribution of normal quasar properties. In addition, we find no reliable evidence ($\lesssim$ 1$\sigma$) to support that the CLQs are a subset of EVQs with less efficient accretion. Finally, we also confirm the anti-breathing of C IV (i.e., line width increases as luminosity increases) in EVQs, and find that in addition to $\sim$ 0.4 dex systematic uncertainty in single-epoch C IV virial black hole mass estimates, an extra scatter of $\sim$ 0.3 dex will be introduced by extreme variability., Comment: 25 pages, 15 figures, 6 tables, published in ApJ

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$\beta$ 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$\beta$ 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$\sigma$ 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., Comment: 13 pages, 8 figures, 1 table, accepted to MNRAS

DES Collaboration: et al., Periodically variable quasars have been suggested as close binary supermassive black holes. We present a systematic search for periodic light curves in 625 spectroscopically confirmed quasars with a median redshift of 1.8 in a 4.6 deg2 overlapping region of the Dark Energy Survey Supernova (DES-SN) fields and the Sloan Digital Sky Survey Stripe 82 (SDSS-S82). Our sample has a unique 20-yr long multicolour (griz) light curve enabled by combining DES-SN Y6 observations with archival SDSS-S82 data. The deep imaging allows us to search for periodic light curves in less luminous quasars (down to r ∼23.5 mag) powered by less massive black holes (with masses ≳ 108.5M⊙) at high redshift for the first time. We find five candidates with significant (at >99.74 per cent single-frequency significance in at least two bands with a global p-value of ∼7 × 10−4–3 × 10−3 accounting for the look-elsewhere effect) periodicity with observed periods of ∼3–5 yr (i.e. 1–2 yr in rest frame) having ∼4–6 cycles spanned by the observations. If all five candidates are periodically variable quasars, this translates into a detection rate of ∼0.8+0.5−0.3 per cent or ∼1.1+0.7−0.5 quasar per deg2. Our detection rate is 4–80 times larger than those found by previous searches using shallower surveys over larger areas. This discrepancy is likely caused by differences in the quasar populations probed and the survey data qualities. We discuss implications on the future direct detection of low-frequency gravitational waves. Continued photometric monitoring will further assess the robustness and characteristics of these candidate periodic quasars to determine their physical origins., The DES data management system is supported by the National Science Foundation under grants AST-1138766 and AST-1536171. The DES participants from Spanish institutions are partially supported by MICINN under grants ESP2017-89838, PGC2018-094773, PGC2018-102021, SEV-2016-0588, SEV-2016-0597, and MDM-2015-0509, some of which include ERDF funds from the European Union. IFAE is partially funded by the CERCA program of the Generalitat de Catalunya. Research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Program (FP7/2007-2013) including ERC grants 240672, 291329, and 306478.

Understanding the origin of \feii\ emission is important because it is crucial to construct the main sequence of Active Galactic Nuclei (AGNs). Despite several decades of observational and theoretical effort, the location of the optical iron emitting region and the mechanism responsible for the positive correlation between the \feii\ strength and the black hole accretion rate remain open questions as yet. In this letter, we report the optical \feii\ response to the central outburst in PS1-10adi, a candidate tidal disruption event (TDE) taking place in an AGN at $z = 0.203$ that has aroused extensive attention. For the first time, we observe that the \feii\ response in the rising phase of its central luminosity is significantly more prominent than that in the decline phase, showing a hysteresis effect. We interpret this hysteresis effect as a consequence of the gradual sublimation of the dust grains situating at the inner surface of the torus into gas when the luminosity of the central engine increases. It is the iron element released from the sublimated dust that contributes evidently to the observed \feii\ emission. This interpretation, together with the weak response of the \hb\ emission as we observe, naturally explains the applicability of relative \feii\ strength as a tracer of the Eddington ratio. In addition, optical iron emission of this origin renders the \feii\ time lag a potential "standard candle" with cosmological implications., Comment: 7 pages, 4 figures, 1 table, Accepted for Publication in The Astrophysical Journal Letters

Quasars whose broad emission lines show temporal, bulk radial velocity (RV) shifts have been proposed as candidate sub-parsec (sub-pc), binary supermassive black holes (BSBHs). We identified a sample of 16 BSBH candidates based on two-epoch spectroscopy among 52 quasars with significant RV shifts over a few rest-frame years. The candidates showed consistent velocity shifts independently measured from two broad lines (H$\beta$ and H$\alpha$ or Mg${\rm \,II}$) without significant changes in the broad-line profiles. Here in the third paper of the series, we present further third- and fourth-epoch spectroscopy for 12 of the 16 candidates for continued RV tests, spanning $\sim$5--15 yr in the quasars' rest frames. Cross-correlation analysis of the broad H$\beta$ calibrated against [O${\rm\,III}]\,\lambda 5007$ suggests that 5 of the 12 quasars remain valid as BSBH candidates. They show broad H$\beta$ RV curves that are consistent with binary orbital motion without significant changes in the broad line profiles. Their broad H$\alpha$ (or Mg${\rm \,II}$) lines display RV shifts that are either consistent with or smaller than those seen in broad H$\beta$. The RV shifts can be explained by a $\sim$0.05--0.1 pc BSBH with an orbital period of $\sim$40--130 yr, assuming a mass ratio of 0.5--2 and a circular orbit. However, the parameters are not well constrained given the few epochs that sample only a small portion of the hypothesized binary orbital cycle. The apparent occurrence rate of sub-pc BSBHs is $\lesssim$13$\pm$5% among all SDSS quasars, with no significant difference in the subsets with and without single-epoch broad line velocity offsets. Dedicated long-term spectroscopic monitoring is still needed to further confirm or reject these BSBH candidates., Comment: 38 pages, 11 figures, MNRAS accepted

We present a study of a sample of 223 radio-loud quasars (up to redshift $, Comment: 11 pages, 11 Figures, accepted for publication in A&A

We report the discovery of SDSS J0849+1114 as the first known triple Type 2 Seyfert nucleus. It represents three active black holes that are identified from new spatially resolved optical slit spectroscopy using the Dual Imaging Spectrograph on the 3.5 m telescope at the Apache Point Observatory. We also present new complementary observations including the Hubble Space Telescope Wide Field Camera 3 U- and Y-band imaging, Chandra Advanced CCD Imaging Spectrometer S-array X-ray 0.5--8 keV imaging spectroscopy, and NSF Karl G. Jansky Very Large Array radio 9.0 GHz imaging in its most extended A configuration. These comprehensive multiwavelength observations, when combined together, strongly suggest that all three nuclei are active galactic nuclei. While they are now still at kiloparsec-scale separations, where the host-galaxy gravitational potential dominates, the black holes may evolve into a bound triple system in $\lesssim$2 Gyr. These triple merger systems may explain the overly massive stellar cores that have been observed in some elliptical galaxies such as M87, which are expected to be unique gravitational wave sources. Similar systems may be more common in the early universe, when galaxy mergers are thought to have been more frequent., 19 pages, 8 figures, 6 tables; ApJ published version

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$\alpha$/H$\beta$ 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$\alpha$/H$\beta$ 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$\alpha$/H$\beta$, and the rare population of broad Mg II emitters in the spectra of massive inactive galaxies., Comment: 15 pages, 10 figures, accepted for publication in ApJ

Changing look active-galatic-nuclei (CL AGNs) can yield considerable insight into accretion physics as well as the co-evolution of black holes and their host galaxies. A large sample of these CL AGNs is essential to achieve the latter goal. We propose an effective method to select CL candidates from spectroscopic quasar catalogs using the mid-infrared (MIR) variability information derived from ALL-WISE/NEOWISE data releases. Our primary selection criteria include both a large amplitude flux variation and a transition of MIR color from an AGN to a normal galaxy. A pilot spectroscopic follow-up of 7 candidates among about 300 candidates selected from Sloan Digital Sky Survey low-redshift (z, 7 pages, 3 figures, 1 table

Periodic quasars are candidates for binary supermassive black holes (BSBHs) efficiently emitting low frequency gravitational waves. Recently, $\sim$150 candidates were identified from optical synoptic surveys. However, they may be false positives caused by stochastic quasar variability given the few cycles covered (typically 1.5). To independently test the binary hypothesis, we search for evidence of truncated or gapped circumbinary accretion discs (CBDs) in their spectral energy distributions (SEDs). Our work is motivated by CBD simulations that predict flux deficits as cutoffs from central cavities opened by secondaries or notches from minidiscs around both BHs. We find that candidate periodic quasars show SEDs similar to those of control quasars matched in redshift and luminosity. While seven of 138 candidates show a blue cutoff in the IR-optical-UV SED, six of which may represent CBDs with central cavities, the red SED fraction is similar to that in control quasars, suggesting no correlation between periodicity and SED anomaly. Alternatively, dust reddening may cause red SEDs. The fraction of extremely radio-loud quasars, e.g., blazars (with $R>100$), is tentatively higher than that in control quasars (at 2.5$\sigma$). Our results suggest that, assuming most periodic candidates are robust, IR-optical-UV SEDs of CBDs are similar to those of accretion discs of single BHs, if the periodicity is driven by BSBHs; the higher blazar fraction may signal precessing radio jets. Alternatively, most current candidate periodic quasars identified from few-cycle light curves may be false positives. Their tentatively higher blazar fraction and lower Eddington ratios may both be caused by selection biases., Comment: 14 pages, 10 figures, 2 tables, accepted for publication in MNRAS

The black hole in the quasar SDSS J140821.67+025733.2 has been reported to have a mass of 1.96 × 1011 M ⊙ based on measurements from the SDSS DR12 Quasar Catalog. As a result, references to this object as the most massive known black hole in the universe have recently appeared in scientific articles and in popular media including Wikipedia. We show that this extremely high mass estimate is incorrect, resulting from a faulty measurement of the C iv line width. From a new measurement of the Mg ii line width, we derive an estimate of M BH ≈ 8 × 109 M ⊙.

The rest-frame UV/optical variability of the quasars in the Sloan Digital Sky Survey (SDSS) Stripe 82 is used to test the Corona-Heated Accretion-disk Reprocessing (CHAR) model of Sun et al. 2020. We adopt our CHAR model and the observed black-hole masses ($M_{\mathrm{BH}}$) and luminosities ($L$) to generate mock light curves that share the same measurement noise and sampling as the real observations. Without any fine-tuning, our CHAR model can satisfactorily reproduce the observed ensemble structure functions for different $M_{\mathrm{BH}}$, $L$, and rest-frame wavelengths. Our analyses reveal that a luminosity-dependent bolometric correction is disfavored over the constant bolometric correction for UV/optical luminosities. Our work demonstrates the possibility of extracting quasar properties (e.g., the bolometric correction or the dimensionless viscosity parameter) by comparing the physical CHAR model with quasar light curves., Comment: 16 pages, 9 figures, accepted to ApJ

We explore the evolution of the time variability (in the optical $g$-band and on timescales of weeks to years) of SDSS Stripe 82 quasars along the quasar main sequence. A parent sample of $1004$ quasars within $0.5\leq z \leq 0.89$ are used for our statistical studies, we then make subsamples from our parent sample: a subsample of $246$ quasars with similar luminosities, and a subsample of $399$ quasars with similar Rfe (i.e., the ratio of the equivalent width of FeII within $4435$--$4685\ \mathrm{\AA}$ to that of \Hbeta). We find the variability amplitude decreases with luminosity ($L_{\mathrm{bol}}$). The anti-correlation between the variability amplitude and Rfe is weak but statistically significant. The characteristic timescale, $\tau$, correlates mostly with quasar luminosity, its dependence on Rfe is statistically insignificant. After controlling luminosity and Rfe, the high- and low-FWHM samples have similar structure functions. These results support the framework that Rfe is governed by Eddington ratio and FWHM of Hbeta is mostly determined by orientation. We then provide new empirical relations between variability parameters and quasar properties (i.e., luminosity and Rfe). Our new relations are consistent with the scenario that quasar variability is driven by the thermal fluctuations in the accretion disk, $\tau$ seems to correspond to the thermal timescale. From our new relations, we find the short-term variability is mostly sensitive to $L_{\mathrm{bol}}$. Basing on this, we propose that quasar short-term (a few months) variability might be a new type of "Standard Candle" and can be adopted to probe cosmology., Comment: 13 pages, 8 figures, Accepted for publication in ApJ

NGC 5548 has been intensively monitored by the AGN Space Telescope and Optical Reverberation Mapping collaboration. Approximately after half of the light curves, the correlation between the broad emission lines and the lag-corrected ultraviolet continua becomes weak. This anomalous behavior is accompanied by an increase of soft X-ray emission. We propose a simple model to understand this anomalous behavior, i.e., the corona might fall down, thereby increasing the covering fraction of the inner disk. Therefore, X-ray and extreme ultraviolet emission suffer from spectral variations. The ultraviolet continua variations are driven by both X-ray and extreme ultraviolet variations. Consequently, the spectral variability induced by the falling corona would dilute the correlation between the broad emission lines and the ultraviolet continua. Our model can explain many additional observational facts, including the dependence of the anomalous behavior on velocity and ionization energy. We also show that the time lag and correlation between the X-ray and the ultraviolet variations change as NGC 5548 displays the anomalous behavior. The time lag is dramatically longer than the expectation from disk reprocessing if the anomalous behavior is properly excluded. During the anomalous state, the time lag approaches the light-travel timescale of disk reprocessing albeit with a much weaker correlation. We speculate that the time lag in the normal state is caused by reprocessing of the broad line region gas. As NGC 5548 enters the abnormal state, the contribution of the broad line region gas is smaller; the time lag reflects disk reprocessing. We also discuss alternative scenarios., Comment: 10 pages, 6 figures, Accepted for publication in ApJ

Kiloparsec-scale dual active galactic nuclei (AGNs) are active supermassive black hole pairs co-rotating in galaxies with separations of less than a few kpc. Expected to be a generic outcome of hierarchical galaxy formation, their frequency and demographics remain uncertain. We have carried out an imaging survey with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) of AGNs with double-peaked narrow [O III] emission lines. HST/WFC3 offers high image quality in the near-infrared (NIR) to resolve the two stellar nuclei, and in the optical to resolve [O III] from ionized gas in the narrow-line regions. This combination has proven to be key in sorting out alternative scenarios. With HST/WFC3 we are able to explore a new population of close dual AGNs at more advanced merger stages than can be probed from the ground. Here we show that the AGN SDSS J0924+0510, which had previously shown two stellar bulges, contains two spatially distinct [O III] regions consistent with a dual AGN. While we cannot completely exclude cross-ionization from a single central engine, the nearly equal ratios of [O III] strongly suggest a dual AGN with a projected angular separation of 0."4, corresponding to a projected physical separation of $r_p$ = 1 kpc at redshift z = 0.1495. This serves as a proof of principle for combining high-resolution NIR and optical imaging to identify close dual AGNs. Our result suggests that studies based on low-resolution and/or low-sensitivity observations may miss close dual AGNs and thereby may underestimate their occurrence rate on $\lesssim$ kpc scales., 11 pages, 8 figures; ApJ in press

Studies have shown that UV/optical light curves of quasars can be described with the prevalent damped random walk (DRW, also known as Ornstein-Uhlenbeck process) model. A white noise power spectral density (PSD) is expected at low frequency in this model, however, direct observational constraint to the low frequency PSD slope is hard due to limited lengths of the light curves available. Meanwhile, quasars show too large scatter in their DRW parameters to be attributed to the uncertainties in the measurements and the dependence of variation to known physical factors. In this work we present simulations showing that, if the low frequency PSD deviates from DRW, the red noise leakage can naturally produce large scatter in variation parameters measured from simulated light curves. The steeper the low frequency PSD slope is, the larger scatter we expect. Based on the observations of SDSS Stripe 82 quasars, we find the low frequency PSD slope should be no steeper than -1.3. The actual slope could be flatter, which consequently requires that quasar variabilities should be influenced by other unknown factors. We speculate that magnetic field and/or metallicity could be such additional factors., 27 pages, 10 figures, Accepted for publication in the Astrophysical Journal

We present Chandra ACIS-S X-ray imaging spectroscopy for five dual active galactic nucleus (AGN) candidates. Our targets were drawn from a sample of 1286 [O III]-selected AGN pairs systematically selected from the Sloan Digital Sky Survey Seventh Data Release. Each of the targets contains two nuclei separated by $\sim$3-9 kpc in projection, both of which are optically classified as Type 2 (obscured) AGNs based on diagnostic ratios of the narrow emission lines. Combined with independent, empirical star formation rate estimates based on the host-galaxy stellar continua, the new Chandra X-ray observations allow us to evaluate the dual-AGN hypothesis for each merging system. We confirm two (SDSS J0907+5203 and SDSS J1544+0446) of the five targets as bona-fide dual AGNs. For the other three targets, the existing data are consistent with the dual-AGN scenario, but we cannot rule out the possibility of stellar/shock heating and/or one AGN ionizing both gaseous components in the merger. The average X-ray-to-[O III] luminosity ratio in our targets seems to be systematically smaller than that observed in single AGNs but is higher than that seen in dual AGNs selected from AGNs with double-peaked narrow emission lines. We suggest that the systematically smaller X-ray-to-[O III] luminosity ratio observed in dual AGNs than in single AGNs is due to a high nuclear gas column likely from strong merger-induced inflows. Unlike double-peaked-[O III]-selected dual AGNs, the new sample selected from resolved galaxy pairs are not subject to the orientation bias caused by the double-peak line-of-sight velocity splitting selection, which also contributes to lowering the X-ray-to-[O III] luminosity ratio., Comment: 18 pages, 8 figures. Accepted for publication in ApJ. Comments welcome. arXiv admin note: text overlap with arXiv:1209.5418

We have collected near-infrared to X-ray data of 20 multi-epoch heavily reddened SDSS quasars to investigate the physical mechanism of reddening. Of these, J2317+0005 is found to be a UV cutoff quasar. Its continuum, which usually appears normal, decreases by a factor 3.5 at 3000{\AA}, compared to its more typical bright state during an interval of 23 days. During this sudden continuum cut-off, the broad emission line fluxes do not change, perhaps due to the large size of the Broad Line Region (BLR), r > 23 / (1+z) days. The UV continuum may have suffered a dramatic drop out. However, there are some difficulties with this explanation. Another possibility is that the intrinsic continuum did not change, but was temporarily blocked out, at least towards our line of sight. As indicated by X-ray observations, the continuum rapidly recovers after 42 days. A comparison of the bright state and dim states would imply an eclipse by a dusty cloud with a reddening curve having a remarkably sharp rise shortward of 3500{\AA}. Under the assumption of being eclipsed by a Keplerian dusty cloud, we characterized the cloud size with our observations, however, which is a little smaller than the 3000\AA\ continuum-emitting size inferred from accretion disk models. Therefore, we speculate this is due to a rapid outflow or inflow with a dusty cloud passing through our line-of-sight to the center., Comment: 11 pages, 6 figures, ApJ accepted

The optical/ultraviolet (UV) variability of quasars has been discovered to be correlated with other quasar properties, such as luminosity, black hole mass and rest-frame wavelength. However,the origin of variability has been a puzzle so far. In this work, we upgrade the accretion disc model (Li & Cao 2008), which assumed the variability is caused by the change of global mass accretion rate, by constraining the disc size to match the viscous timescale of accretion disc to the variability timescale observed and by including the irradiation/X-ray reprocessing to make the emitted spectrum become steeper. We find this hybrid model can reproduce the observed bluer-when-brighter trend quite well, which is used to validate the theoretical model by several works recently. The traditional correlation between the variability amplitude and rest-frame wavelength can also be well fitted by our model. In addition, a weak positive correlation between variability amplitude and black hole mass is present, qualitatively consistent with recent observations., 5 pages,5 figures,accepted by MNRAS Letters

We investigated the optical/ultraviolet (UV) color variations for a sample of 2169 quasars based on multi-epoch spectroscopy in the Sloan Digital Sky Survey (SDSS) data release seven (DR7) and data release nine (DR9). To correct the systematic difference between DR7 and DR9 due to the different instrumental setup, we produced a correction spectrum by using a sample of F-stars observed both in DR7 and DR9. The correction spectrum was then applied to quasars when comparing the spectra of DR7 with DR9. In each object, the color variation was explored by comparing the spectral index of the continuum power-law fit on the brightest spectrum with the faintest one, and also by the shape of their difference spectrum. In 1876 quasars with consistent color variations from two methods, we found that most sources (1755, $\sim 94\%$) show bluer-when-brighter (BWB) trend, and the redder-when-brighter (RWB) trend is only detected in 121 objects ($\sim 6\%$). The common BWB trend is supported by the bluer composite spectrum constructed from bright spectra than that from faint spectra, and also by the blue composite difference spectrum. The correction spectrum is proved to be highly reliable by comparing the composite spectrum from corrected DR9 and original DR7 spectra. Assuming that the optical/UV variability is triggered by fluctuations, RWB trend can likely be explained if the fluctuations occur firstly at outer disk region, and the inner disk region has not fully responded yet when the fluctuation being propagated inward. In contrast, the common BWB trend implies that the fluctuations are likely more often happening firstly in inner disk region., 24 pages, 9 figures, accepted for publication in ApJ

We present the results of extensive multi-band intra-night optical monitoring of BL Lacertae during 2010--2012. BL Lacertae was very active in this period and showed intense variability in almost all wavelengths. We extensively observed it for a total for 38 nights; on 26 of them observations were done quasi-simultaneously in B, V, R and I bands (totaling 113 light curves), with an average sampling interval of around 8 minutes. BL Lacertae showed significant variations on hour-like timescales in a total of 19 nights in different optical bands. We did not find any evidence for periodicities or characteristic variability time-scales in the light curves. The intranight variability amplitude is generally greater at higher frequencies and decreases as the source flux increases. We found spectral variations in BL Lacertae in the sense that the optical spectrum becomes flatter as the flux increases but in several flaring states deviates from the linear trend suggesting different jet components contributing to the emission at different times., Comment: 12 Pages, 5 figures, 3 Tables, Accepted for Publication in MNRAS

We compiled a sample of 60 quasars with spectroscopy on at least six epochs from the Sloan Digital Sky Survey (SDSS) to study the variabilities of the spectral shape, the continuum and the emission lines luminosity. In this paper, we present the results of SDSS J030639.57+000343.1. We found a strong anti-correlation between the continuum luminosity at $5100 \AA$ and the spectral index, implying a bluer-when-brighter trend. The luminosity of the broad $\rm H_\alpha$ line is proportion to the continuum luminosity at $5100 \AA$. Correspondingly, we did not find strong correlation between the equivalent width of broad $\rm H\alpha$ and the continuum luminosity, i.e. no baldwin effect of broad $\rm H\alpha$ in this source., Comment: 4 pages, 3 figures, accepted for publication in JAA

In a sample of 60 quasars selected from Sloan Digital Sky Survey (SDSS) with at least six-epoch spectroscopy, we investigate the variability of emission lines and continuum luminosity at various aspects. A strong anti-correlation between the variability and continuum luminosity at 2500 {\AA} is found for the sample, which is consistent with previous works. In individual sources, we find that half of the sample objects have a trend of being bluer-when-brighter trend (BWB), while the remaining half exhibit redder-when-brighter trend (RWB). Although the mechanism for RWB is unclear, the effects of host galaxy contribution due to seeing variations can not be completely ruled out. As expected from photoionization model, the positive correlations between the broad emission line and continuum luminosity are found in most individual sources, as well as for the whole sample. We confirm the Baldwin effect in most individual objects and the whole sample, while a negative Baldwin effect is also found in several quasars, which can be at least partly (if not all) due to the host galaxy contamination. We find positive correlations between the broad emission line luminosity and line width in most individual quasars, as well as the whole sample, implying a more variable line base than the line core., Comment: 30 pages,10 figures; Accepted, The Astrophysical Journal,2014