Search

We present the identification of changing-look active galactic nuclei (CL-AGNs) from the Dark Energy Spectroscopic Instrument First Data Release and Sloan Digital Sky Survey Data Release 16 at z \leq 0.9. To confirm the CL-AGNs, we utilize spectral flux calibration assessment via an [O\,{\sc iii}]-based calibration, pseudo-photometry examination, and visual inspection. This rigorous selection process allows us to compile a statistical catalog of 561 CL-AGNs, encompassing 527 $\rm H\beta$, 149$\rm H\alpha$, and 129 Mg II CL behaviors. In this sample, we find 1) a 283:278 ratio of turn-on to turn-off CL-AGNs. 2) the critical value for CL events is confirmed around Eddington ratio \sim 0.01. 3) a strong correlation between the change in the luminosity of the broad emission lines (BEL) and variation in the continuum luminosity, with Mg II and $\rm H\beta$ displaying similar responses during CL phases. 4) the Baldwin-Phillips-Terlevich diagram for CL-AGNs shows no statistically difference from the general AGN catalog. 5) five CL-AGNs are associated with asymmetrical mid-infrared flares, possibly linked to tidal disruption events. Given the large CL-AGNs and the stochastic sampling of spectra, we propose that some CL events are inherently due to typical AGN variability during low accretion rates, particularly for CL events of the singular BEL. Finally, we introduce a Peculiar CL phase, characterized by a gradual decline over decades in the light curve and the complete disappearance of entire BEL in faint spectra, indicative of a real transition in the accretion disk., Comment: Submitted to ApJS, comments welcome

Coupling between relativistic jets launched by accreting supermassive black holes and the surrounding gaseous media is a vital ingredient in galaxy evolution models. To constrain the environments in which this feedback takes place over cosmic time, we study the host halo properties of luminous low-frequency radio galaxies ($L_{150 \ \mathrm{MHz}} \gtrsim$ 25.25 W/Hz) selected with the International LOFAR Telescope out to $z \sim 2$ through tomographic clustering and cosmic microwave background lensing measurements. We find that these systems occupy halos characteristic of galaxy groups ($M_h = 10^{13} - 10^{14} h^{-1} M_{\odot}$), evolving at a rate consistent with the mean growth rate of halos over the past $\sim$10 Gyr. The coevolution of the clustering and the luminosity function reveals that the duty cycle of these systems is of order $\sim 10\%$ but has been mildly increasing since $z\sim 2$, while the duty cycle of quasars has been declining. We estimate the characteristic kinetic heating power injected by powerful jets per halo as a function of mass, and compare to the same quantity injected by quasar winds. We find that powerful jet heating dominates over quasar winds in halos $M_h \gtrsim 10^{13} h^{-1} M_{\odot}$ at $z < 2$. These results conform to the paradigm of galaxy evolution in which mechanical jet power feedback is the dominant heating mechanism of the gas content of groups and clusters., Comment: Accepted to ApJ

While reinforcement learning (RL) holds great potential for decision making in the real world, it suffers from a number of unique difficulties which often need specific consideration. In particular: it is highly non-stationary; suffers from high degrees of plasticity loss; and requires exploration to prevent premature convergence to local optima and maximize return. In this paper, we consider whether learned optimization can help overcome these problems. Our method, Learned Optimization for Plasticity, Exploration and Non-stationarity (OPEN), meta-learns an update rule whose input features and output structure are informed by previously proposed solutions to these difficulties. We show that our parameterization is flexible enough to enable meta-learning in diverse learning contexts, including the ability to use stochasticity for exploration. Our experiments demonstrate that when meta-trained on single and small sets of environments, OPEN outperforms or equals traditionally used optimizers. Furthermore, OPEN shows strong generalization across a distribution of environments and a range of agent architectures., Comment: AutoRL Workshop at ICML 2024

We study the impact of black hole nuclear activity on both the global and radial star formation rate (SFR) profiles in X-ray-selected active galactic nuclei (AGN) in the field of miniJPAS, the precursor of the much wider J-PAS project. Our sample includes 32 AGN with z < 0.3 detected via the XMM-Newton and Chandra surveys. For comparison, we assembled a control sample of 71 star-forming (SF) galaxies with similar magnitudes, sizes, and redshifts. To derive the global properties of both the AGN and the control SF sample, we used CIGALE to fit the spectral energy distributions derived from the 56 narrowband and 4 broadband filters from miniJPAS. We find that AGN tend to reside in more massive galaxies than their SF counterparts. After matching samples based on stellar mass and comparing their SFRs and specific SFRs (sSFRs), no significant differences appear. This suggests that the presence of AGN does not strongly influence overall star formation. However, when we used miniJPAS as an integral field unit (IFU) to dissect galaxies along their position angle, a different picture emerges. We find that AGN tend to be more centrally concentrated in mass with respect to SF galaxies. Moreover, we find a suppression of the sSFR up to 1Re and then an enhancement beyond 1Re , strongly contrasting with the decreasing radial profile of sSFRs in SF galaxies. This could point to an inside-out quenching of AGN host galaxies. These findings suggest that the reason we do not see differences on a global scale is because star formation is suppressed in the central regions and enhanced in the outer regions of AGN host galaxies. While limited in terms of sample size, this work highlights the potential of the upcoming J-PAS as a wide-field low-resolution IFU for thousands of nearby galaxies and AGN., Comment: 21 pages, 20 figures, Accepted for Publication in A&A

