Showing 1–13 of 13 results for author: Laloux, B

Genuine Retrieval of the AGN Host Stellar Population (GRAHSP)

Authors: Johannes Buchner, Hattie Starck, Mara Salvato, Hagai Netzer, Zsofi Igo, Brivael Laloux, Antonis Georgakakis, Isabelle Gauger, Anna Olechowska, Nicolas Lopez, Suraj D Shankar, Junyao Li, Kirpal Nandra, Andrea Merloni

Abstract: The assembly and co-evolution of supermassive black holes (SMBH) and their host galaxy stellar population is a key open questions in galaxy evolution. Stellar mass ($M_\star$) and star formation rate (SFR), are inferred by modeling the spectral energy distribution (SED). For galaxies triggering SMBH activity, the active galactic nucleus (AGN) contaminates the light at all wavelengths, hampering th… ▽ More The assembly and co-evolution of supermassive black holes (SMBH) and their host galaxy stellar population is a key open questions in galaxy evolution. Stellar mass ($M_\star$) and star formation rate (SFR), are inferred by modeling the spectral energy distribution (SED). For galaxies triggering SMBH activity, the active galactic nucleus (AGN) contaminates the light at all wavelengths, hampering the inference of galaxy parameters. Incomplete AGN templates can lead to systematic overestimates of the stellar mass, biasing our understanding of AGN-galaxy co-evolution. This challenge has gained further impetus with the advent of sensitive wide-area surveys with millions of luminous AGN, including by eROSITA, Euclid and LSST. We aim to estimate the accuracy and bias of AGN host galaxy parameters and improve upon existing techniques. This work makes two contributions: 1) a new SED fitting code, GRAHSP, with a flexible, empirically motivated AGN model including a power law continuum emission lines, a FeII forest and a flexible infrared torus. We verify that our model reproduces published X-ray to infrared SEDs of AGN to better than 20\% accuracy. A fully Bayesian fit with nested sampling includes uncertainties in the model and the data, making the inference highly robust. 2) we created a benchmark photometric dataset where pure quasars are merged with non-AGN pure galaxies into a hybrid (Chimera) object but with known galaxy and AGN properties. Comparing the true and retrieved $M_\star$, SFR and AGN luminosities shows that previous codes systematically over-estimate $M_\star$ and SFR by 0.5 dex with a wide scatter of 0.7 dex, at AGN luminosities above 10^44 erg/s. In contrast, GRAHSP shows no bias on $M_\star$ and SFR. GRAHSP also estimates more realistic uncertainties. GRAHSP enables characterization of the environmental conditions conducive to black hole growth. (abridged) △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: version resubmitted to A&A after a first positive referee report

The miniJPAS Survey: The radial distribution of star formation rates in faint X-ray active galactic nuclei

Authors: Nischal Acharya, Silvia Bonoli, Mara Salvato, Ariana Cortesi, M. Rosa González Delgado, Ivan Ezequiel Lopez, Isabel Marquez, Ginés Martínez-Solaeche, Abdurro'uf, David Alexander, Marcella Brusa, Jonás Chaves-Montero, Juan Antonio Fernández Ontiveros, Brivael Laloux, Andrea Lapi, George Mountrichas, Cristina Ramos Almeida, Julio Esteban Rodríguez Martín, Francesco Shankar, Roberto Soria, M. José Vilchez, Raul Abramo, Jailson Alcaniz, Narciso Benitez, Saulo Carneiro , et al. (13 additional authors not shown)

Abstract: 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) galax… ▽ More 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. △ Less

Submitted 2 July, 2024; v1 submitted 9 May, 2024; originally announced May 2024.

Comments: 21 pages, 20 figures, Accepted for Publication in A&A

IRX-CIGALE: a tailored module for Low-Luminosity AGN

Authors: I. E. López, G. Yang, G. Mountrichas, M. Brusa, D. M. Alexander, R. D. Baldi, E. Bertola, S. Bonoli, A. Comastri, F. Shankar, N. Acharya, A. V. Alonso Tetilla, A. Lapi, B. Laloux, X. López López, I. Muñoz Rodríguez, B. Musiimenta, N. Osorio Clavijo, L. Sala, D. Sengupta

