Sloan Digital Sky Survey Quasar Research Articles

Stellar Black Holes Can “Stretch” Supermassive Black Hole Accretion Disks

Stellar black holes (sBHs) are widely believed to exist in the accretion disks of active galactic nuclei (AGNs). Previous studies often focus on the transient emission produced by embedded sBHs. Here, we explore the possible observational consequences of an AGN accretion disk that contains a population of accreting sBHs. Embedded accreting sBHs change the effective temperature distribution of the AGN accretion disk by heating gas in the outer regions. Two possible observational consequences are presented. First, the spectral energy distribution has a turnover feature at ∼4700 Å when the supermassive black hole mass is ∼108 M ⊙, which can help explain the observed shallow spectral shape at wavelengths >5000 Å for the Sloan Digital Sky Survey quasar composite spectrum. Second, the half-light radius of a given relatively long wavelength is significantly larger than for an AGN disk without sBHs, which can be tested by microlensing observations. With appropriate sBH distributions, the model can be reconciled with quasar microlensing disk sizes. We propose that the half-light radius–wavelength relation can be utilized to investigate the distributions of embedded sBHs in AGN accretion disks.

Quasars with flare/eclipse-like variability identified in ZTF

ABSTRACT Active galactic nuclei (AGNs) are known to exhibit optical/ultraviolet variability and most of them can be well modelled by the damped random walks. Physical processes that are not related to the accretion disc, such as tidal disruption events (TDEs) or moving foreground dusty clouds, can cause flare-like and eclipse-like features in the optical light curve. Both long-term and high-cadence monitoring are needed to identify such features. By combining the Sloan Digital Sky Survey (SDSS), Panoramic Survey Telescope and Rapid Response System, and the Zwicky Transient Facility (ZTF) survey, we are able to identify a rare sample (11) out of the SDSS quasar catalogue ($\sim\!\! 83\, 000$). These quasars exhibit more or less constant brightness but show rapid optical variation in the epochs of ZTF Public Data Release 2. To investigate the possible origins of these flare/eclipse-like variabilities, we propose the second-epoch spectroscopic observations with the Gran Telescopio CANARIAS. We find that the change in accretion rate plays a significant role in these quasar variabilities. Among them, we identify two changing-look AGN candidates: SDSS J1427+2930 and SDSS J1420+3757. The luminosity change of the former may be caused by the enhanced supermassive black hole’s accretion or the TDE, while the latter is more related to the change in the accretion rate.

Absence of radio-bright dominance in a near-infrared selected sample of red quasars

Context. The dichotomy between red and blue quasars is still an open question. It is debated whether red quasars are simply blue quasars that are observed at certain inclination angles or if they provide insight into a transitional phase in the evolution of quasars. Aims. We investigate the relation between quasar colors and radio-detected fraction because radio observations of quasars provide a powerful tool in distinguishing between quasar models. Methods. We present the eHAQ+GAIA23 sample, which contains quasars from the High A(V) Quasar (HAQ) Survey, the Extended High A(V) Quasar (eHAQ) Survey, and the Gaia quasar survey. All quasars in this sample have been found using a near-infrared color selection of target candidates that have otherwise been missed by the Sloan Digital Sky Survey (SDSS). We implemented a redshift-dependent color cut in g* − i* to select red quasars in the sample and divided them into redshift bins, while using a nearest-neighbors algorithm to control for luminosity and redshift differences between our red quasar sample and a selected blue sample from the SDSS. Within each bin, we cross-matched the quasars to the Faint Images of the Radio Sky at Twenty centimeters (FIRST) survey and determined the radio-detection fraction. Results. For redshifts 0.8 < z ≤ 1.5, the red and blue quasars have a radio-detection fraction of 0.153−0.032+0.037 and 0.132−0.030+0.034, respectively. The red and blue quasars with redshifts 1.5 < z ≤ 2.4 have radio-detection fractions of 0.059−0.016+0.019 and 0.060−0.016+0.019, respectively, and the red and blue quasars with redshifts z > 2.4 have radio-detection fractions of 0.029−0.012+0.017 and 0.058−0.019+0.024, respectively. For the WISE color-selected red quasars, we find a radio-detection fraction of 0.160−0.034+0.038 for redshifts 0.8 < z ≤ 1.5, 0.063−0.017+0.020 for redshifts 1.5 < z ≤ 2.4, and 0.051−0.022+0.030 for redshifts z > 2.4. In other words, we find similar radio-detection fractions for red and blue quasars within < 1σ uncertainty, independent of redshift. This disagrees with what has been found in the literature for red quasars in SDSS. It should be noted that the fraction of broad absorption line (BAL) quasars in red SDSS quasars is about five times lower. BAL quasars have been observed to be more frequently radio quiet than other quasars, therefore the difference in BAL fractions could explain the difference in radio-detection fraction. Conclusions. The dusty torus of a quasar is transparent to radio emission. When we do not observe a difference between red and blue quasars, it leads us to argue that orientation is the main cause of quasar redness. Moreover, the observed higher proportion of BAL quasars in our dataset relative to the SDSS sample, along with the higher rate of radio detections, indicates an association of the redness of quasars and the inherent BAL fraction within the overall quasar population. This correlation suggests that the redness of quasars is intertwined with the inherent occurrence of BAL quasars within the entire population of quasars. In other words, the question why some quasars appear red or exhibit BAL characteristics might not be isolated; it could be directly related to the overall prevalence of BAL quasars in the quasar population. This finding highlights the need to explore the underlying factors contributing to both the redness and the frequency of BAL quasars, as they appear to be interconnected phenomena.