We present a catalog of hard X-ray serendipitous sources detected in the first 80 months of observations by the Nuclear Spectroscopic Telescope Array (NuSTAR). The NuSTAR serendipitous survey 80-month (NSS80) catalog has an unprecedented $\sim$ 62 Ms of effective exposure time over 894 unique fields (a factor of three increase over the 40-month catalog), with an areal coverage of $\sim $36 deg$^2$, larger than all NuSTAR extragalactic surveys. NSS80 provides 1274 hard X-ray sources in the $3-24$ keV band (822 new detections compared to the previous 40-month catalog). Approximately 76% of the NuSTAR sources have lower-energy ($<10$ keV) X-ray counterparts from Chandra, XMM-Newton, and Swift-XRT. We have undertaken an extensive campaign of ground-based spectroscopic follow-up to obtain new source redshifts and classifications for 427 sources. Combining these with existing archival spectroscopy provides redshifts for 550 NSS80 sources, of which 547 are classified. The sample is primarily composed of active galactic nuclei (AGN), detected over a large range in redshift ($z$ = 0.012-3.43), but also includes 58 spectroscopically confirmed Galactic sources. In addition, five AGN/galaxy pairs, one dual AGN system, one BL Lac candidate, and a hotspot of 4C 74.26 (radio quasar) have been identified. The median rest-frame $10-40$ keV luminosity and redshift of the NSS80 are $\langle{L_\mathrm{10-40 keV}}\rangle$ = 1.2 $\times$ 10$^{44}$ erg s$^{-1}$ and $\langle z \rangle = 0.56$. We investigate the optical properties and construct composite optical spectra to search for subtle signatures not present in the individual spectra, finding an excess of redder BL AGN compared to optical quasar surveys predominantly due to the presence of the host-galaxy and, at least in part, due to dust obscuration., Comment: Accepted for publication in ApJ:S. 57 pages, 32 figures