Abstract: The spectral energy distribution (SED) of low-luminosity active galactic nuclei (LLAGN) presents unique challenges due to their comparable radiation output to their host galaxies and complex accretion dynamics. We introduce a novel module within the CIGALE framework specifically designed for SED fitting of LLAGN, incorporating both empirical relationships like $L_\mathrm{X}$--$L_\mathrm{12μみゅーm}$ and… ▽ More The spectral energy distribution (SED) of low-luminosity active galactic nuclei (LLAGN) presents unique challenges due to their comparable radiation output to their host galaxies and complex accretion dynamics. We introduce a novel module within the CIGALE framework specifically designed for SED fitting of LLAGN, incorporating both empirical relationships like $L_\mathrm{X}$--$L_\mathrm{12μみゅーm}$ and physically-based accretion models such as advection-dominated accretion flows (ADAFs) and truncated accretion disks. This allows for more accurate depiction of LLAGN central emissions. Using this module, we analyzed a set of 52 X-ray-detected local galaxies, primarily LINERs and Seyferts, and compared its performance to higher-luminosity AGN from the COSMOS and SDSS datasets. Our results show that the module adeptly estimates bolometric luminosities with high precision, despite significant galaxy contamination. It also introduces a versatile X-ray bolometric correction formula covering a vast range of luminosities. Further, our study explored the $αあるふぁ_\mathrm{ox}$ index, which measures the UV to X-ray emission slope, showing that unlike quasars, LLAGN display either stable or only slightly varying $αあるふぁ_\mathrm{ox}$ values, indicating differing accretion and photon production processes in the low luminosity regime. Additionally, we observed a significant drop of 1.4 dex in specific star formation rates when moving from whole galaxies to a central 9-arcsecond aperture in LLAGN, suggesting potential feedback mechanisms at play. Overall, our findings underscore the importance of a multiwavelength approach in AGN studies, highlighting distinct behaviors of LLAGN compared to quasars, thus enhancing our understanding of LLAGN and providing a framework for future comprehensive AGN population studies. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: Submitted to A&A

The Rise of Faint, Red AGN at $z>4$: A Sample of Little Red Dots in the JWST Extragalactic Legacy Fields

Authors: Dale D. Kocevski, Steven L. Finkelstein, Guillermo Barro, Anthony J. Taylor, Antonello Calabrò, Brivael Laloux, Johannes Buchner, Jonathan R. Trump, Gene C. K. Leung, Guang Yang, Mark Dickinson, Pablo G. Pérez-González, Fabio Pacucci, Kohei Inayoshi, Rachel S. Somerville, Elizabeth J. McGrath, Hollis B. Akins, Micaela B. Bagley, Laura Bisigello, Rebecca A. A. Bowler, Adam Carnall, Caitlin M. Casey, Yingjie Cheng, Nikko J. Cleri, Luca Costantin , et al. (32 additional authors not shown)

Abstract: We present a sample of 341 "little red dots" (LRDs) spanning the redshift range $z\sim2-11$ using data from the CEERS, PRIMER, JADES, UNCOVER and NGDEEP surveys. These sources are likely heavily-reddened AGN that trace a previously-hidden phase of dust-obscured black hole growth in the early Universe. Unlike past use of color indices to identify LRDs, we employ continuum slope fitting using shifti… ▽ More We present a sample of 341 "little red dots" (LRDs) spanning the redshift range $z\sim2-11$ using data from the CEERS, PRIMER, JADES, UNCOVER and NGDEEP surveys. These sources are likely heavily-reddened AGN that trace a previously-hidden phase of dust-obscured black hole growth in the early Universe. Unlike past use of color indices to identify LRDs, we employ continuum slope fitting using shifting bandpasses to sample the same rest-frame emission blueward and redward of the Balmer break. This approach allows us to identify LRDs over a wider redshift range and is less susceptible to contamination from galaxies with strong breaks that otherwise lack a rising red continuum. The redshift distribution of our sample increases at $z<8$ and then undergoes a rapid decline at $z\sim4.5$, which may tie the emergence, and obscuration, of these sources to the inside-out growth that galaxies experience during this epoch. We find that LRDs are 2-3 dex more numerous than bright quasars at $z\sim5-7$, but their number density is only 0.6-1 dex higher than X-ray and UV selected AGN at these redshifts. Within our sample, we have identified the first X-ray detected LRDs at $z=3.1$ and $z=4.66$. An X-ray spectral analysis confirms that these AGN are moderately obscured with $\log\,(N_{\rm H}/{\rm cm}^{2}$) of $23.3^{+0.4}_{-1.3}$ and $22.72^{+0.13}_{-0.16}$. Our analysis reveals that reddened AGN emission dominates their rest-optical light, while the rest-UV originates from their host galaxies. We also present NIRSpec follow-up spectroscopy of 17 LRDs that show broad emission lines consistent with AGN activity. The confirmed AGN fraction of our sample is $71\%$ for sources with F444W$<26.5$. In addition, we find three LRDs with narrow blue-shifted Balmer absorption features in their spectra, suggesting an outflow of high-density, low ionization gas from near the central engine of these faint, red AGN. △ Less