Comprehensive Connection among the Quasars with Different Types of Outflow Absorption Lines

It is commonly accepted that outflows from the central regions of quasars play a substantial role in regulating the global properties of the host galaxy. These outflows are typically detected through blueshifted absorption lines. However, the question remains whether outflows observed with different absorption line types indeed reflect the same environmental or evolutionary stage of the host galaxy. In this study, we use the Sloan Digital Sky Survey quasar catalog and employ the flux ratio of [O II] and [Ne V] emission lines as indicators to compare star formation rates (SFRs) within host galaxies of quasars exhibiting various outflow absorption line types: low-ionization broad absorption line (LoBAL), low-ionization Mini-BAL (LoMini-BAL), low-ionization narrow absorption line (LoNAL), high-ionization broad absorption line (HiBAL), high-ionization Mini-BAL (HiMini-BAL), and high-ionization narrow absorption line (HiNAL). Our findings indicate that the SFR of LoMini-BAL quasars is comparable to that of LoNAL quasars, somewhat less than that of LoBAL quasars, but markedly greater than that of HiBAL quasars. Furthermore, the SFR of HiMini-BAL quasars mirrors that of HiNAL or Non-abs (no associated absorption lines) quasars, but is significantly higher than that of HiBAL quasars. If we consider that differing absorption line types are indicative of the quasar evolution stage, our results propose an inclusive evolution sequence: LoBALs evolve into LoMini-BALs/LoNALs, then progress to HiBALs, and ultimately morph into HiMini-BALs/HiNALs/Non-abs. Concomitantly, the SFR within the host galaxies of quasars appears to decline noticeably nearing the LoNAL phase’s end and rejuvenates before the HiMini-BAL phase.

Multiwavelength Spectral Energy Distribution Analysis of X-Ray Selected Active Galactic Nuclei at z = 0.2–0.8 in the Stripe 82 Region