We adopt a Bayesian X-ray spectral approach to investigate the accretion properties of unobscured ($20<\log(N_{\rm H}/{\rm cm}^{-2}<22$) and obscured ($22< \log(N_{\rm H}/{\rm cm}^{-2}<24$) active galactic nuclei (AGN) to shed light on the orientation vs evolution scenarios for the origin of the obscuring material. For a sample of 3882 X-ray-selected AGN from the Chandra COSMOS Legacy, AEGIS and CDFS extragalactic surveys, we constrain their stellar masses, $M_\star$, intrinsic X-ray luminosities, $L_{\rm X}$, obscuring column densities, $N_{\rm H}$, and specific accretion rates $\lambda\propto L_{\rm X}/M_\star$. By combining these observables within a Bayesian non-parametric approach, we infer, for the first time, the specific accretion rate distribution (SARD) of obscured and unobscured AGN to $z\approx3$, i.e. the probability of a galaxy with mass $M_\star$ at redshift $z$ hosting an AGN with column density $N_{\rm H}$ and specific accretion rate $\lambda$. Our findings indicate that (1) both obscured and unobscured SARDs share similar shapes, shifting towards higher accretion rates with redshift, (2) unobscured SARDs exhibit a systematic offset towards higher $\lambda$ compared to obscured SARD for all redshift intervals, (3) the obscured AGN fraction declines sharply at $\log\lambda_{\rm break} \sim-2$ for $z <0.5$, but shifts to higher $\lambda$ values with increasing redshift, (4) the incidence of AGN within the theoretically unstable blow-out region of the $\lambda-N_{\rm H}$ plane increases with redshift. These observations provide compelling evidence for AGN "downsizing" and radiation-regulated nuclear-scale obscuration with an increasing host galaxy contribution towards higher redshifts., Comment: Submitted to MNRAS

Evolutionary models suggest that the initial growth phases of active galactic nuclei (AGN) are dust-enshrouded, and characterised by jet/wind outflows that should gradually clear the interstellar medium (ISM) in the host by heating and/or expelling the surrounding gas. eFEDSJ091157.4$+$014327 (z$\sim$0.6) was selected from X-ray samples for its characteristics that are similar to sources with outflows which include red, obscured and X-ray luminous. We aim to explore the environment around this red quasar and characterise kinematics within the system. We used spatially resolved spectroscopic data from Multi Unit Spectroscopic Explorer (MUSE) with an average seeing of 0.6" to construct flux, velocity and velocity dispersions maps. We found that the quasar is embedded in an interacting and merging system with three other galaxies $\sim$ 50 kpc from its nucleus. Spatially resolved kinematics reveal that the quasar has extended ionised outflows of up to 9.2 kpc with positive and negative velocities up to 1000 km s$^{-1}$ and -1200 km s$^{-1}$, respectively. The velocity dispersion (W$_{80}$) ranges from 600-1800 km s$^{-1}$. We associate the presence of turbulent and high-velocity components with the outflow. The total mass outflow rate is estimated to be $\sim$ 10 M$_{\odot}$ yr$^{-1}$ and kinetic power of 2x10$^{42}$ erg s$^{-1}$. The kinetic coupling efficiencies range from 0.01%-0.03% and the momentum boosts are $\sim$ 0.2. These low values indicate that the ionised outflow is not energetically relevant. These values don't align with the theoretical predictions of both radiation-pressure-driven outflows and energy-conserving mechanisms. However, note that our results are based only on the ionised phase while theoretical predictions are multiphase. Moreover, the mass loading factor of $\sim$ 5 is an indication that these outflows are more likely AGN-driven than star formation-driven., Comment: Version after addressing referee comments, accepted in Astronomy and Astrophysics(A&A), 17 pages, 12 figures, 1 table

Recent deep and wide radio surveys extend the studies for radio-excess active galactic nuclei (radio-AGNs) to lower luminosities and higher redshifts, providing new insights into the abundance and physical origin of radio-AGNs. Here we focus on the cosmic evolution, physical properties and AGN-host galaxy connections of radio-AGNs selected from a sample of ~ 500,000 galaxies at 0 < z < 4 in GOODS-N, GOODS-S, and COSMOS fields. Combining deep radio data with multi-band, de-blended far-infrared (FIR) and sub-millimeter data, we identify 1162 radio-AGNs through radio excess relative to the FIR-radio relation. We study the cosmic evolution of 1.4 GHz radio luminosity functions (RLFs) for star-forming galaxies (SFGs) and radio-AGNs, which are well described by a pure luminosity evolution of $L_*\propto (1+z)^{-0.31z+3.41}$ and a pure density evolution of $\Phi_*\propto (1+z)^{-0.80z+2.88}$, respectively. We derive the turnover luminosity above which the number density of radio-AGNs surpasses that of SFGs. This crossover luminosity increases as increasing redshift, from $10^{22.9}$ W Hz$^{-1}$ at z ~ 0 to $10^{25.2}$ W Hz$^{-1}$ at z ~ 4. At full redshift range (0 < z < 4), we further derive the probability ($p_{radio}$) of SFGs and quiescent galaxies (QGs) hosting a radio-AGN as a function of stellar mass ($M_*$), radio luminosity ($L_R$), and redshift (z), which yields $p_{radio}\propto (1+z)^{3.54}M_*^{1.02}L_R^{-0.90}$ for SFGs, and $p_{radio}\propto (1+z)^{2.38}M_*^{1.39}L_R^{-0.60}$ for QGs, respectively. It indicates that radio-AGNs in QGs prefer to reside in more massive galaxies with larger $L_R$ than those in SFGs, and radio-AGN fraction increases towards higher redshift in both SFGs and QGs with a more rapid increase in SFGs. Further, we find that the radio-AGN fraction depends on accretion states of BHs and redshift in SFGs, while in QGs it also depends on BH (or galaxy) mass., Comment: Accepted by Astronomy & Astrophysics (A&A), 25 pages, 17 figures, 5 tables

Since the discovery of the cosmic X-ray background (CXB), astronomers have strived to understand the accreting supermassive black holes (SMBHs) contributing to its peak in the 10-40 keV band. Existing soft X-ray telescopes could study this population up to only 10 keV, and, while NuSTAR (focusing on 3--24 keV) made great progress, it also left significant uncertainties in characterizing the hard X-ray population, crucial for calibrating current population synthesis models. This paper presents an in-depth analysis of simulations of two extragalactic surveys (deep and wide) with the High-Energy X-ray Probe (HEX-P), each observed for 2 Ms. Applying established source detection techniques, we show that HEX-P surveys will reach a flux of $\sim$10$^{-15}$ erg s$^{-1}$ cm$^{-2}$ in the 10-40 keV band, an order of magnitude fainter than current NuSTAR surveys. With the large sample of new hard X-ray detected sources ($\sim2000$), we showcase HEX-P's ability to resolve more than 80% of the CXB up to 40 keV into individual sources. The expected precision of HEX-P's resolved background measurement will allow us to distinguish between population synthesis models of SMBH growth. HEX-P leverages accurate broadband (0.5-40 keV) spectral analysis and the combination of soft and hard X-ray colors to provide obscuration constraints even for the fainter sources, with the overall objective of measuring the Compton-thick fraction. With unprecedented sensitivity in the 10--40 keV band, HEX-P will explore the hard X-ray emission from AGN to flux limits never reached before, thus expanding the parameter space for serendipitous discoveries. Consequently, it is plausible that new models will be needed to capture the population HEX-P will unveil., Comment: 28 pages, 13 Figures, 2 tables; submitted to Frontiers in Astronomy and Space Sciences; minor edits in the text, abstract and figure labels

The origin of obscuration in Active Galactic Nuclei (AGN) is still an open debate. In particular, it is unclear what drives the relative contributions to the line-of-sight column densities from galaxy-scale and torus-linked obscuration. The latter source is expected to play a significant role in Unification Models, while the former is thought to be relevant in both Unification and Evolutionary Models. In this work, we make use of a combination of cosmological semi-analytic models and semi-empirical prescriptions for the properties of galaxies and AGN, to study AGN obscuration. We consider a detailed object-by-object modelling of AGN evolution, including different AGN light curves (LCs), gas density profiles, and also AGN feedback-induced gas cavities. Irrespective of our assumptions on specific AGN LC or galaxy gas fractions, we find that, on the strict assumption of an exponential profile for the gas component, galaxy-scale obscuration alone can hardly reproduce the fraction of $\log (N_{\rm H}/$cm$^{-2}) \geq 24$ sources at least at $z\lesssim3$. This requires an additional torus component with a thickness that decreases with luminosity to match the data. The torus should be present in all evolutionary stages of a visible AGN to be effective, although galaxy-scale gas obscuration may be sufficient to reproduce the obscured fraction with $22<\log (N_{\rm H}/$cm$^{-2})<24$ (Compton-thin, CTN) if we assume extremely compact gas disc components. The claimed drop of CTN fractions with increasing luminosity does not appear to be a consequence of AGN feedback, but rather of gas reservoirs becoming more compact with decreasing stellar mass., Comment: MNRAS, accepted, 19 pages, 15 figures, 3 appendices

In the standard quasar model, the accretion disk obscuration is due to the canonical dusty torus. Here, we argue that a substantial part of the quasar obscuration can come from the interstellar medium (ISM) when the quasars are embedded in compact starbursts. We use an obscuration-unbiased sample of 578 infrared (IR) quasars at $z\approx 1-3$ and archival ALMA submillimeter host galaxy sizes to investigate the ISM contribution to the quasar obscuration. We calculate SFR and ISM column densities for the IR quasars and a control sample of submillimeter galaxies (SMGs) not hosting quasar activity and show that: (1) the quasar obscured fraction is constant up to $\rm SFR\approx 300 \: M_{\odot} \: yr^{-1}$, and then increases towards higher SFR, suggesting that the ISM obscuration plays a significant role in starburst host galaxies, and (2) at $\rm SFR\gtrsim 300 \: M_{\odot} \: yr^{-1}$, the SMGs and IR quasars have similarly compact submillimeter sizes ($R_{\rm e}\approx 0.5-3\rm \: kpc$) and, consequently, the ISM can heavily obscure the quasar, even reaching Compton-thick ($N_{\rm H}>10^{24} \rm \: cm^{-2}$) levels in extreme cases. Based on our results, we infer that $\approx 10-30\%$ of the IR quasars with $\rm SFR\gtrsim 300 \: M_{\odot} \: yr^{-1}$ are obscured solely by the ISM., Comment: Accepted for publication in MNRAS Letters

In the context of an evolutionary model, the outflow phase of an Active Galactic Nuclei (AGN) occurs at the peak of its activity, once the central SMBH is massive enough to generate sufficient power to counterbalance the potential well of the host galaxy. This phase plays a vital role in galaxy evolution. We aim to apply various selection methods to isolate powerful AGNs in the feedback phase, trace and characterise their outflows, and explore the link between AGN luminosity and outflow properties. We applied a combination of methods to the eROSITA Final Equatorial Depth survey (eFEDS) catalogue and isolated ~1400 candidates at z>0.5 out of ~11750 AGNs (~12\%). We tested the robustness of our selection on the small subsample of 50 sources with available good quality SDSS spectra at 0.51\%. We find that the majority of the outflows have momentum flux ratios lower than 20 which rules out an energy-conserving nature. Our present work points to the unequivocal existence of a rather short AGN outflow phase, paving the way towards a new avenue to dissect AGN outflows in large samples within eROSITA and beyond., Comment: Accepted for publication in Astronomy and Astrophysics(A&A), 24 pages, 18 figures, 5 tables

Changing-look Active Galactic Nuclei (CL AGN) can be generally confirmed by the emergence (turn-on) or disappearance (turn-off) of broad emission lines, associated with a transient timescale (about $100\sim5000$ days) that is much shorter than predicted by traditional accretion disk models. We carry out a systematic CL AGN search by cross-matching the spectra coming from the Dark Energy Spectroscopic Instrument and the Sloan Digital Sky Survey. Following previous studies, we identify CL AGN based on $\rm{H}\alpha $, $\rm{H}\beta$, and Mg\,{\sc ii} at $z\leq0.75$ and Mg\,{\sc ii}, C\,{\sc iii}], and C\,{\sc iv} at $z>0.75$. We present 130 CL AGN based on visual inspection and three selection criteria, including 2 $\rm{H}\alpha$, 45 $\rm{H}\beta$, 38 Mg\,{\sc ii}, 61 C\,{\sc iii}], and 10 C\,{\sc iv} CL AGN. Twenty cases show simultaneous appearances/disappearances of two broad emission lines while three AGN exhibit the concurrent appearance of three broad emission lines. We also present 91 CL AGN candidates with significant flux variation of broad emission lines but remaining strong broad components. In the confirmed CL AGN, 42 cases show additional CL candidate features for different lines. In this paper, we find 1) a 95:35 ratio of a turn-on to turn-off CL AGN; 2) the highest redshift CL AGN ($z=3.56$) ever discovered; 3) an upper limit transition timescale ranging from 244 to 5762 days in the rest-frame; 4) the majority of CL AGN follow the bluer-when-brighter trend. Our results greatly increase the current CL census ($30\sim50\%$) and would be conducive to explore the underlying physical mechanism.