Submitted 19 April, 2024; v1 submitted 4 April, 2024; originally announced April 2024.

Comments: 23 pages, 17 figures, submitted to ApJ

Accretion properties of X-ray AGN: Evidence for radiation-regulated obscuration with redshift-dependent host galaxy contribution

Authors: Brivael Laloux, Antonis Georgakakis, David M. Alexander, Johannes Buchner, Carolina Andonie, Nischal Acharya, James Aird, Alba V. Alonso-Tetilla, Angela Bongiorno, Ryan C. Hickox, Andrea Lapi, Blessing Musiimenta, Cristina Ramos Almeida, Carolin Vellforth, Francesco Shankar

Abstract: 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 C… ▽ More 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 $λらむだ\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 $λらむだ$. 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 $λらむだ$ compared to obscured SARD for all redshift intervals, (3) the obscured AGN fraction declines sharply at $\logλらむだ_{\rm break} \sim-2$ for $z <0.5$, but shifts to higher $λらむだ$ values with increasing redshift, (4) the incidence of AGN within the theoretically unstable blow-out region of the $λらむだ-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. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: Submitted to MNRAS

Ionised AGN outflows in the Goldfish galaxy -- The illuminating and interacting red quasar eFEDSJ091157.4+014327 at z $\sim$ 0.6

Authors: Blessing Musiimenta, Giovanna Speranza, Tanya Urrutia, Marcella Brusa, Cristina Ramos Almeida, Michele Perna, Ivàn Ezequiel López, David M. Alexander, Brivael Laloux, Francesco Shankar, Andrea Lapi, Mara Salvato, Yoshiki Toba, Carolina Andonie, Ivàn Munoz Rodríguez

Abstract: 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 outf… ▽ More 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. △ Less

Submitted 6 May, 2024; v1 submitted 30 January, 2024; originally announced January 2024.

Comments: Version after addressing referee comments, accepted in Astronomy and Astrophysics(A&A), 17 pages, 12 figures, 1 table

Journal ref: A&A 687, A111 (2024)

AGN X-ray luminosity function and absorption function in the Early Universe ($3\leq z \leq 6$)

Authors: E. Pouliasis, A. Ruiz, I. Georgantopoulos, F. Vito, R. Gilli, C. Vignali, Y. Ueda, E. Koulouridis, M. Akiyama, S. Marchesi, B. Laloux, T. Nagao, S. Paltani, M. Pierre, Y. Toba, M. Habouzit, B. Vijarnwannaluk, C. Garrel

Abstract: The XLF of AGN offers a robust tool to study the evolution and the growth of SMBHs over cosmic time. Owing to the limited area probed by X-ray surveys, optical surveys are routinely used to probe the accretion in the high redshift Universe $z\geq 3$. However, optical surveys may be incomplete because they are strongly affected by dust redenning. In this work, we derive the XLF and its evolution at… ▽ More The XLF of AGN offers a robust tool to study the evolution and the growth of SMBHs over cosmic time. Owing to the limited area probed by X-ray surveys, optical surveys are routinely used to probe the accretion in the high redshift Universe $z\geq 3$. However, optical surveys may be incomplete because they are strongly affected by dust redenning. In this work, we derive the XLF and its evolution at high redshifts using a large sample of AGNs selected in different fields with various areas and depths covering a wide range of luminosities. Additionally, we put the tightest yet constraints on the absorption function in this redshift regime. In particular, we use more than 600 soft X-ray selected high-z sources in the Chandra Deep fields, the Chandra COSMOS Legacy survey and the XMM-XXL northern field. We derive the X-ray spectral properties for all sources via spectral fitting, using a consistent technique and model. For modeling the parametric form of the XLF and the absorption function, we use a Bayesian methodology allowing us to correctly propagate the uncertainties for the observed X-ray properties of our sources and also the absorption effects. The evolution of XLF is in agreement with a pure density evolution model similar to what is witnessed at optical wavelengths, although a luminosity dependent density evolution model cannot be securely ruled out. A large fraction ($60\%)$ of our sources are absorbed by column densities of $\rm N_H \geq 10^{23} cm^{-2} $, while $17$\% of the sources are CTK. Our results favor a scenario where both the ISM of the host and the AGN torus contribute to the obscuration. The derived BHAD is in agreement with the simulations, if one takes into account that the X-ray AGN are hosted by massive galaxies, while it differs from the one derived using JWST data. The latter could be due to the differences in the AGN and host-galaxy properties. △ Less

