"Changing-look" active galactic nuclei (CL-AGNs) challenge our basic ideas about the physics of accretion flows and circumnuclear gas around supermassive black holes. Using first-year Sloan Digital Sky Survey V (SDSS-V) repeated spectroscopy of nearly 29,000 previously known AGNs, combined with dedicated follow-up spectroscopy, and publicly available optical light curves, we have identified 116 CL-AGNs where (at least) one broad emission line has essentially (dis-)appeared, as well as 88 other extremely variable systems. Our CL-AGN sample, with 107 newly identified cases, is the largest reported to date, and includes $\sim0.4\%$ of the AGNs reobserved in first-year SDSS-V operations. Among our CL-AGNs, 67% exhibit dimming while 33% exhibit brightening. Our sample probes extreme AGN spectral variability on months to decades timescales, including some cases of recurring transitions on surprisingly short timescales ($\lesssim 2$ months in the rest frame). We find that CL events are preferentially found in lower-Eddington-ratio ($f_{Edd}$) systems: Our CL-AGNs have a $f_{Edd}$ distribution that significantly differs from that of a carefully constructed, redshift- and luminosity-matched control sample (Anderson-Darling test yielding $p_{\rm AD}\approx 6\times10^{-5}$; median $f_{Edd}\approx0.025$ vs. $0.043$). This preference for low $f_{Edd}$ strengthens previous findings of higher CL-AGN incidence at lower $f_{Edd}$, found in smaller samples. Finally, we show that the broad MgII emission line in our CL-AGN sample tends to vary significantly less than the broad H$\beta$ emission line. Our large CL-AGN sample demonstrates the advantages and challenges in using multi-epoch spectroscopy from large surveys to study extreme AGN variability and physics.
ABSTRACT Highly radio-loud quasars (HRLQs; log R > 2.5) at z ≳ 4 show apparent enhanced X-ray emission compared to matched HRLQs at lower redshifts, perhaps due to a redshift-dependent fractional contribution to the X-ray luminosity from inverse-Compton scattering of cosmic microwave background photons (IC/CMB). Using new Chandra observations and archival X-ray data, we investigate this phenomenon with an optically flux-limited sample of 41 HRLQs at z = 4–5.5 all with sensitive X-ray coverage, the largest sample utilized to date by a wide margin. X-ray enhancements are assessed using X-ray-to-optical flux ratios and spectral energy distributions. We confirm the presence of X-ray enhancements at a 4.9–5.3σ significance level, finding that the median factor of enhancement is ≈1.8 at our sample median redshift of z ≈ 4.4. Under a fractional IC/CMB model, the expected enhancement at lower redshifts is modest; e.g. ≈4 per cent at z ≈ 1.5. We also investigate a sample of seven radio-loud quasars (RLQs; log R > 1) at even higher redshifts of z = 5.6–6.8, using new and archival X-ray data. These RLQs also show evidence for X-ray enhancements by a median factor of ≈2.7 at a 3.7–4.9σ significance level. The X-ray spectral and other properties of these z = 5.6–6.8 RLQs, however, pose challenges for a straightforward fractional IC/CMB interpretation of their enhancements.
High-quality Extragalactic Legacy-field Monitoring (HELM) is a long-term observing program that photometrically monitors several well-studied extragalactic legacy fields with the Dark Energy Camera (DECam) imager on the CTIO 4 m Blanco telescope. Since 2019 February, HELM has been monitoring regions within COSMOS, XMM-LSS, CDF-S, S-CVZ, ELAIS-S1, and SDSS Stripe 82 with few-day cadences in the (u)gri(z) bands, over a collective sky area of similar to 38 deg2. The main science goal of HELM is to provide high-quality optical light curves for a large sample of active galactic nuclei (AGNs), and to build decades-long time baselines when combining past and future optical light curves in these legacy fields. These optical images and light curves will facilitate the measurements of AGN reverberation mapping lags, as well as studies of AGN variability and its dependencies on accretion properties. In addition, the time-resolved and coadded DECam photometry will enable a broad range of science applications from galaxy evolution to time-domain science. We describe the design and implementation of the program and present the first data release that includes source catalogs and the first similar to 3.5 yr of light curves during 2019A-2022A.
High-quality Extragalactic Legacy-field Monitoring (HELM) is a long-term observing program that photometrically monitors several well-studied extragalactic legacy fields with the Dark Energy Camera (DECam) imager on the CTIO 4m Blanco telescope. Since Feb 2019, HELM has been monitoring regions within COSMOS, XMM-LSS, CDF-S, S-CVZ, ELAIS-S1, and SDSS Stripe 82 with few-day cadences in the (u)gri(z) bands, over a collective sky area of ∼ 38 deg^2. The main science goal of HELM is to provide high-quality optical light curves for a large sample of active galactic nuclei (AGNs), and to build decades-long time baselines when combining past and future optical light curves in these legacy fields. These optical images and light curves will facilitate the measurements of AGN reverberation mapping lags, as well as studies of AGN variability and its dependences on accretion properties. In addition, the time-resolved and coadded DECam photometry will enable a broad range of science applications from galaxy evolution to time-domain science. We describe the design and implementation of the program and present the first data release that includes source catalogs and the first ∼ 3.5 years of light curves during 2019A–2022A.
We present a routinized and reliable method to obtain source catalogs from the $\textit{Nuclear Spectroscopic Telescope Array}$ ($\textit{NuSTAR}$) extragalactic surveys of the Extended $\textit{Chandra}$ Deep Field-South (E-CDF-S) and $\textit{Chandra}$ Deep Field-North (CDF-N). The $\textit{NuSTAR}$ E-CDF-S survey covers a sky area of $\approx30'\times30'$ to a maximum depth of $\sim$ 230 ks corrected for vignetting in the 3--24 keV band, with a total of 58 sources detected in our E-CDF-S catalog; the $\textit{NuSTAR}$ CDF-N survey covers a sky area of $\approx7'\times10'$ to a maximum depth of $\sim$ 440 ks corrected for vignetting in the 3--24 keV band, with a total of 42 sources detected in our CDF-N catalog that is produced for the first time. We verify the reliability of our two catalogs by crossmatching them with the relevant catalogs from the $\textit{Chandra}$ X-ray observatory, and find that the fluxes of our $\textit{NuSTAR}$ sources are generally consistent with that of their $\textit{Chandra}$ counterparts. Our two catalogs are produced following the exactly same method and made publicly available, thereby providing a uniform platform that facilitates further studies involving these two fields. Our source-detection method provides a systematic approach for source cataloging in other $\textit{NuSTAR}$ extragalactic surveys.
We study the ensemble X-ray variability properties of Active Galactic Nuclei (AGN) over a large range of timescales (20 ks $\leq T\leq$ 14 yrs), redshift ($0\leq z \lesssim 3$), luminosities ($10^{40}$ erg s$^{-1}\leq L_X\leq 10^{46}$ erg s$^{-1}$) and black hole (BH) masses ($10^6 \leq $M$_\odot \leq 10^9$). We propose the use of the variance-frequency diagram, as a viable alternative to the study of the power spectral density (PSD), which is not yet accessible for distant, faint and/or sparsely sampled AGN. We show that the data collected from archival observations and previous literature studies are fully consistent with a universal PSD form which does not show any evidence for systematic evolution of shape or amplitude with redshift or luminosity, even if there may be differences between individual AGN at a given redshift or luminosity. We find new evidence that the PSD bend frequency depends on BH mass and, possibly, on accretion rate. We finally discuss the implications for current and future AGN population and cosmological studies.
The Hot Universe Baryon Surveyor (HUBS) is a proposed space-based X-ray telescope for detecting X-ray emissions from the hot gas content in our universe. With its unprecedented spatially-resolved high-resolution spectroscopy and large field of view, the HUBS mission will be uniquely qualified to measure the physical and chemical properties of the hot gas in the interstellar medium, the circumgalactic medium, the intergalactic medium, and the intracluster medium. These measurements will be valuable for two key scientific goals of HUBS, namely to unravel the AGN and stellar feedback physics that governs the formation and evolution of galaxies, and to probe the baryon budget and multi-phase states from galactic to cosmological scales. In addition to these two goals, the HUBS mission will also help us solve some problems in the fields of galaxy clusters, AGNs, diffuse X-ray backgrounds, supernova remnants, and compact objects. This paper discusses the perspective of advancing these fields using the HUBS telescope.
We present 850um imaging of the XMM-LSS field observed for 170 hours as part of the James Clerk Maxwell Telescope SCUBA-2 Large eXtragalactic Survey (S2LXS). S2LXS XMM-LSS maps an area of 9 square degrees, reaching a moderate depth of 1-sigma 4 mJy/beam. This is the largest contiguous area of extragalactic sky mapped by JCMT at 850um to date. The wide area of the S2LXS XMM-LSS survey allows us to probe the ultra-bright (S_850um > 15 mJy), yet rare submillimetre population. We present the S2LXS XMM-LSS catalogue, which comprises 40 sources detected at >5-sigma significance, with deboosted flux densities in the range of 7 mJy to 48 mJy. We robustly measure the bright-end of the 850um number counts at flux densities >7 mJy, reducing the Poisson errors compared to existing measurements. The S2LXS XMM-LSS observed number counts show the characteristic upturn at bright fluxes, expected to be motivated by local sources of submillimetre emission and high-redshift strongly lensed galaxies. We find that the observed 850um number counts are best reproduced by model predictions that include either strong lensing or source blending from a 15 arcsec beam, indicating that both may make an important contribution to the observed over-abundance of bright single dish 850um selected sources. We make the S2LXS XMM-LSS 850um map and >5-sigma catalogue presented here publicly available.
We present a study of the relation between the [OIII] 5007A emission line width (sigma_{[OIII]}) and stellar velocity dispersion (sigma_{*}), utilizing a sample of 740 type 1 active galactic nuclei (AGNs) with high-quality spectra at redshift z < 1.0. We find the broad correlation between the core component of [OIII] emission line width (sigma_{[OIII,core]}) and sigma_{*} with a scatter of 0.11~dex for the low redshift (z < 0.1) sample; for redshift (0.3 < z < 1.0) AGNs, the scatter is larger, being 0.16~dex. We also find that the Eddington ratio (L_{bol}/L_{Edd}) may play an important role in the discrepancies between sigma_{[OIII,core]} and sigma_{*}. As the L_{bol}/L_{Edd} increases, sigma_{[OIII,core]} tends to be larger than sigma_{*}. By classifying our local sample with different minor-to-major axis ratios, we find that sigma_{*} is larger than sigma_{[OIII,core]} for those edge-on spiral galaxies. In addition, we also find that the effects of outflow strength properties such as maximum outflow velocity (V_{max}) and the broader component of [OIII] emission line width and line shift (sigma_{[OIII,out]} and V_{[OIII,out]}) may play a major role in the discrepancies between sigma_{[OIII,core]} and sigma_{*}. The discrepancies between sigma_{[OIII,core]} and sigma_{*} are larger when V_{max}, V_{[OIII,out]}, and sigma_{[OIII,out]} increase. Our results show that the outflow strengths may have significant effects on the differences between narrow-line region gas and stellar kinematics in AGNs. We suggest that caution should be taken when using sigma_{[OIII,core]} as a surrogate for sigma_{*}. In addition, the substitute of sigma_{[OIII,core]} for sigma_{*} could be used only for low luminosity AGNs.
The coevolution between supermassive black holes (SMBHs) and their host galaxies has been proposed for more than a decade, albeit with little direct evidence about black hole accretion activities regulating galaxy star formation at z > 1. In this paper, we study the lifetimes of X-ray active galactic nuclei (AGNs) in UV-selected red sequence (RS), blue cloud (BC) and green valley (GV) galaxies, finding that AGN accretion activities are most prominent in GV galaxies at z ∼ 1.5–2, compared with RS and BC galaxies. We also compare AGN accretion timescales with typical color transition timescales of UV-selected galaxies. We find that the lifetime of GV galaxies at z ∼ 1.5–2 is very close to the typical timescale when the AGNs residing in them stay in the high-accretion-rate mode at these redshifts; for BC galaxies, the consistency between the color transition timescale and the black hole strong accretion lifetime is more likely to happen at lower redshifts (z < 1). Our results support the scenario where AGN accretion activities govern UV color transitions of host galaxies, making galaxies and their central SMBHs coevolve with each other.
Context.Extragalactic fast X-ray transients (FXRTs) are short flashes of X-ray photons of unknown origin that last a few seconds to hours.Aims.Our ignorance about their physical mechanisms and progenitor systems is due in part to the lack of clear multiwavelength counterparts in most cases, because FXRTs have only been identified serendipitously.Methods.We develop a systematic search for FXRTs in theChandraSource Catalog (Data Release 2.0; 169.6 Ms over 592.4 deg2, using only observations with |b|> 10° and before 2015), using a straightforward X-ray flare search algorithm and incorporating various multiwavelength constraints to rule out Galactic contamination and characterize the candidates.Results.We report the detection of 14 FXRT candidates from a parent sample of 214 701 sources. Candidates have peak 0.5–7 keV fluxes between 1 × 10−13and 2 × 10−10erg cm−2s−1andT90values from 4 to 48 ks. The sample can be subdivided into two groups: six “nearby” FXRTs that occurred withind ≲ 100 Mpc and eight “distant” FXRTs with likely redshifts ≳0.1. Three distant FXRT candidates exhibit light curves with a plateau (≈1–3 ks duration) followed by a power-law decay and X-ray spectral softening, similar to what was observed for the previously reported FXRT CDF-S XT2, a proposed magnetar-powered binary neutron star merger event. After applying completeness corrections, we calculate event rates for the nearby and distant samples of 53.7−15.1+22.6and 28.2−6.9+9.8deg−2yr−1, respectively.Conclusions.This novel sample ofChandra-detected extragalactic FXRT candidates, although modest in size, breaks new ground in terms of characterizing the diverse properties, nature, and possible progenitors of these enigmatic events.
The shape of the ambient circumnuclear medium (ACM) density profile can probe the history of accretion onto the central supermassive black holes in galaxies and the circumnuclear environment. However, due to the limitations of instrument resolution, the density profiles of the ACM for most galaxies remain largely unknown. In this work, we propose a novel method to measure the ACM density profile of active galactic nuclei (AGNs) by the equilibrium between the radiation pressure on the warm absorbers (WAs, a type of AGN outflow) and the drag pressure from the ACM. We study the correlation between the outflow velocity and ionization parameter of WAs in each of the five Seyfert 1 galaxies (NGC 3227, NGC 3783, NGC 4051, NGC 4593, and NGC 5548), inferring that the density profile of the ACM is between n ∝ r −1.7 and n ∝ r −2.15 (n is number density and r is distance) from 0.01 pc to parsec scales in these five AGNs. Our results indicate that the ACM density profile in Seyfert 1 galaxies is steeper than the prediction by the spherically symmetric Bondi accretion model and the simulated results of the hot accretion flow, but more in line with the prediction by the standard thin-disk model.
The relationship between quasars and their host galaxies provides clues on how supermassive black holes (SMBHs) and massive galaxies are jointly assembled. To elucidate this connection, we measure the structural and photometric properties of the host galaxies of ∼5000 Sloan Digital Sky Survey quasars at 0.2 < z < 1 using five-band (grizy) optical imaging from the Hyper Suprime-Cam Subaru Strategic Program. An automated analysis tool is used to forward-model the blended emission of the quasar as characterized by the point-spread function and the underlying host galaxy as a two-dimensional Sérsic profile. We find that quasars are preferentially hosted by massive star-forming galaxies with disklike light profiles. Furthermore, the size distribution of quasar hosts is broad at a given stellar mass and the average values exhibit a size–stellar mass relation as seen with inactive galaxies. In contrast, the sizes of quasar hosts are more compact than those of inactive star-forming galaxies on average, but not as compact as those of quiescent galaxies of similar stellar masses. This is true irrespective of quasar properties, including the bolometric luminosity, Eddington ratio, and black hole mass. These results are consistent with a scenario in which galaxies are concurrently fueling an SMBH and building their stellar bulge from a centrally concentrated gas reservoir. Alternatively, quasar hosts may be experiencing a compaction process in which stars from the disk and inflowing gas are responsible for growing the bulge. In addition, we confirm that the host galaxies of type 1 quasars have a bias of being closer to face-on systems, suggesting that galactic-scale dust can contribute to obscuring the broad-line region.
ABSTRACT Recent studies show that a universal relation between black hole (BH) growth and stellar mass (M⋆) or star formation rate (SFR) is an oversimplification of BH–galaxy coevolution, and that morphological and structural properties of host galaxies must also be considered. Particularly, a possible connection between BH growth and host-galaxy compactness was identified among star-forming (SF) galaxies. Utilizing ≈6300 massive galaxies with I814W < 24 at z < 1.2 in the Cosmic Evolution Survey (COSMOS) field, we perform systematic partial correlation analyses to investigate how sample-averaged BH accretion rate ($\rm \overline{BHAR}$) depends on host-galaxy compactness among SF galaxies, when controlling for morphology and M⋆ (or SFR). The projected central surface mass density within 1 kpc, Σ1, is utilized to represent host-galaxy compactness in our study. We find that the $\rm \overline{BHAR}$–Σ1 relation is stronger than either the $\rm \overline{BHAR}$–M⋆ or $\rm \overline{BHAR}$–SFR relation among SF galaxies, and this $\rm \overline{BHAR}$–Σ1 relation applies to both bulge-dominated galaxies and galaxies that are not dominated by bulges. This $\rm \overline{BHAR}$–Σ1 relation among SF galaxies suggests a link between BH growth and the central gas density of host galaxies on the kpc scale, which may further imply a common origin of the gas in the vicinity of the BH and in the central ∼kpc of the galaxy. This $\rm \overline{BHAR}$–Σ1 relation can also be interpreted as the relation between BH growth and the central velocity dispersion of host galaxies at a given gas content (i.e. gas mass fraction), indicating the role of the host-galaxy potential well in regulating accretion on to the BH.
We present single-epoch black hole mass (MBH) estimators based on the rest-frame ultraviolet (UV) Mg II 2798 angstrom and optical H beta 4861 angstrom emission lines. To enlarge the luminosity range of active galactic nuclei (AGNs), we combine the 31 reverberation-mapped AGNs with relatively low luminosities from Bahk et al., 47 moderateluminosity AGNs from Woo et al., and 425 high-luminosity AGNs from the Sloan Digital Sky Survey. The combined sample has a monochromatic luminosity at 5100 A in the range log lambda L-5100 similar to 41.3-46.5 erg s(-1), over the range 5.5 < log(M-BH/M-circle dot)<9.5. Based on the fiducial mass from the line dispersion or FWHM of H beta paired with the continuum luminosity at 5100 A, we calibrate the best-fit parameters in the black hole mass estimators using the Mg II line. We find that the differences in the line profiles between Mg II and H beta have significant effects on calibrating the UV M-BH estimators. By exploring the systematic discrepancy between the UV and optical M-BH estimators as a function of AGN properties, we suggest adding a correction term Delta M = -1.44 log(FWHMMgII/sigma(MgII)) + 0.33 in the equation for the UV mass estimator. We also find a similar to 0.1 dex bias in the MBH estimation due to the difference in the spectral slope in the range 2800-5200 angstrom. Depending on whether the selection of M-BH estimator is based on either line dispersion or FWHM and either continuum or line luminosity, the derived UV mass estimators show greater than or similar to 0.1 dex intrinsic scatter with respect to the fiducial H beta-based M-BH.
Active galactic nuclei (AGNs) have long been observed to "twinkle" (i.e., their brightness varies with time) on timescales from days to years in the UV/optical bands. Such AGN UV/optical variability is essential for probing the physics of supermassive black holes (SMBHs), the accretion disk, and the broad-line region. Here we show that the temperature fluctuations of an AGN accretion disk, which is magnetically coupled with the corona, can account for observed high-quality AGN optical light curves. We calculate the temperature fluctuations by considering the gas physics of the accreted matter near the SMBH. We find that the resulting simulated AGN UV/optical light curves share the same statistical properties as the observed ones as long as the dimensionless viscosity parameter $\alpha$, which is widely believed to be controlled by magnetohydrodynamic (MHD) turbulence in the accretion disk, is about $0.01$---$0.2$. Moreover, our model can simultaneously explain the larger-than-expected accretion disk sizes and the dependence of UV/optical variability upon wavelength for NGC 5548. Our model also has the potential to explain some other observational facts of AGN UV/optical variability, including the timescale-dependent bluer-when-brighter color variability and the dependence of UV/optical variability on AGN luminosity and black hole mass. Our results also demonstrate a promising way to infer the black-hole mass, the accretion rate, and the radiative efficiency, thereby facilitating understanding of the gas physics and MHD turbulence near the SMBH and its cosmic mass growth history by fitting the AGN UV/optical light curves in the era of time-domain astronomy.
The unified model of active galactic nuclei (AGNs) proposes that different AGN optical spectral types are caused by different viewing angles with respect to an obscuring "torus." Therefore, this model predicts that type 1 and type 2 AGNs should have similar host-galaxy properties. We investigate this prediction with 2463 X-ray-selected AGNs in the COSMOS field. We divide our sample into type 1 and type 2 AGNs based on their spectra, morphologies, and variability. We derive their host-galaxy stellar masses (M-*) through spectral energy distribution (SED) fitting, and we find that the hosts M-* of type 1 AGNs tend to be slightly smaller than those of type 2 AGNs by Delta(log M-*) over bar approximate to 0.2 dex (approximate to 4 sigma significance). Besides deriving star formation rates (SFRs) from SED fitting, we also utilize far-infrared (FIR) photometry and a stacking method to obtain FIR-based SFRs. We find that the SFRs of type 1 and type 2 sources are similar once their redshifts and X-ray luminosities are controlled. We also investigate the cosmic environment, and we find that the surface number densities (sub-Mpc) and cosmic-web environments (approximate to 1-10 Mpc) are similar for both populations. In summary, our analyses show that the host galaxies of type 1 and type 2 AGNs have similar SFRs and cosmic environments in general, but the former tend to have a lower M-* than the latter. The difference in M-* indicates that the AGN unification model is not strictly correct, and both host galaxy and torus may contribute to the optical obscuration of AGNs.
Two bright X-ray transients were reported from the Chandra Deep Field South (CDF-S) archival data, namely CDF-S XT1 and XT2. Whereas the nature of the former is not identified, the latter was suggested as an excellent candidate for a rapidly spinning magnetar born from a binary neutron star (BNS) merger. Here we propose a unified model to interpret both transients within the framework of the BNS merger magnetar model. According to our picture, CDF-S XT2 is observed from the “free zone” where the magnetar spindown powered X-ray emission escapes freely, whereas CDF-S XT1 originates from the “trapped zone” where the X-ray emission is initially blocked by the dynamical ejecta and becomes transparent after the ejecta is pushed to a distance where Thomson optical depth drops below unity. We fit the magnetar model to the light curves of both transients and derived consistent parameters for the two events, with magnetic field, initial spin period, and X-ray emission efficiency being ( B p = 10 16 G, P i = 1.2 ms, η = 0.001) and ( B p = 10 15.8 G, P i = 4.4 ms, η = 0.001) for XT1 and XT2, respectively. The “isotropic equivalent” ejecta mass of XT1 is M ej ∼ 10 −3 M ⊙ , while it is not constrained for XT2. Our results suggest that more extreme magnetar parameters are required to have XT1 detected from the trapped zone. The model parameters for both events are generally consistent with those derived from short gamma-ray burst (SGRB) X-ray plateau observations. The host-galaxy properties of both transients are also consistent with those of SGRBs. The event rate densities of both XT1 and XT2 are consistent with that of BNS mergers.