We present findings of the detection of Magnesium II (Mg II, {\lambda} = 2796, 2803 {\AA}) absorbers from the early data release of the Dark Energy Spectroscopic Instrument (DESI). DESI is projected to obtain spectroscopy of approximately 3 million quasars (QSOs), of which over 99% are anticipated to be at redshifts greater than z > 0.3, such that DESI would be able to observe an associated or intervening Mg II absorber illuminated by the background QSO. We have developed an autonomous supplementary spectral pipeline that detects these systems through an initial line-fitting process and then confirms the line properties using a Markov Chain Monte Carlo sampler. Based upon a visual inspection of the resulting systems, we estimate that this sample has a purity greater than 99%. We have also investigated the completeness of our sample in regard to both the signal-to-noise properties of the input spectra and the rest-frame equivalent width (W0) of the absorber systems. From a parent catalog containing 83,207 quasars, we detect a total of 23,921 Mg II absorption systems following a series of quality cuts. Extrapolating from this occurrence rate of 28.8% implies a catalog at the completion of the five-year DESI survey that will contain over eight hundred thousand Mg II absorbers. The cataloging of these systems will enable significant further research because they carry information regarding circumgalactic medium environments, the distribution of intervening galaxies, and the growth of metallicity across the redshift range 0.3 < z < 2.5., Comment: 17 pages, 8 figures