Submitted 24 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A. 25 pages, 15 figures (+4 in Appendix), 2 tables (+1 in Appendix)

Journal ref: A&A 685, A97 (2024)

Probing the roles of orientation and multi-scale gas distributions in shaping the obscuration of Active Galactic Nuclei through cosmic time

Authors: Alba V. Alonso-Tetilla, Francesco Shankar, Fabio Fontanot, Nicola Menci, Milena Valentini, Johannes Buchner, Brivael Laloux, Andrea Lapi, Annagrazia Puglisi, David M. Alexander, Viola Allevato, Carolina Andonie, Silvia Bonoli, Michaela Hirschmann, Ivan E. Lopez, Sandra I. Raimundo, Cristina Ramos Almeida

Abstract: 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 Mod… ▽ More 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. △ Less

Submitted 24 October, 2023; originally announced October 2023.

Comments: MNRAS, accepted, 19 pages, 15 figures, 3 appendices

Obscuration beyond the nucleus: infrared quasars can be buried in extreme compact starbursts

Authors: Carolina Andonie, David M. Alexander, Claire Greenwell, Annagrazia Puglisi, Brivael Laloux, Alba V. Alonso-Tetilla, Gabriela Calistro Rivera, Chris Harrison, Ryan C. Hickox, Melanie Kaasinen, Andrea Lapi, Iván E. López, Grayson Petter, Cristina Ramos Almeida, David J. Rosario, Francesco Shankar, Carolin Villforth

Abstract: 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 in… ▽ More 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. △ Less

Submitted 10 October, 2023; v1 submitted 3 October, 2023; originally announced October 2023.

Comments: Accepted for publication in MNRAS Letters

The miniJPAS survey: AGN & host galaxy co-evolution of X-ray selected sources

Authors: I. E. López, M. Brusa, S. Bonoli, F. Shankar, N. Acharya, B. Laloux, K. Dolag, A. Georgakakis, A. Lapi, C. Ramos Almeida, M. Salvato, J. Chaves-Montero, P. Coelho, L. A. Díaz-García, J. A. Fernández-Ontiveros, A. Hernán-Caballero, R. M. González Delgado, I. Marquez, M. Pović, R. Soria, C. Queiroz, P. T. Rahna, R. Abramo, J. Alcaniz, N. Benitez , et al. (13 additional authors not shown)

Abstract: Studies indicate strong evidence of a scaling relation in the local Universe between the supermassive black hole mass ($M_\rm{BH}$) and the stellar mass of their host galaxies ($M_\star$). They even show similar histories across cosmic times of their differential terms: star formation rate (SFR) and black hole accretion rate (BHAR). However, a clear picture of this coevolution is far from being un… ▽ More Studies indicate strong evidence of a scaling relation in the local Universe between the supermassive black hole mass ($M_\rm{BH}$) and the stellar mass of their host galaxies ($M_\star$). They even show similar histories across cosmic times of their differential terms: star formation rate (SFR) and black hole accretion rate (BHAR). However, a clear picture of this coevolution is far from being understood. We select an X-ray sample of active galactic nuclei (AGN) up to $z=2.5$ in the miniJPAS footprint. Their X-ray to infrared spectral energy distributions (SEDs) have been modeled with CIGALE, constraining the emission to 68 bands. For a final sample of 308 galaxies, we derive their physical properties (e.g., $M_\star$, $\rm{SFR}$, $\rm{SFH}$, and $L_\rm{AGN}$). We also fit their optical spectra for a subsample of 113 sources to estimate the $M_\rm{BH}$. We calculate the BHAR depending on two radiative efficiency regimes. We find that the Eddington ratios ($λらむだ$) and its popular proxy ($L_\rm{X}$/$M_\star$) have 0.6 dex of difference, and a KS-test indicates that they come from different distributions. Our sources exhibit a considerable scatter on the $M_\rm{BH}$-$M_\star$ relation, which can explain the difference between $λらむだ$ and its proxy. We also model three evolution scenarios to recover the integral properties at $z=0$. Using the SFR and BHAR, we show a notable diminution in the scattering between $M_\rm{BH}$-$M_\star$. For the last scenario, we consider the SFH and a simple energy budget for the AGN accretion, obtaining a relation similar to the local Universe. Our study covers $\sim 1$ deg$^2$ in the sky and is sensitive to biases in luminosity. Nevertheless, we show that, for bright sources, the link between SFR and BHAR, and their decoupling based on an energy limit is the key that leads to the local $M_\rm{BH}$-$M_\star$ scaling relation. △ Less

