Search for: All records

Abstract Dust-obscured galaxies (DOGs) with extremely red optical-to-infrared colors are often associated with intense starburst and active galactic nucleus (AGN) activity. Studying DOGs can provide insights into the processes that drive the growth of galaxies and their central supermassive black holes. However, the general DOG population is heterogeneous, spanning a wide range of evolutionary stages, and has X-ray obscuring column densities (N_H) covering low to high levels. In this work, we focus on seven high Eddington ratio DOGs ( $\log {\lambda }_{\mathrm{Edd}}\gtrsim -0.5$) to examine their X-ray obscuration properties using new and archival X-ray observations. We confirm that these systems are generally heavily obscured, with six out of seven havingN_H ≳ 10²³cm⁻²and three out of seven havingN_H ≳ 10²⁴cm⁻². Based on the observed similarity with the rare Hot DOG population, we argue that both high-λ_EddDOGs and Hot DOGs likely trace the postmerger phase, during which AGNs are enshrouded by large columns of dust-rich material. 

Abstract The coevolution of supermassive black holes and their host galaxies represents a fundamental question in astrophysics. One approach to investigating this question involves comparing the star formation rates (SFRs) of active galactic nuclei (AGNs) with those of typical star-forming galaxies. At relatively low redshifts (z≲ 1), radio AGNs manifest diminished SFRs, indicating suppressed star formation, but their behavior at higher redshifts is unclear. To examine this, we leveraged galaxy and radio-AGN data from the well-characterized W-CDF-S, ELAIS-S1, and XMM-LSS fields. We established two mass-complete reference star-forming galaxy samples and two radio-AGN samples, consisting of 1763 and 6766 radio AGNs, the former being higher in purity and the latter more complete. We subsequently computed star-forming fractions (f_SF; the fraction of star-forming galaxies to all galaxies) for galaxies and radio-AGN host galaxies and conducted a robust comparison between them up toz≈ 3. We found that the tendency for radio AGNs to reside in massive galaxies primarily accounts for their lowf_SF, which also shows a strong negative dependence uponM_⋆and a strong positive evolution withz. To investigate further the star formation characteristics of those star-forming radio AGNs, we constructed the star-forming main sequence (MS) and investigated the behavior of the position of AGNs relative to the MS atz≈ 0–3. Our results reveal that radio AGNs display lower SFRs than star-forming galaxies in the low-zand high-M_⋆regime and, conversely, exhibit comparable or higher SFRs than MS star-forming galaxies at higher redshifts or lowerM_⋆. 

Abstract Dust-obscured galaxies (DOGs) are enshrouded by dust and many are believed to host accreting supermassive black holes (SMBHs), which makes them unique objects for probing the coevolution of galaxies and SMBHs. We select and characterize DOGs in the 13 deg²XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS), leveraging the superb multiwavelength data—from X-rays to radio. We select 3738 DOGs atz≈ 1.6–2.1 in XMM-SERVS, while maintaining good data quality without introducing significant bias. This represents the largest DOG sample with thorough multiwavelength source characterization. Spectral energy distribution modeling shows DOGs are a heterogeneous population consisting of both normal galaxies and active galactic nuclei (AGNs). Our DOGs are massive ( $\log M_{\star }/M_{☉}\approx 10.7-11.3$), 174 are detected in X-rays, and they are generally radio-quiet systems. X-ray detected DOGs are luminous and are moderately to heavily obscured in X-rays. Stacking analyses for the X-ray undetected DOGs show highly significant average detections. Critically, we compare DOGs with matched galaxy populations. DOGs have similar AGN fractions compared with typical galaxy populations. X-ray detected DOGs have higherM_⋆and higher X-ray obscuration, but they are not more star-forming than typical X-ray AGNs. Our results potentially challenge the relevance of the merger-driven galaxy-SMBH coevolution framework for X-ray detected DOGs. 