Millions of quasar spectra will be collected by the Dark Energy Spectroscopic Instrument (DESI), leading to a four-fold increase in the number of known quasars. High accuracy quasar classification is essential to tighten constraints on cosmological parameters measured at the highest redshifts DESI observes ($z>2.0$). We present the spectral templates for identification and redshift estimation of quasars in the DESI Year 1 data release. The quasar templates are comprised of two quasar eigenspectra sets, trained on spectra from the Sloan Digital Sky Survey. The sets are specialized to reconstruct quasar spectral variation observed over separate yet overlapping redshift ranges and, together, are capable of identifying DESI quasars from $0.05 < z <7.0$. The new quasar templates show significant improvement over the previous DESI quasar templates regarding catastrophic failure rates, redshift precision and accuracy, quasar completeness, and the contamination fraction in the final quasar sample., Comment: to be published in Astronomical Journal; 21 pages, 8 figures

We analyze the physical properties of 8 X-ray selected active galactic nuclei (AGN) and one candidate protoquasar system (ADF22A1) in the $z = 3.09$ SSA22 protocluster by fitting their X-ray-to-IR spectral energy distributions (SEDs) using our SED fitting code, Lightning. We recover star formation histories (SFH) for 7 of these systems which are well-fit by composite stellar population plus AGN models. We find indications that 4/9 of the SSA22 AGN systems we study have host galaxies below the main sequence, with $\rm SFR/SFR_{MS} \leq -0.4$. The remaining SSA22 systems, including ADF22A1, are consistent with obscured supermassive black hole (SMBH) growth in star forming galaxies. We estimate the SMBH accretion rates and masses, and compare the properties and SFH of the 9 protocluster AGN systems with X-ray detected AGN candidates in the Chandra Deep Fields (CDF), finding that the distributions of SMBH growth rates, star formation rates, SMBH masses, and stellar masses for the protocluster AGN are consistent with field AGN. We constrain the ratio between the sample-averaged SSA22 SMBH mass and CDF SMBH mass to $<1.41$. While the AGN are located near the density peaks of the protocluster, we find no statistically significant trends between the AGN or host galaxy properties and their location in the protocluster. We interpret the similarity of the protocluster and field AGN populations together with existing results as suggesting that the protocluster and field AGN co-evolve with their hosts in the same ways, while AGN-triggering events are more likely in the protocluster., Comment: 27 pages, 14 figures, 8 tables. Updated with published version

Active dwarf galaxies are important because they contribute to the evolution of dwarf galaxies and can reveal their hosted massive black holes. However, the sample size of such sources beyond the local universe is still highly limited. In this work, we search for active dwarf galaxies in the recently completed XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS). XMM-SERVS is currently the largest medium-depth X-ray survey covering 13 $\mathrm{deg}^2$ in three extragalactic fields, which all have well-characterized multi-wavelength information. After considering several factors that may lead to misidentifications, we identify 73 active dwarf galaxies at $z<1$, which constitutes the currently largest X-ray-selected sample beyond the local universe. Our sources are generally less obscured than predictions based on the massive-AGN (active galactic nucleus) X-ray luminosity function and have a low radio-excess fraction. We find that our sources reside in similar environments to inactive dwarf galaxies. We further quantify the accretion distribution of the dwarf-galaxy population after considering various selection effects and find that it decreases with X-ray luminosity, but redshift evolution cannot be statistically confirmed. Depending upon how we define an AGN, the active fraction may or may not show a strong dependence on stellar mass. Their Eddington ratios and X-ray bolometric corrections significantly deviate from the expected relation, which is likely caused by several large underlying systematic biases when estimating the relevant parameters for dwarf galaxies. Throughout this work, we also highlight problems in reliably measuring photometric redshifts and overcoming strong selection effects for distant active dwarf galaxies., Comment: 34 pages, 21 figures, 6 tables, accepted for publication in ApJ

We perform X-ray spectral analyses to derive characteristics (e.g., column density, X-ray luminosity) of $\approx$10,200 active galactic nuclei (AGNs) in the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS), which was designed to investigate the growth of supermassive black holes across a wide dynamic range of cosmic environments. Using physical torus models (e.g., Borus02) and a Bayesian approach, we uncover 22 representative Compton-thick (CT; $N_{\rm H} \;>\; 1.5\times10^{24}\; \rm cm^{-2}$) AGN candidates with good signal-to-noise ratios as well as a large sample of 136 heavily obscured AGNs. We also find an increasing CT fraction (\fct ) from low ($z<0.75$) to high ($z>0.75$) redshift. Our CT candidates tend to show hard X-ray spectral shapes and dust extinction in their SED fits, which may shed light on the connection between AGN obscuration and host-galaxy evolution., Comment: 12 pages, 9 figures, accepted for publication in ApJ