Submitted 2 February, 2023; originally announced February 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 672, A137 (2023)

The demographics of obscured AGN from X-ray spectroscopy guided by multiwavelength information

Authors: Brivael Laloux, Antonis Georgakakis, Carolina Andonie, David M. Alexander, Angel Ruiz, David J. Rosario, James Aird, Johannes Buchner, Francisco J. Carrera, Andrea Lapi, Cristina Ramos Almeida, Mara Salvato, Francesco Shankar

Abstract: 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 pri… ▽ More 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μみゅーm$ spectral energy distributions with galaxy and AGN templates to determine the mid-infrared ($6μみゅーm$) luminosity of the AGN component. Empirical correlations between X-ray and $6μみゅー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. △ Less

Submitted 7 November, 2022; originally announced November 2022.

Comments: Accepted for publication at MNRAS, 21 pages, 20 figures, 6 tables

A panchromatic view of infrared quasars: excess star formation and radio emission in the most heavily obscured systems

Authors: Carolina Andonie, David M. Alexander, David Rosario, Brivael Laloux, Antonis Georgakakis, Leah K. Morabito, Carolin Villforth, Mathilda Avirett-Mackenzie, Gabriela Calistro Rivera, Agnese Del Moro, Sotiria Fotopoulou, Chris Harrison, Andrea Lapi, James Petley, Grayson Petter, Francesco Shankar

Abstract: 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$,… ▽ More 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\: μみゅー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. △ Less

Submitted 27 September, 2022; originally announced September 2022.

Comments: Accepted for publication in MNRAS. 22 pages, 15 figures, 4 tables

Towards modelling Ghostly DLAs

Authors: Brivael Laloux, Patrick Petitjean, Camille Noûs

Abstract: We use simple models of the spatial structure of the quasar broad line region (BLR) to investigate the properties of so-called ghostly damped Lyman-αあるふぁ (DLA) systems detected in SDSS data. These absorbers are characterized by the presence of strong metal lines but no Hi Lyman-αあるふぁ trough is seen in the quasar spectrum indicating that, although the region emitting the quasar continuum is covered by an a… ▽ More We use simple models of the spatial structure of the quasar broad line region (BLR) to investigate the properties of so-called ghostly damped Lyman-αあるふぁ (DLA) systems detected in SDSS data. These absorbers are characterized by the presence of strong metal lines but no Hi Lyman-αあるふぁ trough is seen in the quasar spectrum indicating that, although the region emitting the quasar continuum is covered by an absorbing cloud, the BLR is only partially covered. One of the models has a spherical geometry, another one is the combination of two wind flows whereas the third model is a Keplerian disk. The models can reproduce the typical shape of the quasar Lyman-αあるふぁ emission and different ghostly configurations. We show that the DLA Hi column density can be recovered precisely independently of the BLR model used. The size of the absorbing cloud and its distance to the centre of the AGN are correlated. However it may be possible to disentangle the two using an independent estimate of the radius from the determination of the particle density. Comparison of the model outputs with SDSS data shows that the wind and disk models are more versatile than the spherical one and can be more easily adapted to the observations. For all the systems we derive logN(Hi)(cm^{-2})>20.5. With higher quality data it may be possible to distinguish between the models. △ Less

Submitted 20 January, 2021; originally announced January 2021.

Comments: Accepted for publication in MNRAS