Abstract Supermassive black holes (SMBHs) can grow through both accretion and mergers. It is still unclear how SMBHs evolve under these two channels from high redshifts to the SMBH population we observe in the local Universe. Observations can directly constrain the accretion channel but cannot effectively constrain mergers yet, while cosmological simulations provide galaxy merger information but can hardly return accretion properties consistent with observations. In this work, we combine the observed accretion channel and the simulated merger channel, taking advantage of observations and cosmological simulations, to depict a realistic evolution pattern of the SMBH population. With this methodology, we can derive the scaling relation between the black hole mass (M_BH) and host-galaxy stellar mass (M_⋆), and the local black hole mass function (BHMF). Our scaling relation is lower than those based on dynamically measuredM_BH, supporting the claim that dynamically measured SMBH samples may be biased. We show that the scaling relation has little redshift evolution. The BHMF steadily increases fromz= 4 toz= 1 and remains largely unchanged fromz= 1 toz= 0. The overall SMBH growth is generally dominated by the accretion channel, with possible exceptions at high mass (M_BH≳ 10⁸M_⊙orM_⋆≳ 10¹¹M_⊙) and low redshift (z≲ 1). We also predict that around 25% of the total SMBH mass budget in the local Universe may be locked within long-lived, wandering SMBHs, and the wandering mass fraction and wandering SMBH counts increase withM_⋆. 

Abstract Dust-obscured galaxies (DOGs) containing central supermassive black holes (SMBHs) that are rapidly accreting (i.e., having high Eddington ratios,λ_Edd) may represent a key phase closest to the peak of both the black hole and galaxy growth in the coevolution framework for SMBHs and galaxies. In this work, we present a 68 ks XMM-Newton observation of the high-λ_EddDOG J1324+4501 atz∼ 0.8, which was initially observed by Chandra. We analyze the XMM-Newton spectra jointly with archival Chandra spectra. In performing a detailed X-ray spectral analysis, we find that the source is intrinsically X-ray luminous with $\log (L_{\mathrm{X}}$/erg ${\mathrm{s}}^{-1})={44.71}_{-0.12}^{+0.08}$and heavily obscured with $\log \left(N_{\mathrm{H}}/{\mathrm{cm}}^{-2}\right)={23.43}_{-0.13}^{+0.09}$. We further utilize UV-to-IR archival photometry to measure and fit the source’s spectral energy distribution to estimate its host-galaxy properties. We present a supplementary comparison sample of 21 X-ray luminous DOGs from the XMM-SERVS survey with sufficient (>200) 0.5–10 keV counts to perform a similarly detailed X-ray spectral analysis. Of the X-ray luminous DOGs in our sample, we find that J1324+4501 is the most remarkable, possessing one of the highest X-ray luminosities, column densities, and star formation rates. We demonstrate that J1324+4501 is in an extreme evolutionary stage where SMBH accretion and galaxy growth are at their peaks. 

Abstract We present the rest-frame ultraviolet−optical spectral properties of 65 broad absorption line (BAL) quasars from the Gemini Near Infrared Spectrograph−Distant Quasar Survey (GNIRS-DQS). These properties are compared with those of 195 non-BAL quasars from GNIRS-DQS in order to identify the drivers for the appearance of BALs in quasar spectra. In particular, we compare equivalent widths and velocity widths, as well as velocity offsets from systemic redshifts, of principal emission lines. In spite of the differences between their rest-frame ultraviolet spectra, we find that luminous BAL quasars are generally indistinguishable from their non-BAL counterparts in the rest-frame optical band at redshifts 1.55 ≲z≲ 3.50. We do not find any correlation between BAL trough properties and the Hβ-based supermassive black hole masses and normalized accretion rates in our sample. Considering the Sloan Digital Sky Survey quasar sample, which includes the GNIRS-DQS sample, we find that a monochromatic luminosity at rest-frame 2500 Å of ≳10⁴⁵erg s⁻¹is a necessary condition for launching BAL outflows in quasars. We compare our findings with other BAL quasar samples and discuss the roles that accretion rate and orientation play in the appearance of BAL troughs in quasar spectra. 

Abstract We present an observational study of wind acceleration based on four low-ionization broad absorption line (BAL) quasars (J0136, J1238, J1259, and J1344). J0136 and J1344 (group 1) are radio-quiet and show large BAL-velocity shifts as opposed to stable line-locking associated absorption lines (AALs). Notably, J1344 displays a linear relation between BAL-velocity shift and time interval over three consecutive epochs, characteristic of compelling evidence for BAL acceleration. J1238 and J1259 (group 2) exhibit small BAL-velocity shifts along with steep-spectrum, weak radio emission at 3.0 and 1.4 GHz. All four quasars have spectral energy distributions (SEDs) with a peak atλ_rest∼ 10μm, suggesting a link between the BAL acceleration and hot dust emission. The group-2 quasars are redder than group-1 quasars and have a steeper rise at 1μm <λ_rest< 3μm in their SEDs. All but J1238 exhibit a steep rise followed by a plateau-like time evolution in BAL-velocity shift. Our investigations, combined with previous studies of BAL acceleration, indicate that (1) the coupling process between the BALs and the interstellar medium (ISM) is one of the major avenues for the origin of quasar reddening and patchy obscuration, (2) AAL outflows are ubiquitous and likely signify large-scale remnants of BAL winds coupled to the ISM, and (3) wind deceleration that is closely linked to the BAL–ISM coupling process may produce weak radio emission in otherwise radio-quiet quasars. 