We present the first result in exploring the gaseous halo and galaxy correlation using the Dark Energy Spectroscopic Instrument (DESI) survey validation data in the Cosmic Evolution Survey (COSMOS) and Hyper Suprime-Cam (HSC) field. We obtain the multiphase gaseous halo properties in the circumgalactic medium (CGM) by using 115 quasar spectra (S/N > 3). We detect MgII absorption at redshift 0.6 < z < 2.5, CIV absorption at 1.6 < z < 3.6, and HI absorption associated with the MgII and CIV. By cross-matching the COSMOS2020 catalog, we identify the MgII and CIV host galaxies in ten quasar fields at 0.9 < z < 3.1. We find that within the impact parameter of 250 kpc, a tight correlation is seen between strong MgII equivalent width and the host galaxy star formation rate. The covering fraction fc of strong MgII selected galaxies, which is the ratio of absorbing galaxy in a certain galaxy population, shows significant evolution in the main-sequence galaxies and marginal evolution in all the galaxy populations within 250 kpc at 0.9 < z < 2.2. The fc increase in the main-sequence galaxies likely suggests the co-evolution of strong MgII absorbing gas and the main-sequence galaxies at the cosmic noon. Furthermore, several MgII and CIV absorbing gas is detected out of the galaxy virial radius, tentatively indicating the feedback produced by the star formation and/or the environmental effects., Comment: 24 pages, 12 figures, 6 tables, accepted for publication in ApJ

We report the first results of a high-redshift ($z$ >~ 5) quasar survey using the Dark Energy Spectroscopic Instrument (DESI). As a DESI secondary target program, this survey is designed to carry out a systematic search and investigation of quasars at $z$ >~ 5, up to redshift 6.8. The target selection is based on the DESI Legacy Imaging Surveys (the Legacy Surveys) DR9 photometry, combined with the Pan-STARRS1 data and $J$-band photometry from public surveys. A first quasar sample has been constructed from the DESI Survey Validation 3 (SV3) and first-year observations until May 2022. This sample includes more than 400 new quasars at redshift 4.7 <= $z$ < 6.6, down to 21.5 magnitude in the $z$ band, discovered from 35% of the entire target sample. Remarkably, there are 220 new quasars identified at $z$ >= 5, more than one third of existing quasars previously published at this redshift. The observations so far result in an average success rate of 23% at $z$ > 4.7. The current spectral dataset has already allowed analysis of interesting individual objects (e.g., quasars with damped Ly$\alpha$ absorbers and broad absorption line features), and statistical analysis will follow the survey's completion. A set of science projects will be carried out leveraging this program, including quasar luminosity function, quasar clustering, intergalactic medium, quasar spectral properties, intervening absorbers, and properties of early supermassive black holes. Additionally, a sample of 38 new quasars at $z$ ~ 3.8-5.7 discovered from a pilot survey in the DESI SV1 is also published in this paper., Comment: 28 pages, 9 figures, and 2 tables; Published in ApJS

Obscuration in quasars may arise from steep viewing angles along the dusty torus, or instead may represent a distinct phase of supermassive black hole growth. We test these scenarios by probing the host dark matter halo environments of $\sim 1.4$ million WISE-selected obscured and unobscured quasars at $\langle z \rangle = 1.4$ using angular clustering measurements as well as cross-correlation measurements of quasar positions with the gravitational lensing of the cosmic microwave background (CMB). We interpret these signals within a halo occupation distribution (HOD) framework to conclude that obscured systems reside in more massive effective halos ($ \sim 10^{12.9} h^{-1} M_{\odot}$) than their unobscured counterparts ($ \sim 10^{12.6} h^{-1} M_{\odot}$), though we do not detect a difference in the satellite fraction. We find excellent agreement between the clustering and lensing analyses and show that this implies the observed difference is robust to uncertainties in the obscured quasar redshift distribution, highlighting the power of combining angular clustering and weak lensing measurements. This finding appears in tension with models that ascribe obscuration exclusively to orientation of the dusty torus along the line-of-sight, and instead may be consistent with the notion that some obscured quasars are attenuated by galaxy-scale or circumnuclear material during an evolutionary phase., Comment: Accepted to The Astrophysical Journal

A complete census of Active Galactic Nuclei (AGN) is a prerequisite for understanding the growth of supermassive black holes across cosmic time. A significant challenge toward this goal is the whereabouts of heavily obscured AGN that remain uncertain. This paper sets new constraints on the demographics of this population by developing a methodology that combines X-ray spectral information with priors derived from multiwavelength observations. We select X-ray AGN in the Chandra COSMOS Legacy survey and fit their $2.2-500\mu m$ spectral energy distributions with galaxy and AGN templates to determine the mid-infrared ($6\mu m$) luminosity of the AGN component. Empirical correlations between X-ray and $6\mu m$ luminosities are then adopted to infer the intrinsic accretion luminosity at X-rays for individual AGN. This is used as prior information in our Bayesian X-ray spectral analysis to estimate physical properties, such as line-of-sight obscuration. Our approach breaks the degeneracies between accretion luminosity and obscuration that affect X-ray spectral analysis, particularly for the most heavily obscured (Compton-Thick) AGN with low photon counts X-ray spectra. The X-ray spectral results are then combined with the selection function of the Chandra COSMOS Legacy survey to derive the AGN space density and a Compton-Thick fraction of $21.0^{+16.1}_{-9.9}\%$ at redshifts $z<0.5$. At higher redshift, our analysis suggests upper limits to the Compton-Thick AGN fraction of $\le 40\%$. These estimates are at the low end of the range of values determined in the literature and underline the importance of multiwavelength approaches for tackling the challenge of heavily obscured AGN demographics., Comment: Accepted for publication at MNRAS, 21 pages, 20 figures, 6 tables