We perform a systematic, multiwavelength spectral energy distribution (SED) analysis of X-ray detected active galactic nuclei (AGNs) at z = 0.2–0.8 with Sloan Digital Sky Survey (SDSS) counterparts in the Stripe 82 region, consisting of 60 type 1 and 137 type 2 AGNs covering a 2–10 keV luminosity range of 41.6<logLx<44.7 . The latest CIGALE code, where dusty polar components are included, is employed. To obtain reliable host and AGN parameters in type 1 AGNs, we utilize the image-decomposed optical SEDs of host galaxies by Li et al. based on the Subaru Hyper-Suprime Cam images. The mean ratio of black hole masses (M BH) and stellar masses (M stellar) of our X-ray detected type 1 AGN sample, log(MBH/Mstellar)=−2.7±0.5 , is close to the local relation between BH and stellar masses as reported by Li et al. for SDSS quasars. This ratio is slightly lower than that found for more luminous ( logLbol>45 ) type 1 AGNs at z ∼ 1.5. This can be explained by the AGN luminosity dependence of log(MBH/Mstellar) , which little evolves with redshift. We confirm the trend that the UV-to-X-ray slope (α OX) or X-ray-to-bolometric correction factor (κ 2–10) increases with AGN luminosity or Eddington ratio. We find that type 1 and type 2 AGNs with the same luminosity ranges share similar host stellar mass distributions, whereas type 2s tend to show smaller AGN luminosities than type 1s. This supports the luminosity-dependent (or Eddington-ratio-dependent) unified scheme.

Estimators of Bolometric Luminosity and Black Hole Mass with Mid-infrared Continuum Luminosities for Dust-obscured Quasars: Prevalence of Dust-obscured SDSS Quasars

We present bolometric luminosity (L bol) and black hole (BH) mass (M BH) estimators based on mid-infrared (MIR) continuum luminosity (hereafter, L MIR) that are measured from infrared (IR) photometric data. The L MIR-based estimators are relatively immune from dust extinction effects, hence they can be used for dust-obscured quasars. To derive the L bol and M BH estimators, we use unobscured quasars selected from the Sloan Digital Sky Survey (SDSS) quasar catalog, which have wide ranges of L bol (1044.62–1046.16 erg s−1) and M BH (107.14–109.69 M ⊙). We find empirical relations between (i) continuum luminosity at 5100 Å (hereafter, L5100) and L MIR; (ii) L bol and L MIR. Using these relations, we derive the L MIR-based L bol and M BH estimators. We find that our estimators allow the determination of L bol and M BH at an accuracy of ∼0.2 dex against the fiducial estimates based on the optical properties of the unobscured quasars. We apply the L MIR-based estimators to SDSS quasars at z ≲ 0.5 including obscured ones. The ratios of L bol from the L MIR-based estimators to those from the optical luminosity-based estimators become larger with the amount of the dust extinction, and a non-negligible fraction (∼15%) of the SDSS quasars exhibits ratios greater than 1.5. This result suggests that dust extinction can significantly affect physical parameter derivations even for SDSS quasars, and that dust extinction needs to be carefully taken into account when deriving quasar properties.

QUOTAS: A New Research Platform for the Data-driven Discovery of Black Holes