ABSTRACT Quasar winds can shock and sweep up ambient interstellar medium (ISM) gas, contributing to galactic quenching. We combine and extend past models of energy-conserving shock bubbles around quasars, investigate model implications from an observational standpoint, and test model predictions using new high-resolution spectroscopic observations of the broad absorption-line quasar SDSS J030000.56+004828.0 (J0300). Even with constant energy input from the wind, a bubble’s expansion decelerates over time as more ISM gas is swept up. Our new observations enable a direct search for this deceleration. We obtain the tightest reported 3σ limit on the average rest-frame deceleration (or acceleration) of a quasar outflow: |a| < 0.1 km s−1 yr−1 (<3 × 10−4 cm s−2) in the relatively low-velocity Ca ii outflow of J0300 over 9.65 rest-frame years. We can satisfy these limits with certain parameter choices in our model, but the large velocity range of the Ca ii absorption in J0300 rules out the hypothesis that such gas shares the velocity of the swept-up ISM gas in a self-similar shock bubble. We investigate the possibility of ram-pressure acceleration of preexisting ISM clouds and conclude that the velocity range seen in Ca ii in J0300 is potentially consistent with such an explanation. The Ca ii-absorbing gas clouds in J0300 have been inferred to have high densities by Choi et al., in which case they can only have been accelerated to their current speeds if they were originally at least an order of magnitude less dense than they are today. 

Abstract Periodic signatures in time-domain observations of quasars have been used to search for binary supermassive black holes (SMBHs). These searches, across existing time-domain surveys, have produced several hundred candidates. The general stochastic variability of quasars, however, can masquerade as a false-positive periodic signal, especially when monitoring cadence and duration are limited. In this work, we predict the detectability of binary SMBHs in the upcoming Rubin Observatory Legacy Survey of Space and Time (LSST). We apply computationally inexpensive sinusoidal curve fits to millions of simulated LSST Deep Drilling Field light curves of both single, isolated quasars and binary quasars. The period and phase of simulated binary signals can generally be disentangled from quasar variability. Binary amplitude is overestimated and poorly recovered for two-thirds of potential binaries due to quasar accretion variability. Quasars with strong intrinsic variability can obscure a binary signal too much for recovery. We also find that the most luminous quasars mimic current binary candidate light curves and their properties: The false-positive rates are 60% for these quasars. The reliable recovery of binary period and phase for a wide range of input binary LSST light curves is promising for multi-messenger characterization of binary SMBHs. However, pure electromagnetic detections of binaries using photometric periodicity with amplitude greater than 0.1 mag will result in samples that are overwhelmed by false positives. This paper represents an important and computationally inexpensive way forward for understanding the true and false-positive rates for binary candidates identified by Rubin. 

Abstract We present the results of an investigation of a highly variable Civbroad absorption line (BAL) feature in spectra of the quasar SBS 1408+544 (z= 2.337) that shows a significant shift in velocity over time. This source was observed as a part of the Sloan Digital Sky Survey (SDSS) Reverberation Mapping project and the SDSS-V Black Hole Mapper Reverberation Mapping project, and has been included in two previous studies, both of which identified significant variability in a high-velocity CivBAL on timescales of just a few days in the quasar rest frame. Using ∼130 spectra acquired over 8 yr of spectroscopic monitoring with SDSS, we have determined that this BAL is not only varying in strength, but is also systematically shifting to higher velocities. Using cross-correlation methods, we measure the velocity shifts (and corresponding acceleration) of the BAL over a wide range of timescales, measuring an overall velocity shift of $\mathrm{\Delta }v=-{683}_{-84}^{+89}$km s⁻¹over the 8 yr monitoring period. This corresponds to an average rest-frame acceleration ofa= 1.04 ${}_{-0.13}^{+0.14}$cm s⁻², though the magnitude of the acceleration on shorter timescales is not constant throughout. We place our measurements in the context of BAL-acceleration models and examine various possible causes of the observed velocity shift.