Recent studies have revealed a strong relation between sample-averaged black-hole (BH) accretion rate (BHAR) and star formation rate (SFR) among bulge-dominated galaxies, i.e., "lockstep" BH-bulge growth, in the distant universe. This relation might be closely related to the BH-bulge mass correlation observed in the local universe. To understand further BH-bulge coevolution, we present ALMA CO(2-1) or CO(3-2) observations of 7 star-forming bulge-dominated galaxies at z=0.5-2.5. Using the ALMA data, we detect significant ($>3\sigma$) CO emission from 4 objects. For our sample of 7 galaxies, we measure (or constrain with upper limits) their CO line fluxes and estimate molecular gas masses ($M_{gas}$). We also estimate their stellar masses ($M_{star}$) and SFRs by modelling their spectral energy distributions (SEDs). Using these physical properties, we derive the gas-depletion timescales ($t_{dep} = M_{gas}/SFR$) and compare them with the bulge/BH growth timescales ($t_{grow} = M_{star}/SFR \sim M_{BH}/BHAR$). Our sample generally has $t_{dep}$ shorter than $t_{grow}$ by a median factor of $\gtrsim 4$, indicating that the cold gas will be depleted before significant bulge/BH growth takes place. This result suggests that the BH-bulge lockstep growth is mainly responsible for maintaining their mass relation, not creating it. We note that our sample is small and limited to $z<2.5$; JWST and ALMA will be able to probe to higher redshifts in the near future., Comment: 15 pages, 7 figures, and 3 tables. Accepted by ApJ

Our current understanding of the cosmic star formation history at z>3 is primarily based on UV-selected galaxies (i.e., LBGs). Recent studies of H-dropouts have revealed that we may be missing a large proportion of star formation that is taking place in massive galaxies at z>3. In this work, we extend the H-dropout criterion to lower masses to select optically dark/faint galaxies (OFGs), in order to complete the census between LBGs and H-dropouts. Our criterion (H> 26.5 mag & [4.5] < 25 mag) combined with a de-blending technique is designed to select not only extremely dust-obscured massive galaxies but also normal star-forming galaxies. In total, we identified 27 OFGs at z_phot > 3 (z_med=4.1) in the GOODS-ALMA field, covering a wide distribution of stellar masses with log($M_{\star}$/$M_{\odot}$) = 9.4-11.1. We find that up to 75% of the OFGs with log($M_{\star}$/$M_{\odot}$) = 9.5-10.5 were neglected by previous LBGs and H-dropout selection techniques. After performing stacking analyses, the OFGs exhibit shorter gas depletion timescales, slightly lower gas fractions, and lower dust temperatures than typical star-forming galaxies. Their SFR_tot (SFR_ IR+SFR_UV) is much larger than SFR_UVcorr (corrected for dust extinction), with SFR_tot/SFR_UVcorr = $8\pm1$, suggesting the presence of hidden dust regions in the OFGs that absorb all UV photons. The average dust size measured by a circular Gaussian model fit is R_e(1.13 mm)=1.01$\pm$0.05 kpc. We find that the cosmic SFRD at z>3 contributed by massive OFGs is at least two orders of magnitude higher than the one contributed by equivalently massive LBGs. Finally, we calculate the combined contribution of OFGs and LBGs to the cosmic SFRD at z=4-5 to be 4 $\times$ 10$^{-2}$ $M_{\odot}$ yr$^{-1}$Mpc$^{-3}$, which is about 0.15 dex (43%) higher than the SFRD derived from UV-selected samples alone at the same redshift., Comment: 26 pages, 14 figures, 8 tables. Accepted for publication in A&A

To understand the Active Galactic Nuclei (AGN) phenomenon and their impact on the evolution of galaxies, a complete AGN census is required; however, finding heavily obscured AGNs is observationally challenging. Here we use the deep and extensive multi-wavelength data in the COSMOS field to select a complete sample of 578 infrared (IR) quasars ($L_{\rm AGN,IR}>10^{45}\rm \: erg\: s^{-1}$) at $z<3$, with minimal obscuration bias, using detailed UV-to-far IR spectral energy distribution (SED) fitting. We complement our SED constraints with X-ray and radio observations to further investigate the properties of the sample. Overall, 322 of the IR quasars are detected by Chandra and have individual X-ray spectral constraints. From a combination of X-ray stacking and $L_{\rm 2-10\rm keV}$ - $L_{\rm 6\: \mu m}$ analyses, we show that the majority of the X-ray faint and undetected quasars are heavily obscured (many are likely Compton thick), highlighting the effectiveness of the mid-IR band to find obscured AGNs. We find that 355 ($\approx$61%) IR quasars are obscured ($N_{\rm H}>10^{22}\rm \: cm^{-2}$) and identify differences in the average properties between the obscured and unobscured quasars: (1) obscured quasars have star-formation rates $\approx 3$ times higher than unobscured systems for no significant difference in stellar mass and (2) obscured quasars have stronger radio emission than unobscured systems, with a radio-loudness parameter $\approx 0.2 \rm \: dex$ higher. These results are inconsistent with a simple orientation model but in general agreement with either extreme host-galaxy obscuration towards the obscured quasars or a scenario where obscured quasars are an early phase in the evolution of quasars., Comment: Accepted for publication in MNRAS. 22 pages, 15 figures, 4 tables