We present QUOTAS, a novel research platform for the data-driven investigation of supermassive black hole (SMBH) populations. While SMBH data—observations and simulations—have grown in complexity and abundance, our computational environments and tools have not matured commensurately to exhaust opportunities for discovery. To explore the BH, host galaxy, and parent dark matter halo connection—in this pilot version—we assemble and colocate the high-redshift, z > 3 quasar population alongside simulated data at the same cosmic epochs. As a first demonstration of the utility of QUOTAS, we investigate correlations between observed Sloan Digital Sky Survey (SDSS) quasars and their hosts with those derived from simulations. Leveraging machine-learning algorithms (ML), to expand simulation volumes, we show that halo properties extracted from smaller dark-matter-only simulation boxes successfully replicate halo populations in larger boxes. Next, using the Illustris-TNG300 simulation that includes baryonic physics as the training set, we populate the larger LEGACY Expanse dark-matter-only box with quasars, and show that observed SDSS quasar occupation statistics are accurately replicated. First science results from QUOTAS comparing colocated observational and ML-trained simulated data at z3 are presented. QUOTAS demonstrates the power of ML, in analyzing and exploring large data sets, while also offering a unique opportunity to interrogate theoretical assumptions that underpin accretion and feedback models. QUOTAS and all related materials are publicly available at the Google Kaggle platform. (The full data set—observational data and simulation data—are available at: https://www.kaggle.com/ and the codes are available at:https://www.kaggle.com/datasets/quotasplatform/quotas)

Cosmological-scale Lyα Forest Absorption around Galaxies and AGNs Probed with the HETDEX and SDSS Spectroscopic Data

We present cosmological-scale three-dimensional neutral hydrogen (H i) tomographic maps at z = 2–3 over a total of 837 deg2 in two blank fields that are developed with Lyα forest absorptions of 14,736 background Sloan Digital Sky Survey (SDSS) quasars at z = 2.08–3.67. Using the tomographic maps, we investigate the large-scale (≳10 h −1 cMpc) average H i radial profiles and two-direction profiles of the line-of-sight (LOS) and transverse directions around galaxies and active galactic nuclei (AGNs) at z = 2–3 identified by the Hobby–Eberly Telescope Dark Energy eXperiment survey and SDSS, respectively. The peak of the H i radial profile around galaxies is lower than the one around AGNs, suggesting that the dark matter halos of galaxies are less massive on average than those of AGNs. The LOS profile of AGNs is narrower than the transverse profile, indicating the Kaiser effect. There exist weak absorption outskirts at ≳30 h −1 cMpc beyond H i structures of galaxies and AGNs found in the LOS profiles that can be explained by the H i gas at ≳30 h −1 cMpc falling toward the source position. Our findings indicate that the H i radial profile of AGNs has transitions from proximity zones (≲a few h −1 cMpc) to the H i structures (∼1–30 h −1 cMpc) and the weak absorption outskirts (≳30 h −1 cMpc). Although there is no significant dependence of AGN types (type 1 vs. type 2) on the H i profiles, the peaks of the radial profiles anticorrelate with AGN luminosities, suggesting that AGNs’ ionization effects are stronger than the gas mass differences.

Testing AGN outflow and accretion models with C iv and He ii emission line demographics in z ≈ 2 quasars

ABSTRACT Using ≈190 000 spectra from the 17th data release of the Sloan Digital Sky Survey (SDSS), we investigate the ultraviolet emission line properties in z ≈ 2 quasars. Specifically, we quantify how the shape of C iv λ1549 and the equivalent width (EW) of He ii λ1640 depend on the black hole mass and Eddington ratio inferred from Mg ii λ2800. Above L/LEdd ≳ 0.2, there is a strong mass dependence in both C iv blueshift and He ii EW. Large C iv blueshifts are observed only in regions with both high mass and high accretion rate. Including X-ray measurements for a subsample of 5000 objects, we interpret our observations in the context of AGN accretion and outflow mechanisms. The observed trends in He ii and 2 keV strength are broadly consistent with theoretical qsosed models of AGN spectral energy distributions (SEDs) for low spin black holes, where the ionizing SED depends on the accretion disc temperature and the strength of the soft excess. High spin models are not consistent with observations, suggesting SDSS quasars at z ≈ 2 may in general have low spins. We find a dramatic switch in behaviour at L/LEdd ≲ 0.1: the ultraviolet emission properties show much weaker trends, and no longer agree with qsosed predictions, hinting at changes in the structure of the broad line region. Overall, the observed emission line trends are generally consistent with predictions for radiation line driving where quasar outflows are governed by the SED, which itself results from the accretion flow and hence depends on both the SMBH mass and accretion rate.

ALMA High-Level Data Products: submillimetre counterparts of SDSS quasars in the ALMA footprint

ABSTRACT The Atacama Large Millimetre/submillimetre Array (ALMA) is the world’s most advanced radio interferometric facility, producing science data with an average rate of about 1 TB per day. After a process of calibration, imaging and quality assurance, the scientific data are stored in the ALMA Science Archive (ASA), along with the corresponding raw data, making the ASA an invaluable resource for original astronomical research. Due to their complexity, each ALMA data set has the potential for scientific results that go well beyond the ideas behind the original proposal that led to each observation. For this reason, the European ALMA Regional Centre initiated the High-Level Data Products initiative to develop science-oriented data products derived from data sets publicly available in the ASA, that go beyond the formal ALMA deliverables. The first instance of this initiative is the creation of a catalogue of submillimetre (submm) detections of Sloan Digital Sky Survey (SDSS) quasars from the SDSS Data Release 14 that lie in the aggregate ALMA footprint observed since ALMA Cycle 0. The ALMA fluxes are extracted in an automatic fashion, using the ALMA Data Mining Toolkit. All extractions above a signal-to-noise cut of 3.5 are considered, they have been visually inspected and the reliable detections are presented in a catalogue of 376 entries, corresponding to 275 unique quasars. Interesting targets found in the process, i.e. lensed or jetted quasars as well as quasars with nearby submm counterparts are highlighted, to facilitate further studies or potential follow up observations.

A deep 1.4 GHz survey of the J1030 equatorial field: A new window on radio source populations across cosmic time

We present deep L-Band observations of the equatorial field centered on the z = 6.3 Sloan Digital Sky Survey (SDSS) quasar (QSO). This field is rich of multiwavelength photometry and spectroscopy data, making it an ideal laboratory for galaxy evolution studies. Our observations reach a 1σ sensitivity of ~2.5 µJy at the center of the field. We extracted a catalog of 1489 radio sources down to a flux density of ~12.5 µJy (5σ) over a field of view of ~ 30′ diameter. We derived the source counts accounting for catalog reliability and completeness, and compared them with others available in the literature. Our source counts are among the deepest available so far, and, overall, are consistent with recent counts’ determinations and models. They show a slight excess at flux densities ~50 µJy, possibly associated with the presence of known overdensities in the field. We detected for the first time in the radio band the SDSS J1030+0524 QSO (26 ± 5 µJy, 8σ significance level). For this object, we derived an optical radio loudness RO = 0.62±0.12, which makes it the most radio quiet among active galactic nuclei (AGN) discovered so far at z ≳ 6 and detected at radio wavelengths. We unveiled extended diffuse radio emission associated with the lobes of a bright Fanaroff-Riley type II (FRII) radio galaxy located close to the center of the J1030 field, which is likely to become the future brightest cluster galaxy of a protocluster at z = 1.7. The lobes’ complex morphology, coupled with the presence of X-ray diffuse emission detected around the FRII galaxy lobes, may point toward an interaction between the radio jets and the external medium. We also investigated the relation between radio and X-ray luminosity for a sample of 243 X-ray-selected objects obtained from 500 ks Chandra observations of the same field, and spanning a wide redshift range (0 ≲ z ≲ 3). Focused on sources with a spectroscopic redshift and classification, we found that sources hosted by early-type galaxies and AGN follow log(LR)/log(LX) linear correlations with slopes of ~0.6 and ~0.8, respectively. This is interpreted as a likely signature of different efficiency in the accretion process. Finally, we found that most of these sources (≳87%) show a radio-to-X-ray radio loudness RX ≲ −3.5, classifying these objects as radio quiet.

The optical variability analysis of quasars based on fractal Brownian motion in Sloan digital sky survey

<p indent="0mm">The optical variability analysis of quasars based on fractal Brownian motion in the Sloan Digital Sky Survey (SDSS). The violent activities of the high energy variability can be observed in multiple wavelength bands of quasars, and studying quasar variability can contribute to a deeper understanding of the physical mechanism of the emission from the central region. Many studies have shown that the great majority of quasar variability is considered to originate in the accretion disk’s thermal emission, but some people do not agree with this view. The Hurst exponent of SDSS quasar variability is calculated using the rescaled range analysis method to investigate the fractal properties of variability in quasars from the SDSS stripe 82. The Hurst exponent results of variability show that the quasar variability is self-similar and long-range correlated. The variability obeys the biased random walk model, and its variation is similar to fractal Brownian motion. The power spectral exponent obtained from the Hurst exponent conversion shows that most of the variabilities can be approximately fitted by the damped random walk model, but a small number of samples have large deviations.

On the Kinematics of Cold, Metal-enriched Galactic Fountain Flows in Nearby Star-forming Galaxies

We use medium-resolution Keck/Echellette Spectrograph and Imager spectroscopy of bright quasars to study cool gas traced by Ca ii λλ3934, 3969 and Na i λλ5891, 5897 absorption in the interstellar/circumgalactic media of 21 foreground star-forming galaxies at redshifts 0.03 < z < 0.20 with stellar masses 7.4 ≤ log M */M ⊙ ≤ 10.6. The quasar–galaxy pairs were drawn from a unique sample of Sloan Digital Sky Survey quasar spectra with intervening nebular emission, and thus have exceptionally close impact parameters (R ⊥ < 13 kpc). The strength of this line emission implies that the galaxies’ star formation rates (SFRs) span a broad range, with several lying well above the star-forming sequence. We use Voigt profile modeling to derive column densities and component velocities for each absorber, finding that column densities N(Ca ii) > 1012.5 cm−2 (N(Na i) > 1012.0 cm−2) occur with an incidence f C(Ca ii) = 0.63+0.10 −0.11 (f C(Na i) = 0.57+0.10 −0.11). We find no evidence for a dependence of f C or the rest-frame equivalent widths W r (Ca ii K) or W r (Na i 5891) on R ⊥ or M *. Instead, W r (Ca ii K) is correlated with local SFR at >3σ significance, suggesting that Ca ii traces star formation-driven outflows. While most of the absorbers have velocities within ±50 km s−1 of the host redshift, their velocity widths (characterized by Δv 90) are universally 30–177 km s−1 larger than that implied by tilted-ring modeling of the velocities of interstellar material. These kinematics must trace galactic fountain flows and demonstrate that they persist at R ⊥ > 5 kpc. Finally, we assess the relationship between dust reddening and W r (Ca ii K) (W r (Na i 5891)), finding that 33% (24%) of the absorbers are inconsistent with the best-fit Milky Way E(B−V)-W r relations at >3σ significance.

The relation between quasars’ optical spectra and variability

Abstract Brightness variation is an essential feature of quasars, but its mechanism and relationship to other physical quantities are not understood well. We aimed to find the relationship between the optical variability and spectral features to reveal the regularity behind the random variation. It is known that a quasar’s Fe ii/Hβ flux ratio and equivalent width of [O iii]5007 are negatively correlated; this is called Eigenvector 1. In this work, we visualized the relationship between the position on this Eigenvector 1 (EV1) plane and how the brightness of the quasars had changed after ∼10 yr. We conducted three analyses, using a different quasar sample in each. The first analysis showed the relation between the quasars’ distributions on the EV1 plane and how much they had changed brightness, using 13438 Sloan Digital Sky Survey quasars. This result shows how brightness changes later are clearly related to the position on the EV1 plane. In the second analysis, we plotted the sources reported as “changing-look quasars” (or “changing-state quasars”) on the EV1 plane. This result shows that the position on the EV1 plane corresponds to the activity level of each source, and the bright or dim states of them are distributed on the opposite sides divided by the typical quasar distribution. In the third analysis, we examined the transition vectors on the EV1 plane using sources with multiple-epoch spectra. This result shows that the brightening and dimming sources move on a similar path and they reach a position corresponding to the opposite activity level. We also found this trend is opposite to the empirical rule that $R_{\rm {Fe\, \small {II}}}$ positively correlated with the Eddington ratio, which has been proposed based on the trends of a large number of quasars. From all these analyses, it is indicated that quasars tend to oscillate between both sides of the distribution ridge on the EV1 plane; each of them corresponds to a dim state and a bright state. This trend in optical variation suggests that significant brightness changes, such as changing-look quasars, are expected to repeat.

Concordance between Observations and Simulations in the Evolution of the Mass Relation between Supermassive Black Holes and Their Host Galaxies

We carry out a comparative analysis of the relation between the mass of supermassive black holes (BHs) and the stellar mass of their host galaxies at 0.2 < z < 1.7 using well-matched observations and multiple state-of-the-art simulations (e.g., MassiveBlackII, Horizon-AGN, Illustris, TNG, and a semianalytic model). The observed sample consists of 646 uniformly selected Sloan Digital Sky Survey quasars (0.2 < z < 0.8) and 32 broad-line active galactic nuclei (AGNs; 1.2 < z < 1.7) with imaging from Hyper Suprime-Cam (HSC) for the former and Hubble Space Telescope (HST) for the latter. We first add realistic observational uncertainties to the simulation data and then construct a simulated sample in the same manner as the observations. Over the full redshift range, our analysis demonstrates that all simulations predict a level of intrinsic scatter of the scaling relations comparable to the observations that appear to agree with the dispersion of the local relation. Regarding the mean relation, Horizon-AGN and TNG are in closest agreement with the observations at low and high redshift (z ∼ 0.2 and 1.5, respectively), while the other simulations show subtle differences within the uncertainties. For insight into the physics involved, the scatter of the scaling relation, seen in the SAM, is reduced by a factor of two and closer to the observations after adopting a new feedback model that considers the geometry of the AGN outflow. The consistency in the dispersion with redshift in our analysis supports the importance of both quasar- and radio-mode feedback prescriptions in the simulations. Finally, we highlight the importance of increasing the sensitivity (e.g., using the James Webb Space Telescope), thereby pushing to lower masses and minimizing biases due to selection effects.

Inferences on Relations between Distant Supermassive Black Holes and Their Hosts Complemented by the Galaxy Fundamental Plane

The realization of fundamental relations between supermassive black holes and their host galaxies would have profound implications in astrophysics. To add further context to studies of their coevolution, an investigation is carried out to gain insight as to whether quasars and their hosts at earlier epochs follow the local relation between black hole mass (M BH) and stellar velocity dispersion (σ *). We use 584 Sloan Digital Sky Survey quasars at 0.2 < z < 0.8 with black hole measurements and properties of their hosts from the Hyper Suprime-Cam Subaru Strategic Program. An inference of σ * is achieved for each based on the total stellar mass (M *) and size of the host galaxy by using the galaxy mass fundamental plane for inactive galaxies at similar redshifts. In agreement with past studies, quasars occupy elevated positions from the local M BH−σ * relation which can be considered as a flattening of the relation. Based on a simulated sample, we demonstrate that an evolving intrinsic M BH−σ * relation can match the observations. However, we hypothesize that these changes are simply a consequence of a nonevolving intrinsic relationship between M BH and M *. Reassuringly, there is evidence of migration onto the local M BH−σ * for galaxies that are either massive, quiescent or compact. Thus, the bulges of quasar hosts at high redshift are growing and likely to align onto the mass scaling relation with their black holes at later times.

Spectroscopic redshift determination with Bayesian convolutional networks

Astronomy is facing large amounts of data, so astronomers have to rely on automated methods to analyse them. However, automated methods might produce incorrect values. Therefore, we need to develop different automated methods and perform a consistency check to identify them. If there is a lot of labelled data, convolutional neural networks are a powerful method for any task. We illustrate the consistency check on spectroscopic redshift determination with a method based on a Bayesian convolutional neural network inspired by VGG networks. The method provides predictive uncertainties that enable us to (1.) determine unusual or problematic spectra for visual inspection; (2.) do thresholding that allows us to balance between the error of redshift predictions and coverage. We used the 12th Sloan Digital Sky Survey quasar superset as the training set for the method. We evaluated its generalisation capability on about three-quarters of a million spectra from the 16th quasar superset of the same survey. On the 16th quasar superset, the method performs better in terms of the root-mean-squared error than the most used template fitting method. Using redshift predictions of the proposed method, we identified spectra with incorrectly determined redshifts that are unrecognised quasars or were misclassified as them.

Discovery of a Cool, Metal-rich Gas Reservoir in the Outskirts of z ≈ 0.5 Clusters

We built the first-ever statistically significant sample of ≈80,000 background quasar–foreground cluster pairs to study the cool, metal-rich gas in the outskirts (>R 500) of z ≈ 0.5 clusters with a median mass of ≈1014.2 M ⊙. The sample was obtained by crossmatching the Sloan Digital Sky Survey (SDSS) cluster catalog of Wen & Han and the SDSS quasar catalog of Lyke et al. The median impact parameter (ρ cl) of the clusters from the quasar sightlines is 2.4 Mpc (median ρ cl/R 500 = 3.6). A strong Mg II, along with marginal Fe II, absorption is detected in the mean and median stacked spectra of the quasars with a total Mg II rest-frame equivalent width () of 0.034 ± 0.005 Å (7σ) and 0.010 ± 0.003 Å (3σ), respectively. The shows a declining trend with increasing ρ cl and ρ cl/R 500, but does not show any significant trend with mass (M 500) or redshift (z cl) within the small M 500 and z cl ranges probed here. The Mg II absorption signal and the trends persist even if we exclude the quasar–cluster pairs where the background quasars may be probing the circumgalactic medium of bright galaxies with impact parameters <300 kpc. The Mg II (and Fe II) absorption reported here is the first detection of its kind. It indicates the presence of a cool, metal-rich gas reservoir surrounding galaxy clusters out to several R 500. We suggest that the metal-rich gas in the cluster outskirts arise from stripped materials and that gas stripping may be important out to large clustocentric distances (>3R 500).

MALS SALT-NOT Survey of MIR-selected Powerful Radio-bright AGN at 0 < z < 3.5

We present results of an optical spectroscopic survey using SALT and the Nordic Optical Telescope to build a Wide-field Infrared Survey Explorer mid-infrared color-based, dust-unbiased sample of powerful radio-bright (>200 mJy at 1.4 GHz) active galactic nuclei (AGN) for the MeerKAT Absorption Line Survey (MALS). Our sample has 250 AGN (median z = 1.8) showing emission lines, 26 with no emission lines, and 27 without optical counterparts. Overall, our sample is fainter (Δi = 0.6 mag) and redder (Δ(g−i) = 0.2 mag) than radio-selected quasars, and representative of fainter quasar population detected in optical surveys. About 20% of the sources are narrow-line AGN (NLAGN)–65% of these, at z < 0.5 are galaxies without strong nuclear emission, and 10% at z > 1.9, have emission line ratios similar to radio galaxies. The farthest NLAGN in our sample is M1513-2524 (z em = 3.132), and the largest radio source (size ∼330 kpc) is M0909-3133 (z em = 0.884). We discuss in detail 110 AGN at 1.9 < z < 3.5. Despite representing the radio loudest quasars (median R = 3685), their Eddington ratios are similar to the Sloan Digital Sky Survey quasars having lower R. We detect four C iv broad-absorption line (BAL) QSOs, all among AGN with least R, and highest black hole masses and Eddington ratios. The BAL detection rate (%) is consistent with that seen in extremely powerful (L 1.4GHz > 1025 W Hz−1) quasars. Using optical light curves, radio polarization, and γ-ray detections, we identify seven high-probability BL Lacertae objects. We also summarize the full MALS footprint to search for H i 21 cm and OH 18 cm lines at z < 2.

Fitting AGN/Galaxy X-Ray-to-radio SEDs with CIGALE and Improvement of the Code

Abstract Modern and future surveys effectively provide a panchromatic view for large numbers of extragalactic objects. Consistently modeling these multiwavelength survey data is a critical but challenging task for extragalactic studies. The Code Investigating GALaxy Emission (cigale) is an efficient python code for spectral energy distribution (SED) fitting of galaxies and active galactic nuclei (AGNs). Recently, a major extension of cigale (named x-cigale) has been developed to account for AGN/galaxy X-ray emission and improve AGN modeling at UV-to-IR wavelengths. Here, we apply x-cigale to different samples, including Cosmological Evolution Survey (COSMOS) spectroscopic type 2 AGNs, Chandra Deep Field-South X-ray detected normal galaxies, Sloan Digital Sky Survey quasars, and COSMOS radio objects. From these tests, we identify several weaknesses of x-cigale and improve the code accordingly. These improvements are mainly related to AGN intrinsic X-ray anisotropy, X-ray binary emission, AGN accretion-disk SED shape, and AGN radio emission. These updates improve the fit quality and allow for new interpretation of the results, based on which we discuss physical implications. For example, we find that AGN intrinsic X-ray anisotropy is moderate, and can be modeled as L X ( θ ) ∝ 1 + cos θ , where θ is the viewing angle measured from the AGN axis. We merge the new code into the major branch of cigale, and publicly release this new version as cigale v2022.0 on https://cigale.lam.fr.