In 2021 May, the Dark Energy Spectroscopic Instrument (DESI) began a 5 yr survey of approximately 50 million total extragalactic and Galactic targets. The primary DESI dark-time targets are emission line galaxies (ELGs), luminous red galaxies (LRGs) and quasars (QSOs). In bright time, DESI will focus on two surveys known as the Bright Galaxy Survey (BGS) and the Milky Way Survey (MWS). DESI also observes a selection of "secondary" targets for bespoke science goals. This paper gives an overview of the publicly available pipeline (desitarget) used to process targets for DESI observations. Highlights include details of the different DESI survey targeting phases, the targeting ID (TARGETID) used to define unique targets, the bitmasks used to indicate a particular type of target, the data model and structure of DESI targeting files, and examples of how to access and use the desitarget code base. This paper will also describe "supporting" DESI target classes, such as standard stars, sky locations, and random catalogs that mimic the angular selection function of DESI targets. The DESI target selection pipeline is complex and sizable; this paper attempts to summarize the most salient information required to understand and work with DESI targeting data., Comment: AJ, accepted, 27 pages, 4 figures, 10 tables, one of a suite of 8 papers detailing targeting for DESI. Minor textual updates to better match the final, accepted version. Also added two missing co-authors

A key component of the Dark Energy Spectroscopic Instrument (DESI) survey validation (SV) is a detailed visual inspection (VI) of the optical spectroscopic data to quantify key survey metrics. In this paper we present results from VI of the quasar survey using deep coadded SV spectra. We show that the majority (~70%) of the main-survey targets are spectroscopically confirmed as quasars, with ~16% galaxies, ~6% stars, and ~8% low-quality spectra lacking reliable features. A non-negligible fraction of the quasars are misidentified by the standard spectroscopic pipeline but we show that the majority can be recovered using post-pipeline "afterburner" quasar-identification approaches. We combine these "afterburners" with our standard pipeline to create a modified pipeline to improve the overall quasar yield. At the depth of the main DESI survey both pipelines achieve a good-redshift purity (reliable redshifts measured within 3000 km/s) of ~99%; however, the modified pipeline recovers ~94% of the visually inspected quasars, as compared to ~86% from the standard pipeline. We demonstrate that both pipelines achieve an median redshift precision and accuracy of ~100 km/s and ~70 km/s, respectively. We constructed composite spectra to investigate why some quasars are missed by the standard spectroscopic pipeline and find that they are more host-galaxy dominated (i.e., distant analogs of "Seyfert galaxies") and/or dust reddened than the standard-pipeline quasars. We also show example spectra to demonstrate the overall diversity of the DESI quasar sample and provide strong-lensing candidates where two targets contribute to a single spectrum., Comment: Astronomical journal (in press). 26 pages, 15 figures, 9 tables, one of a suite of 8 papers detailing targeting for DESI. Figure data available from Zenodo (see paper for details)

The Dark Energy Spectroscopic Instrument (DESI) Survey has obtained a set of spectroscopic measurements of galaxies to validate the final survey design and target selections. To assist in these tasks, we visually inspect (VI) DESI spectra of approximately 2,500 bright galaxies, 3,500 luminous red galaxies (LRGs), and 10,000 emission line galaxies (ELGs), to obtain robust redshift identifications. We then utilize the VI redshift information to characterize the performance of the DESI operation. Based on the VI catalogs, our results show that the final survey design yields samples of bright galaxies, LRGs, and ELGs with purity greater than $99\%$. Moreover, we demonstrate that the precision of the redshift measurements is approximately 10 km/s for bright galaxies and ELGs and approximately 40 km/s for LRGs. The average redshift accuracy is within 10 km/s for the three types of galaxies. The VI process also helps improve the quality of the DESI data by identifying spurious spectral features introduced by the pipeline. Finally, we show examples of unexpected real astronomical objects, such as Ly$\alpha$ emitters and strong lensing candidates, identified by VI. These results demonstrate the importance and utility of visually inspecting data from incoming and upcoming surveys, especially during their early operation phases., Comment: 22 pages, 14 figures, 3 tables, one of a suite of 8 papers detailing targeting for DESI. ApJ accepted version with minor textual updates

The Dark Energy Spectroscopic Instrument (DESI) survey will measure large-scale structures using quasars as direct tracers of dark matter in the redshift range 0.92.1. We present several methods to select candidate quasars for DESI, using input photometric imaging in three optical bands (g, r, z) from the DESI Legacy Imaging Surveys and two infrared bands (W1, W2) from the Wide-field Infrared Explorer (WISE). These methods were extensively tested during the Survey Validation of DESI. In this paper, we report on the results obtained with the different methods and present the selection we optimized for the DESI main survey. The final quasar target selection is based on a Random Forest algorithm and selects quasars in the magnitude range 16.599% purity for a nominal effective exposure time of ~1000s. With a 310 per sq. deg. target density, the main selection allows DESI to select more than 200 QSOs per sq. deg. (including 60 quasars with z>2.1), exceeding the project requirements by 20%. The redshift distribution of the selected quasars is in excellent agreement with quasar luminosity function predictions., Comment: 21 pages, 21 figures, typos corrected, references added