Showing 1–50 of 282 results for author: Gilli, R

Fast SMBH growth in the SPT2349--56 protocluster at $z=4.3$

Authors: Fabio Vito, William N. Brandt, Andrea Comastri, Roberto Gilli, Rob J. Ivison, Giorgio Lanzuisi, Bret D. Lehmer, Ivan E. Lopez, Paolo Tozzi, Cristian Vignali

Abstract: Protoclusters at $z>2$ are gas-rich regions characterized by high star-formation activity. The same physical properties that enhance star formation in protoclusters are also thought to boost the growth of SMBHs. We aim to test this scenario by probing the AGN content of SPT2349-56, a massive, gas-rich, and highly star-forming protocluster core at $z=4.3$ discovered as an overdensity of DSFGs, via… ▽ More Protoclusters at $z>2$ are gas-rich regions characterized by high star-formation activity. The same physical properties that enhance star formation in protoclusters are also thought to boost the growth of SMBHs. We aim to test this scenario by probing the AGN content of SPT2349-56, a massive, gas-rich, and highly star-forming protocluster core at $z=4.3$ discovered as an overdensity of DSFGs, via Chandra (200 ks) observations, and comparing the results with the field environment. We detected two protocluster members, corresponding to an AGN fraction among DSFGs of $\approx10\%$. This value is consistent with other protoclusters at $z=2-4$, but higher than the AGN incidence among DSFGs in the field environment. Both AGN are heavily obscured sources and hosted in star-forming galaxies with $\approx3\times10^{10}\,\mathrm{M_\odot}$ stellar masses. We estimate that the ISM in the host galaxies can contribute significantly to the nuclear obscuration. One of the two AGN is highly luminous ($L_X=2\times10^{45}\,\mathrm{erg\,s^{-1}}$) and Compton-thick ($N_H=2\times10^{24}\,\mathrm{cm^{-2}}$), and likely powered by a $M_{BH}>6\times10^8\,\mathrm{M_\odot}$ SMBH. Its high accretion rate suggests that it is in the phase of efficient growth required to explain the presence of extremely massive SMBHs in the centers of local galaxy clusters. Considering SPT2349-56 and DRC, a similar protocuster at $z=4$, we find that gas-rich protocluster cores at $z\approx4$ enhance the triggering of luminous (log$\frac{L_X}{\mathrm{erg\,s^{-1}}}=45-46$) AGN by 3-5 orders of magnitude with respect to the field environment. Our results indicate that gas-rich protoclusters at high redshift boost the growth of SMBHs, which will likely impact the subsequent evolution of the structures, and thus represent key science targets to obtain a complete understanding of the relation between environment and galaxy evolution. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: Submitted to A&A, referee comments included

The XMM-Newton and NuSTAR view of IRASF11119+3257. I Detection of multiple UFO components and a very cold corona

Authors: G. Lanzuisi, G. Matzeu, P. Baldini, E. Bertola, A. Comastri, F. Tombesi, A. Luminari, V. Braito, J. Reeves, G. Chartas, S. Bianchi, M. Brusa, G. Cresci, E. Nardini, E. Piconcelli, L. Zappacosta, R. Serafinelli, M. Gaspari, R. Gilli, M. Cappi, M. Dadina, M. Perna, C. Vignali, S. Veilleux

Abstract: IRASF11119 is an ultra-luminous IR galaxy with post-merger morphology, hosting a type-1 QSO at z=0.189. Its 2013 Suzaku spectrum shows a prominent Ultra Fast Outflow (UFO) absorption feature (v_out~0.25c). In 2021, we obtained the first XMM-Newton long look of the target, coordinated with a simultaneous NuSTAR observation. The new high-quality data allow us to detect at P>99.8% c.l. multiple absor… ▽ More IRASF11119 is an ultra-luminous IR galaxy with post-merger morphology, hosting a type-1 QSO at z=0.189. Its 2013 Suzaku spectrum shows a prominent Ultra Fast Outflow (UFO) absorption feature (v_out~0.25c). In 2021, we obtained the first XMM-Newton long look of the target, coordinated with a simultaneous NuSTAR observation. The new high-quality data allow us to detect at P>99.8% c.l. multiple absorption features associated with the known UFO. Furthermore, an emission plus absorption feature at 1.1-1.3 keV reveals the presence of a blueshifted P-Cygni profile in the soft band. We associate the hard band features with blends of FeXXV and FeXXVI He$α$-Ly$α$ and He$β$-Ly$β$ line pairs and infer a large column (N$_H$~$10^{24}$ cm$^{-2}$) of highly ionized (log$ξ$~5) gas outflowing at v_out=0.27c. The 1 keV feature can be associated with a blend of Fe and Ne transitions, produced by a lower column (N$_H$~$10^{21}$ cm$^{-2}$) and ionization (log$ξ$~2.6) gas component outflowing at the same speed. Using a radiative-transfer disk wind model to fit the highly ionized UFO, we derive a large mass outflow rate, comparable with the mass accretion rate (M$_{out}$=4.25 M$_{Sun}$/yr, ~1.6 M$_{acc}$), and kinetic energy and momentum flux among the highest reported in the literature. We measure an extremely low high-energy cut-off (E$_c$~25 keV). Several other cases in the literature suggest that a steep X-ray continuum may be related to the formation of powerful winds. The lack of a significant momentum boost between the nuclear UFO and the different phases of the large-scale outflow, observed in IRASF11119 and in a growing number of sources with powerful UFOs, can be explained by (i) a momentum-driven expansion, (ii) an inefficient coupling of the UFO with the host ISM, or (iii) by repeated energy-driven expansion episodes with low duty-cycle, that average out on long time-scales. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 21 pages, 13 figures, accepted for publication in A&A

A quasar-galaxy merger at $z\sim 6.2$: rapid host growth via accretion of two massive satellite galaxies

Authors: Roberto Decarli, Federica Loiacono, Emanuele Paolo Farina, Massimo Dotti, Alessandro Lupi, Romain A. Meyer, Marco Mignoli, Antonio Pensabene, Michael A. Strauss, Bram Venemans, Jinyi Yang, Fabian Walter, Julien Wolf, Eduardo Bañados, Laura Blecha, Sarah Bosman, Chris L. Carilli, Andrea Comastri, Thomas Connor, Tiago Costa, Anna-Christina Eilers, Xiaohui Fan, Roberto Gilli, Hyunsung D. Jun, Weizhe Liu , et al. (16 additional authors not shown)

Abstract: We present JWST/NIRSpec Integral Field Spectroscopy in the rest-frame optical bands of the system PJ308-21, a quasar at $z=6.2342$ caught as its host galaxy interacts with companion galaxies. We detect spatially extended emission of several emission lines (H$α$, H$β$, [OIII], [NII], [SII], HeII), which we use to study the properties of the ionized phase of the interstellar medium: the source and h… ▽ More We present JWST/NIRSpec Integral Field Spectroscopy in the rest-frame optical bands of the system PJ308-21, a quasar at $z=6.2342$ caught as its host galaxy interacts with companion galaxies. We detect spatially extended emission of several emission lines (H$α$, H$β$, [OIII], [NII], [SII], HeII), which we use to study the properties of the ionized phase of the interstellar medium: the source and hardness of the photoionizing radiation field, metallicity, dust reddening, electron density and temperature, and star formation. We also marginally detect continuum starlight emission associated with the companion sources. We find that at least two independent satellite galaxies are part of the system. While the quasar host appears highly enriched and obscured, with AGN-like photoionization conditions, the western companion shows minimal dust extinction, low metallicity ($Z\sim0.4$ Z$_\odot$), and star-formation driven photoionization. The eastern companion shows higher extinction and metallicity ($Z\sim0.8$ Z$_\odot$) compared to the western companion, and it is at least partially photoionized by the nearby quasar. We do not find any indication of AGN in the companion sources. Our study shows that while the quasar host galaxy is already very massive ($M_{\rm dyn}>10^{11}$ M$_\odot$), it is still rapidly building up by accreting two relatively massive ($M_{\rm star}\sim 10^{10}$ M$_\odot$) companion sources. This dataset showcases the power of JWST in exposing the build-up of massive galaxies in the first Gyr of the Universe. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 15 pages, 16 figures. Accepted for publication in A&A

JWST meets Chandra: a large population of Compton thick, feedback-free, and X-ray weak AGN, with a sprinkle of SNe

Authors: Roberto Maiolino, Guido Risaliti, Matilde Signorini, Bartolomeo Trefoloni, Ignas Juodzbalis, Jan Scholtz, Hannah Uebler, Francesco D'Eugenio, Stefano Carniani, Andy Fabian, Xihan Ji, Giovanni Mazzolari, Elena Bertola, Marcella Brusa, Andrew J. Bunker, Stephane Charlot, Andrea Comastri, Giovanni Cresci, Christa Noel DeCoursey, Eiichi Egami, Fabrizio Fiore, Roberto Gilli, Michele Perna, Sandro Tacchella, Giacomo Venturi

Abstract: We investigate the X-ray properties of a large sample of 71 broad line and narrow line AGN at 2<z<11 discovered by JWST in the GOODS fields, which have the deepest Chandra observations ever obtained. Despite the widespread presence of AGN signatures in their rest-optical and -UV spectra, the vast majority of them is X-ray undetected. The stacked X-ray data of the non-detected sources also results… ▽ More We investigate the X-ray properties of a large sample of 71 broad line and narrow line AGN at 2<z<11 discovered by JWST in the GOODS fields, which have the deepest Chandra observations ever obtained. Despite the widespread presence of AGN signatures in their rest-optical and -UV spectra, the vast majority of them is X-ray undetected. The stacked X-ray data of the non-detected sources also results in a non-detection. The upper limit on the X-ray emission for many of these AGN is one or even two orders of magnitude lower than expected from a standard AGN SED. Heavy X-ray absorption by clouds with large (Compton thick) column density and low dust content, such as the Broad Line Region (BLR) clouds, can explain the X-ray weakness. In this scenario the BLR covering factor should be much larger than in low-z AGN or luminous quasar; this is supported by the larger equivalent width of the broad component of Halpha in JWST-selected AGN. We also find that the JWST-discovered AGN lack the prominent, fast outflows characterizing low-z AGN and luminous quasars, suggesting that, in JWST-selected AGN, dense gas lingers in the nuclear region, resulting in large covering factors. We also note that a large fraction of JWST-selected AGN match the definition of NLSy1, typically characterized by a steep X-ray spectrum, and this can further contribute to their observed weakness at high-z. Finally, we discuss that the broad Balmer lines used to identify type 1 AGN cannot be ascribed to Very Massive Stars, Tidal Disruption Events, or Supernovae, although we show that a minority of the faintest broad lines could potentially be associated with the echo of superluminous SNe or TDE. Scenarios in which the broad lines are ascribed to galactic outflows are also untenable. We emphasize that confirming any of the scenarios discussed above will require X-ray missions more sensitive than Chandra. (abridged) △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 21 pages, 15 figures, 2 tables, submitted to MNRAS

The neutrino background from non-jetted active galactic nuclei

Authors: P. Padovani, R. Gilli, E. Resconi, C. Bellenghi, F. Henningsen

Abstract: Aims. We calculate the contribution to the neutrino background from the non-jetted active galactic nuclei (AGN) population following the recent IceCube association of TeV neutrinos with NGC 1068. Methods. We exploit our robust knowledge of the AGN X-ray luminosity function and evolution and convert it to the neutrino band by using NGC 1068 as a benchmark and a theoretically motivated neutrino spec… ▽ More Aims. We calculate the contribution to the neutrino background from the non-jetted active galactic nuclei (AGN) population following the recent IceCube association of TeV neutrinos with NGC 1068. Methods. We exploit our robust knowledge of the AGN X-ray luminosity function and evolution and convert it to the neutrino band by using NGC 1068 as a benchmark and a theoretically motivated neutrino spectrum. Results. The resulting neutrino background up to redshift 5 does not violate either the IceCube diffuse flux or the upper bounds for non-jetted AGN, although barely so. This is consistent with a scenario where the latter class makes a substantial contribution mostly below 1 PeV, while jetted AGN, i.e. blazars, dominate above this energy, in intriguing agreement with the dip in the neutrino data at ~ 300 TeV. More and better IceCube data on Seyfert galaxies will allow us to constrain the fraction of neutrino emitters among non-jetted AGN. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: 4 pages, 1 figure, accepted for publication in Astronomy & Astrophysics Letters

Journal ref: A&A 684, L21 (2024)

The ALPINE-ALMA [CII] Survey: Dust emission effective radius up to 3 kpc in the Early Universe

Authors: F. Pozzi, F. Calura, Q. D'Amato, M. Gavarente, M. Bethermin, M. Boquien, V. Casasola, A. Cimatti, R. Cochrane, M. Dessauges-Zavadsky, A. Enia, F. Esposito, A. L. Faisst, R. Gilli, M. Ginolfi, R. Gobat, C. Gruppioni, C. C. Hayward, E. Ibar, A. M. Koekemoer, B. C. Lemaux, G. E. Magdis, J. Molina, M. Talia, L. Vallini , et al. (2 additional authors not shown)

Abstract: Measurements of the size of dust continuum emission are an important tool for constraining the spatial extent of star formation and hence the build-up of stellar mass. Compact dust emission has generally been observed at Cosmic Noon (z~2-3). However, at earlier epochs, toward the end of the Reionization (z~4-6), only the sizes of a handful of IR-bright galaxies have been measured. In this work, we… ▽ More Measurements of the size of dust continuum emission are an important tool for constraining the spatial extent of star formation and hence the build-up of stellar mass. Compact dust emission has generally been observed at Cosmic Noon (z~2-3). However, at earlier epochs, toward the end of the Reionization (z~4-6), only the sizes of a handful of IR-bright galaxies have been measured. In this work, we derive the dust emission sizes of main-sequence galaxies at z~5 from the ALPINE survey. We measure the dust effective radius r_e,FIR in the uv-plane in Band 7 of ALMA for seven ALPINE galaxies with resolved emission and we compare it with rest-frame UV and [CII]158$μ$m measurements. We study the r_e,FIR-L_IR scaling relation by considering our dust size measurements and all the data in literature at z~4-6. Finally, we compare our size measurements with predictions from simulations. The dust emission in the selected ALPINE galaxies is rather extended (r_e,FIR~1.5-3 kpc), similar to [CII]158 um but a factor of ~2 larger than the rest-frame UV emission. Putting together all the measurements at z~5, spanning 2 decades in luminosity from L_IR ~ 10^11 L_sun to L_IR ~ 10^13 L_sun, the data highlight a steeply increasing trend of the r_e,FIR-L_IR relation at L_IR< 10^12 L_sun, followed by a downturn and a decreasing trend at brighter luminosities. Finally, simulations that extend up to the stellar masses of the ALPINE galaxies considered in the present work predict a sub-set of galaxies (~25% at 10^10 M_sun < M_star < 10^11 M_sun) with sizes as large as those measured. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Comments: 8 pages, Accepted by A&A

Journal ref: A&A 686, A187 (2024)

Supermassive Black Hole Winds in X-rays -- SUBWAYS. III. A population study on ultra-fast outflows

Authors: V. E. Gianolli, S. Bianchi, P-O Petrucci, M. Brusa, G. Chartas, G. Lanzuisi, G. A. Matzeu, M. Parra, F. Ursini, E. Behar, M. Bischetti, A. Comastri, E. Costantini, G. Cresci, M. Dadina, B. De Marco, A. De Rosa, F. Fiore, M. Gaspari, R. Gilli, M. Giustini, M. Guainazzi, A. R. King, S. Kraemer, G. Kriss , et al. (22 additional authors not shown)

Abstract: The detection of blue-shifted absorption lines likely associated with ionized Iron K-shell transitions in the X-ray spectra of many Active Galactic Nuclei (AGN) suggests the presence of a highly ionized gas outflowing with mildly relativistic velocities (0.03c-0.6c), named Ultra-Fast Outflow (UFO). Within the SUBWAYS project we characterized these winds starting from a sample of 22 radio-quiet qua… ▽ More The detection of blue-shifted absorption lines likely associated with ionized Iron K-shell transitions in the X-ray spectra of many Active Galactic Nuclei (AGN) suggests the presence of a highly ionized gas outflowing with mildly relativistic velocities (0.03c-0.6c), named Ultra-Fast Outflow (UFO). Within the SUBWAYS project we characterized these winds starting from a sample of 22 radio-quiet quasars at 0.1 < z < 0.4, and compared the results with similar studies in the literature on samples of 42 local radio-quiet Seyfert galaxies and 14 high redshift radio-quiet quasars. The scope of our work is a statistical study of UFO parameters and incidence, considering key physical properties of the sources, e.g. supermassive black hole (SMBH) mass, bolometric luminosity, accretion rates and Spectral Energy Distribution, with the aim of gaining new insights into the UFO launching mechanisms. We find indications that highly luminous AGN with steeper X-ray/UV ratio, are more likely to host UFO. The presence of UFO is not significantly related to any other AGN property in our sample. These findings suggest that the UFO phenomenon may be transient. Focusing on AGN with UFO, other important results are: (1) faster UFO have larger ionization parameters and column densities; (2) X-ray radiation plays a more crucial role in driving highly ionized winds compared to UV; (3) the correlation between outflow velocity and luminosity is significantly flatter than what expected for radiatively driven winds; (4) more massive BH experience higher wind mass-losses, suppressing accretion of matter onto the BH; (5) the UFO launching radius is positively correlated with the Eddington ratio. Furthermore, our analysis suggest the involvement of multiple launching mechanisms, including radiation pressure and magneto-hydrodynamic processes, rather than pointing to a single, universally applicable mechanism. △ Less

Submitted 11 April, 2024; v1 submitted 14 March, 2024; originally announced March 2024.

Comments: 52 pages, 154 figures and 7 tables. Accepted in Astronomy & Astrophysics, in press

COSMOS2020: Investigating the AGN-obscured accretion phase at $z\sim 1$ via [NeV] selection

Authors: L. Barchiesi, C. Vignali, F. Pozzi, R. Gilli, M. Mignoli, C. Gruppioni, A. Lapi, S. Marchesi, F. Ricci, C. M. Urry

Abstract: We investigated the properties of 94 [NeV]3426AA-selected type 2 AGN in COSMOS at z=0.6-1.2, performing optical-to-far-infrared spectral energy distribution fitting. In addition, we analyze the X-ray spectra of the X-ray-detected sources to obtain reliable values of the AGN obscuration and intrinsic luminosity. We found that more than two-thirds of our sample is composed of very obscured sources… ▽ More We investigated the properties of 94 [NeV]3426AA-selected type 2 AGN in COSMOS at z=0.6-1.2, performing optical-to-far-infrared spectral energy distribution fitting. In addition, we analyze the X-ray spectra of the X-ray-detected sources to obtain reliable values of the AGN obscuration and intrinsic luminosity. We found that more than two-thirds of our sample is composed of very obscured sources, with about 20% of the sources being candidate CT-AGN and half being AGNs in a strong phase of accretion. With respect to non-active galaxies, we find a higher fraction of sources within the main sequence and little evidence for AGNs quenching the SF. The comparison with the prediction from the in situ co-evolution model suggests that [NeV] is an effective tool for selecting galaxies in the obscured growth phase of the BH-galaxy co-evolution paradigm. We find that the "quenching phase" is still to come for most of the sample and only few galaxies show evidence of quenched SF activity. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: Accepted in A&A, 16 pages, 14 figures

Winds of change: the nuclear and galaxy-scale outflows and the X-ray variability of 2MASS 0918+2117

Authors: P. Baldini, G. Lanzuisi, M. Brusa, A. Merloni, K. Gkimisi, M. Perna, I. E. Lopez, E. Bertola, Z. Igo, S. Waddell, B. Musiimenta, C. Aydar, R. Arcodia, G. A. Matzeu, A. Luminari, J. Buchner, C. Vignali, M. Dadina, A. Comastri, G. Cresci, S. Marchesi, R. Gilli, F. Tombesi, R. Serafinelli

Abstract: Powerful outflows from active galactic nuclei (AGN) can significantly impact the gas reservoirs of their host galaxies. However, it is still unclear how these outflows can propagate from the very central regions of galaxies to their outskirts, and whether nuclear winds can be driven by and/or be responsible for drastic spectral transitions. In this work we test feedback propagation models on the c… ▽ More Powerful outflows from active galactic nuclei (AGN) can significantly impact the gas reservoirs of their host galaxies. However, it is still unclear how these outflows can propagate from the very central regions of galaxies to their outskirts, and whether nuclear winds can be driven by and/or be responsible for drastic spectral transitions. In this work we test feedback propagation models on the case test of 2MASS 0918+2117 (2M0918), a z=0.149 X-ray variable AGN, which showed tentative evidence for nuclear ultra-fast outflows (UFOs) in a 2005 XMM-Newton observation. We also investigate whether UFOs can be related to the observed X-ray variability. We observed 2M0918 with XMM-Newton and NuSTAR in 2020 to confirm the presence and characterize the UFOs. We perform a kinematic analysis of the 2005 SDSS optical spectrum to reveal and measure the properties of galaxy-scale ionized outflows. Furthermore, we construct 20-year-long lightcurves of observed flux, line-of-sight column density, and intrinsic accretion rate from the spectra of the first 4 SRG/eROSITA all-sky surveys and archival observations from Chandra and XMM-Newton.We significantly detect UFOs with v$\sim$0.16c and galaxy-scale ionized outflows with velocities of $\sim$ 700 km/s. We also find that the drastic X-ray variability (factors >10) can be explained both in terms of variable obscuration and variable intrinsic luminosity.Comparing the energetics of the two outflow phases, 2M0918 is consistent with momentum-driven wind propagation. 2M0918 expands the sample of AGN with both UFOs and ionized gas winds from 5 to 6, and brings the sample of AGN hosting multiscale outflows to 19, contributing to a clearer picture of feedback physics. From the variations in accretion rate, column density, and ionization level of the obscurer, we propose a scenario that connects obscurers, an accretion enhancement, and the emergence of UFOs △ Less

Submitted 26 February, 2024; originally announced February 2024.

Comments: 18 pages, 18 figures

Journal ref: A&A 686, A217 (2024)

A quasar-galaxy merger at $z\sim6.2$: black hole mass and quasar properties from the NIRSpec spectrum

Authors: Federica Loiacono, Roberto Decarli, Marco Mignoli, Emanuele Paolo Farina, Eduardo Bañados, Sarah Bosman, Anna-Christina Eilers, Jan-Torge Schindler, Michael A. Strauss, Marianne Vestergaard, Feige Wang, Laura Blecha, Chris L. Carilli, Andrea Comastri, Thomas Connor, Tiago Costa, Massimo Dotti, Xiaohui Fan, Roberto Gilli, Hyunsung D. Jun, Weizhe Liu, Alessandro Lupi, Madeline A. Marshall, Chiara Mazzucchelli, Romain A. Meyer , et al. (9 additional authors not shown)

Abstract: We present JWST/NIRSpec integral field data of the quasar PJ308-21 at $z=6.2342$. As shown by previous ALMA and HST imaging, the quasar has two companion sources, interacting with the quasar host galaxy. The high-resolution G395H/290LP NIRSpec spectrum covers the $2.87-5.27\ \rm μm$ wavelength range and shows the rest-frame optical emission of the quasar with exquisite quality ($S/N\sim 100-400$ p… ▽ More We present JWST/NIRSpec integral field data of the quasar PJ308-21 at $z=6.2342$. As shown by previous ALMA and HST imaging, the quasar has two companion sources, interacting with the quasar host galaxy. The high-resolution G395H/290LP NIRSpec spectrum covers the $2.87-5.27\ \rm μm$ wavelength range and shows the rest-frame optical emission of the quasar with exquisite quality ($S/N\sim 100-400$ per spectral element). Based on the H$β$ line from the broad line region, we obtain an estimate of the black hole mass $M_{\rm BH,Hβ}\sim 2.7\times 10^{9}\ \rm M_{\odot}$. This value is within a factor $\lesssim 1.5$ of the H$α$-based black hole mass from the same spectrum ($M_{\rm BH, Hα}\sim 1.93\times 10^{9}\ \rm M_{\odot}$) and is consistent with a previous estimate relying on the MgII $λ2799$ ($M_{\rm BH, MgII}\sim 2.65\times 10^{9}\ \rm M_{\odot}$). All these $M_{\rm BH}$ are within the $\sim 0.5$ dex intrinsic scatter of the adopted mass calibrations. The high Eddington ratio of PJ308-21 $λ_{\rm Edd,Hβ}\sim 0.67$ ($λ_{\rm Edd,Hα}\sim 0.96$) is in line with the overall quasar population at $z \gtrsim 6$. The relative strengths of the [OIII], FeII and H$β$ lines are consistent with the empirical "Eigenvector 1" correlations as observed for low redshift quasars. We find evidence for blueshifted [OIII] $λ5007$ emission with a velocity offset $Δv_{\rm [OIII]}=-1922\pm 39$ km s$^{-1}$ from the systemic velocity and a $\rm FWHM([OIII])=2776^{+75}_{-74}$ km s$^{-1}$. This may be the signature of an outflow from the nuclear region, despite the true values of $Δv_{\rm [OIII]}$ and $\rm FWHM([OIII])$ are likely more uncertain due to the blending with H$β$ and FeII lines. Our study demonstrates the unique capabilities of NIRSpec in capturing quasar spectra at cosmic dawn and studying their properties in unprecedented detail. △ Less

Submitted 20 February, 2024; originally announced February 2024.

Comments: 17 pages, 9 figures, accepted for publication in A&A

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

Heavily Obscured AGN detection: a Radio vs X-ray challenge

Authors: Giovanni Mazzolari, Roberto Gilli, Marcella Brusa, Marco Mignoli, Fabio Vito, Isabella Prandoni, Stefano Marchesi, Marco Chiaberge, Giorgio Lanzuisi, Quirino D'Amato, Andrea Comastri, Cristian Vignali, Kazushi Iwasawa, Colin Norman

Abstract: In this work, we study the AGN radio detection effectiveness in the major deep extragalactic surveys, considering different AGN obscuration levels, redshift, and AGN bolometric luminosities. We particularly focus on comparing their radio and X-ray detectability, making predictions for present and future radio surveys. We extrapolate the predictions of AGN population synthesis model of cosmic X-ray… ▽ More In this work, we study the AGN radio detection effectiveness in the major deep extragalactic surveys, considering different AGN obscuration levels, redshift, and AGN bolometric luminosities. We particularly focus on comparing their radio and X-ray detectability, making predictions for present and future radio surveys. We extrapolate the predictions of AGN population synthesis model of cosmic X-ray background (CXB) to the radio band, by deriving the 1.4 GHz luminosity functions of unobscured (i.e. with hydrogen column densities $\log N_{H} <22$), obscured ($22<\log N_{H}<24$) and Compton-thick (CTK, $\log N_{H} >24$) AGN. We then use these functions to forecast the number of detectable AGN based on the area, flux limit, and completeness of a given radio survey and compare it with the AGN number resulting from X-ray predictions. When applied to deep extragalactic fields covered both by radio and X-ray observations, we show that, while X-ray selection is generally more effective in detecting unobscured AGN, the surface density of CTK AGN radio detected is on average $\sim 10$ times larger than the X-ray one, and even greater at high redshifts, considering the current surveys and facilities. Our results suggest that thousands of CTK AGN are already present in current radio catalogs, but most of them escaped any detection in the corresponding X-ray observations. We also present expectations for the number of AGN to be detected by the Square Kilometer Array Observatory (SKAO) in its future deep and wide radio continuum surveys, finding that it will be able to detect more than 2000 AGN at $z>6$ and some tens at $z>10$, more than half of which are expected to be CTK. △ Less

Submitted 31 January, 2024; originally announced February 2024.

Comments: 17 pages, 11 figures, accepted for publication on A&A

Journal ref: A&A 687, A120 (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)

An extended Lyman $α$ outflow from a radio galaxy at z=3.7?

Authors: Miguel Coloma Puga, Barbara Balmaverde, Alessandro Capetti, Francesco Massaro, Cristina Ramos Almeida, George Miley, Roberto Gilli, Alessandro Marconi

Abstract: Spatially resolved observations of AGN host galaxies undergoing feedback processes are one of the most relevant avenues through which galactic evolution can be studied, given the long lasting effects AGN feedback has on gas reservoirs, star formation, and AGN environments at all scales. Within this context we report results from VLT/MUSE integral field optical spectroscopy of TN J1049-1258, one of… ▽ More Spatially resolved observations of AGN host galaxies undergoing feedback processes are one of the most relevant avenues through which galactic evolution can be studied, given the long lasting effects AGN feedback has on gas reservoirs, star formation, and AGN environments at all scales. Within this context we report results from VLT/MUSE integral field optical spectroscopy of TN J1049-1258, one of the most powerful radio sources known, at a redshift of 3.7. We detected extended ($\sim$ 18 kpc) Lyman $α$ emission, spatially aligned with the radio axis, redshifted by 2250 $\pm$ 60 km s$^{-1}$ with respect to the host galaxy systemic velocity, and co-spatial with UV continuum emission. This Lyman $α$ emission could arise from a companion galaxy, although there are arguments against this interpretation. Alternatively, it might correspond to an outflow of ionized gas stemming from the radio galaxy. The outflow would be the highest redshift spatially resolved ionized outflow to date. The enormous amount of energy injected, however, appears to be unable to quench the host galaxy's prodigious star formation, occurring at a rate of $\sim$4500 M$_{\odot} yr^{-1}$, estimated using its far infra-red luminosity. Within the field we also found two companion galaxies at projected distances of $\sim$25 kpc and $\sim$60 kpc from the host, which suggests the host galaxy is harbored within a protocluster. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: Pre-proofs version - Accepted for publication in ApJL

AGN in overdense environments at high-$z$ with AXIS

Authors: Fabio Vito, Paolo Tozzi, Roberto Gilli, Stefano Marchesi, Nico Cappelluti, Adi Foord

Abstract: Overdense regions at high redshift ($z \gtrsim 2$) are perfect laboratories to study the relations between environment and SMBH growth, and the AGN feedback processes on the surrounding galaxies and diffuse gas. In this white paper, we discuss how AXIS will 1) constrain the AGN incidence in protoclusters, as a function of parameters such as redshift, overdensity, mass of the structure; 2) search f… ▽ More Overdense regions at high redshift ($z \gtrsim 2$) are perfect laboratories to study the relations between environment and SMBH growth, and the AGN feedback processes on the surrounding galaxies and diffuse gas. In this white paper, we discuss how AXIS will 1) constrain the AGN incidence in protoclusters, as a function of parameters such as redshift, overdensity, mass of the structure; 2) search for low-luminosity and obscured AGN in the satellite galaxies of luminous QSOs at $z>6$, exploiting the large galaxy density around such biased objects; 3) probe the AGN feedback on the proto-ICM via the measurement of the AGN contribution to the gas ionization and excitation, and the detection of extended X-ray emission from the ionized gas and from radio jets; 4) discover new large-scale structures in the wide and deep AXIS surveys as spikes in the redshift distribution of X-ray sources. These goals can be achieved only with an X-ray mission with the capabilities of AXIS, ensuring a strong synergy with current and future state-of-the-art facilities in other wavelengths. This White Paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS White Papers can be found at http://axis.astro.umd.edu/ with a mission overview at https://arxiv.org/abs/2311.00780. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: 8 pages, 3 figures. This White Paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS White Papers can be found at the AXIS website (http://axis.astro.umd.edu/) with a mission overview at arXiv:2311.00780

An extremely metal poor star complex in the reionization era: Approaching Population III stars with JWST

Authors: E. Vanzella, F. Loiacono, P. Bergamini, U. Mestric, M. Castellano, P. Rosati, M. Meneghetti, C. Grillo, F. Calura, M. Mignoli, M. Bradac, A. Adamo, G. Rihtarsic, M. Dickinson, M. Gronke, A. Zanella, F. Annibali, C. Willott, M. Messa, E. Sani, A. Acebron, A. Bolamperti, A. Comastri, R. Gilli, K. I. Caputi , et al. (9 additional authors not shown)

Abstract: We present JWST/NIRSpec integral field spectroscopy (IFS) of a lensed Population III candidate stellar complex (dubbed Lensed And Pristine 1, LAP1), with a lensing-corrected stellar mass ~<10^4 Msun, absolute luminosity M_UV > -11.2 (m_UV > 35.6), confirmed at redshift 6.639 +/- 0.004. The system is strongly amplified (μ>~ 100) by straddling a critical line of the Hubble Frontier Field galaxy clus… ▽ More We present JWST/NIRSpec integral field spectroscopy (IFS) of a lensed Population III candidate stellar complex (dubbed Lensed And Pristine 1, LAP1), with a lensing-corrected stellar mass ~<10^4 Msun, absolute luminosity M_UV > -11.2 (m_UV > 35.6), confirmed at redshift 6.639 +/- 0.004. The system is strongly amplified (μ>~ 100) by straddling a critical line of the Hubble Frontier Field galaxy cluster MACS J0416. Despite the stellar continuum is currently not detected in the Hubble and JWST/NIRCam and NIRISS imaging, arclet-like shapes of Lyman and Balmer lines, Lya, Hg, Hb and Ha are detected with NIRSpec IFS with signal-to-noise ratios SNR=5-13 and large equivalent widths (>300-2000A), along with a remarkably weak [OIII]4959-5007 at SNR ~ 4. LAP1 shows a large ionizing photon production efficiency, log(ξ_{ion}[erg~Hz^{-1}])>26. From the metallicity indexes R23 = ([OIII]4959-5007 + [OII]3727) / Hb ~< 0.74 and R3 = ([OIII]5007 / Hb) = 0.55 +/- 0.14, we derive an oxygen abundance 12+log(O/H) ~< 6.3. Intriguingly, the Ha emission is also measured in mirrored sub-components where no [OIII] is detected, providing even more stringent upper limits on the metallicity if in-situ star formation is ongoing in this region (12+log(O/H) < 6, or Z < 0.002 Zsun). The formal stellar mass limit of the sub-components would correspond to ~10^{3} Msun or M_UV fainter than -10. Alternatively, such a metal-free pure line emitting region could be the first case of a fluorescing HI gas region, induced by transverse escaping ionizing radiation from a nearby star-complex. The presence of large equivalent-width hydrogen lines and the deficiency of metal lines in such a small region, make LAP1 the most metal poor star-forming region currently known in the reionization era and a promising site that may host isolated, pristine stars. △ Less

Submitted 23 May, 2023; originally announced May 2023.

Comments: 15 pages, 10 figures, 1 table. Submitted to A&A. Comments welcome

Journal ref: A&A 678, A173 (2023)

X-ray properties and obscured fraction of AGN in the J1030 Chandra field

Authors: Matilde Signorini, Stefano Marchesi, Roberto Gilli, Marcella Brusa, Andrea Comastri, Quirino D'Amato, Kazushi Iwasawa, Giorgio Lanzuisi, Giovanni Mazzolari, Marco Mignoli, Alessandro Peca, Isabella Prandoni, Paolo Tozzi, Cristian Vignali, Fabio Vito, Colin Norman

Abstract: The 500ks Chandra ACIS-I observation of the field around the $z=6.31$ quasar SDSS J1030+0524 is currently the 5th deepest extragalactic X-ray survey. The rich multi-band coverage of the field allowed for an effective identification and redshift determination of the X-ray source counterparts: to date a catalog of 243 extragalactic X-ray sources with either a spectroscopic or photometric redshift es… ▽ More The 500ks Chandra ACIS-I observation of the field around the $z=6.31$ quasar SDSS J1030+0524 is currently the 5th deepest extragalactic X-ray survey. The rich multi-band coverage of the field allowed for an effective identification and redshift determination of the X-ray source counterparts: to date a catalog of 243 extragalactic X-ray sources with either a spectroscopic or photometric redshift estimate in the range $z\approx0-6$ is available over a 355 arcmin$^2$ area. Given its depth and the multi-band information, this catalog is an excellent resource to investigate X-ray spectral properties of distant Active Galactic Nuclei (AGN) and derive the redshift evolution of their obscuration. We performed a thorough X-ray spectral analysis for each object in the sample, measuring its nuclear column density $N_{\rm H}$ and intrinsic (de-absorbed) 2-10 keV rest-frame luminosity, $L_{2-10}$. Whenever possible, we also used the presence of the Fe K$_α$ emission line to improve the photometric redshift estimates. We measured the fractions of AGN hidden by column densities in excess of $10^{22}$ and $10^{23}$cm$^{-2}$ ($f_{22}$ and $f_{23}$, respectively) as a function of $L_{2-10}$ and redshift, and corrected for selection effects to recover the intrinsic obscured fractions. At $z\sim 1.2$, we found $f_{22}\sim0.7-0.8$ and $f_{23}\sim0.5-0.6$, respectively, in broad agreement with the results from other X-ray surveys. No significant variations with X-ray luminosity were found within the limited luminosity range probed by our sample (log$L_{2-10}\sim 42.8-44.3$). When focusing on luminous AGN with log$L_{2-10}\sim44$ to maximize the sample completeness up to large cosmological distances, we did not observe any significant change in $f_{22}$ or $f_{23}$ over the redshift range $z\sim0.8-3$. Nonetheless, the obscured fractions we measure are significantly higher than ... △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: A&A, in press

Journal ref: A&A 676, A49 (2023)

Compton-thick AGN in the NuSTAR Era X: Analysing seven local CT-AGN candidates

Authors: Dhrubojyoti Sengupta, Stefano Marchesi, Cristian Vignali, Núria Torres-Albà, Elena Bertola, Andrealuna Pizzetti, Giorgio Lanzuisi, Francesco Salvestrini, Xiurui Zhao, Massimo Gaspari, Roberto Gilli, Andrea Comastri, Alberto Traina, Francesco Tombesi, Ross Silver, Francesca Pozzi, Marco Ajello

Abstract: We present the broad-band X-ray spectral analysis (0.6-50 keV) of seven Compton-Thick active galactic nuclei (CT-AGN; line-of-sight, l.o.s., column density $>10^{24}$ cm$^{-2}$) candidates selected from the Swift-BAT 100-month catalog, using archival NuSTAR data. This work is in continuation of the on-going research of the Clemson-INAF group to classify CT-AGN candidates at redshift $z<0.05$, usin… ▽ More We present the broad-band X-ray spectral analysis (0.6-50 keV) of seven Compton-Thick active galactic nuclei (CT-AGN; line-of-sight, l.o.s., column density $>10^{24}$ cm$^{-2}$) candidates selected from the Swift-BAT 100-month catalog, using archival NuSTAR data. This work is in continuation of the on-going research of the Clemson-INAF group to classify CT-AGN candidates at redshift $z<0.05$, using physically-motivated torus models. Our results confirm that three out of seven targets are \textit{bona-fide} CT-AGN. Adding our results to the previously analysed sources using NuSTAR data, we increase the population of bona-fide CT-AGN by $\sim9\%$, bringing the total number to 35 out of 414 AGN. We also performed a comparative study using MyTorus and borus02 on the spectra in our sample, finding that both physical models are strongly consistent in the parameter space of l.o.s. column density and photon index. Furthermore, the clumpiness of the torus clouds is also investigated by separately computing the line-of-sight and average torus column densities, in each of the seven sources. Adding our results to all the previous 48 CT-AGN candidates analysed by the Clemson-INAF research team having NuSTAR observations: we find $78\%$ of the sources are likely to have a clumpy distribution of the obscuring material surrounding the accreting supermassive black hole. △ Less

Submitted 12 May, 2023; originally announced May 2023.

Comments: Accepted for publication in the Section 4. Extragalactic astronomy of Astronomy and Astrophysics, 20 pages, 10 figures

Journal ref: A&A 676, A103 (2023)

The most luminous blue quasars at 3.0<z<3.3 -- III. LBT spectra and accretion parameters

Authors: Bartolomeo Trefoloni, Elisabeta Lusso, Emanuele Nardini, Guido Risaliti, Giada Bargiacchi, Susanna Bisogni, Francesca M. Civano, Martin Elvis, Giuseppina Fabbiano, Roberto Gilli, Alessandro Marconi, Gordon T. Richards, Andrea Sacchi, Francesco Salvestrini, Matilde Signorini, Cristian Vignali

Abstract: We present the analysis of the rest frame ultraviolet and optical spectra of 30 bright blue quasars at $z\sim3$, selected to examine the suitability of AGN as cosmological probes. In our previous works, we found an unexpectedly high fraction ($\approx 25 \%$) of X-ray weak quasars in the sample. The latter sources also display a flatter UV continuum and a broader and fainter CIV profile in the arc… ▽ More We present the analysis of the rest frame ultraviolet and optical spectra of 30 bright blue quasars at $z\sim3$, selected to examine the suitability of AGN as cosmological probes. In our previous works, we found an unexpectedly high fraction ($\approx 25 \%$) of X-ray weak quasars in the sample. The latter sources also display a flatter UV continuum and a broader and fainter CIV profile in the archival UV data with respect to their X-ray normal counterparts. Here we present new observations with the LBT in both the $zJ$ (rest-frame $\simeq$2300-3100 $\rm \mathring{A}$) and the $K_S$ ($\simeq$4750-5350 $\rm \mathring{A}$) bands. We estimated black hole masses ($M_{\rm BH}$) and Eddington ratios ($λ_{\rm Edd}$) from the from the H$β$ and MgII emission lines, finding that our $z\sim3$ quasars are on average highly accreting ($\langle λ_{\rm Edd} \rangle\simeq 1.2$ and $\langle M_{\rm BH} \rangle\simeq 10^{9.7}M_\odot$), with no difference in $λ_{\rm Edd}$ or $M_{\rm BH}$ between X-ray weak and X-ray normal quasars. From the $zJ$ spectra, we derive flux and equivalent width of MgII and FeII, finding that X-ray weak quasars display higher FeII/MgII ratios with respect to typical quasars. FeII/MgII ratios of X-ray normal quasars are instead consistent with other estimates up to $z\simeq6.5$, corroborating the idea of already chemically mature BLRs at early cosmic time. From the $K_S$ spectra, we find that all the X-ray weak quasars present generally weaker [OIII] emission (EW<10 $\rm \mathring{A}$) than the normal ones. The sample as a whole, however, abides by the known X-ray/[OIII] luminosity correlation, hence the different [OIII] properties are likely due to an intrinsically weaker [OIII] emission in X-ray weak objects, associated to the shape of the spectral energy distribution. We interpret these results in the framework of accretion-disc winds. △ Less

Submitted 12 May, 2023; originally announced May 2023.

Comments: 25 pages, 17 figures, Accepted for publication on Astronomy & Astrophysics

HYPerluminous quasars at the Epoch of ReionizatION (HYPERION). A new regime for the X-ray nuclear properties of the first quasars

Authors: L. Zappacosta, E. Piconcelli, F. Fiore, I. Saccheo, R. Valiante, C. Vignali, F. Vito, M. Volonteri, M. Bischetti, A. Comastri, C. Done, M. Elvis, E. Giallongo, F. La Franca, G. Lanzuisi, M. Laurenti, G. Miniutti, A. Bongiorno, M. Brusa, F. Civano, S. Carniani, V. D'Odorico, C. Feruglio, S. Gallerani, R. Gilli , et al. (18 additional authors not shown)

Abstract: The existence of luminous quasars (QSO) at the Epoch of Reionization (EoR; i.e. z>6) powered by supermassive black holes (SMBH) with masses $\gtrsim10^9~M_\odot$ challenges models of early SMBH formation. To shed light on the nature of these sources we started a multiwavelength programme based on a sample of 18 HYPerluminous quasars at the Epoch of ReionizatION (HYPERION). These are the luminous Q… ▽ More The existence of luminous quasars (QSO) at the Epoch of Reionization (EoR; i.e. z>6) powered by supermassive black holes (SMBH) with masses $\gtrsim10^9~M_\odot$ challenges models of early SMBH formation. To shed light on the nature of these sources we started a multiwavelength programme based on a sample of 18 HYPerluminous quasars at the Epoch of ReionizatION (HYPERION). These are the luminous QSOs whose SMBH must have had the fastest mass growth during the Universe first Gyr. In this paper we present the HYPERION sample and report on the first of the 3 years planned observations of the 2.4 Ms XMM-Newton Multi-Year Heritage program on which HYPERION is based. The goal of this program is to accurately characterize the X-ray nuclear properties of QSOs at the EoR. Through a joint X-ray spectral analysis of 10 sources, in the rest-frame $\sim2-50$ keV range, we report a steep average photon index ($Γ\sim2.4\pm0.1$). Absorption is not required. The average $Γ$ is inconsistent at $\geq4σ$ level with the canonical 1.8-2 value measured in QSO at z<6. This spectral slope is also much steeper than that reported in lower-z QSOs with similar luminosity or accretion rate, thus suggesting a genuine redshift evolution. Alternatively, we can interpret this result as the presence of an unusually low-energy cutoff $E_{cut}\sim20$ keV on a standard $Γ=1.9$ power-law. We also report on mild indications that HYPERION QSOs show higher soft X-ray emission at 2 keV compared to the UV one at 2500A than expected by lower-z luminous AGN. We speculate that a redshift-dependent coupling between the corona and accretion disc or intrinsically different coronal properties may account for the steep spectral slopes, especially in the presence of powerful winds. The reported slopes, if confirmed at lower luminosities, may have an important impact on future X-ray AGN studies in the early Universe. △ Less

Submitted 19 July, 2023; v1 submitted 3 May, 2023; originally announced May 2023.

Comments: 21 pages (including appendix), 12 figures, 4 tables. Accepted for pubblication in A&A

Journal ref: A&A 678, A201 (2023)

Large-scale clustering of buried X-ray AGN: Trends in AGN obscuration and redshift evolution

Authors: Akke Viitanen, Viola Allevato, Alexis Finoguenov, Francesco Shankar, Roberto Gilli, Giorgio Lanzuisi, Fabio Vito

Abstract: In order to test active galactic nucleus (AGN) unification and evolutionary models, we measured the AGN clustering properties as a function of AGN obscuration defined in terms of hydrogen column density, $N_{\rm H}$. In addition to measuring the clustering of unobscured ($N_{\rm H} < 10^{22}\,{\rm cm}^{-2}$) and moderately obscured ($10^{22} \leq N_{\rm H} < 10^{23.5}$) AGNs, we also targeted high… ▽ More In order to test active galactic nucleus (AGN) unification and evolutionary models, we measured the AGN clustering properties as a function of AGN obscuration defined in terms of hydrogen column density, $N_{\rm H}$. In addition to measuring the clustering of unobscured ($N_{\rm H} < 10^{22}\,{\rm cm}^{-2}$) and moderately obscured ($10^{22} \leq N_{\rm H} < 10^{23.5}$) AGNs, we also targeted highly obscured sources ($N_{\rm H}\geq 10^{23.5}$) up to redshifts of $z=3$. We have compiled one of the largest samples of X-ray-selected AGNs from a total of eight deep XMM/Chandra surveys. We measured the clustering as a function of both AGN obscuration and redshift using the projected two-point correlation function, $w_{\rm p}(r_{\rm p})$. We modeled the large-scale clustering signal, measured the AGN bias, $b(z, N_{\rm H})$, and interpreted it in terms of the typical AGN host dark matter halo, $M_{\rm halo}(z, N_{\rm H}$). We find no significant dependence of AGN clustering on obscuration, suggesting similar typical masses of the hosting halos as a function of $N_{\rm H}$. This result matches expectations of AGN unification models, in which AGN obscuration depends mainly on the viewing angle of the obscuring torus. We measured, for the first time, the clustering of highly obscured AGNs and find that these objects reside in halos with typical mass $\log M_{\rm halo} = 12.98_{-0.22}^{+0.17} [h^{-1} M_\odot]$ ($12.28_{-0.19}^{+0.13}$) at low $z \sim 0.7$ (high $z \sim 1.8$) redshifts. We find that irrespective of obscuration, an increase in AGN bias with redshift is slower than the expectation for a constant halo mass and instead follows the growth rate of halos, known as the passive evolution track. This implies that for those AGNs the clustering is mainly driven by the mass growth rate of the hosting halos and galaxies across cosmic time. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Comments: 13 pages, 6 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 674, A214 (2023)

Measuring the Cosmic X-ray Background in 3-20keV with Straylight from NuSTAR

Authors: Steven Rossland, Daniel Wik, Brian Grefenstette, Nico Cappelluti, Francesca Civano, Fabio Gastaldello, Roberto Gilli, Fiona Harrison, Ann Hornschemeier, Ryan Hickox, Roman Krivonos, Kristin Madsen, Silvano Molendi, Andrew Ptak, Daniel Stern, Andreas Zoglauer

Abstract: By characterizing the contribution of stray light to large datasets from the NuSTAR X-ray observatory collected over 2012--2017, we report a measurement of the cosmic X-ray background in the 3--20 keV energy range. These data represent $\sim20\%$ sky coverage while avoiding Galactic Ridge X-ray emission and are less weighted by deep, survey fields than previous measurements with NuSTAR. Images in… ▽ More By characterizing the contribution of stray light to large datasets from the NuSTAR X-ray observatory collected over 2012--2017, we report a measurement of the cosmic X-ray background in the 3--20 keV energy range. These data represent $\sim20\%$ sky coverage while avoiding Galactic Ridge X-ray emission and are less weighted by deep, survey fields than previous measurements with NuSTAR. Images in narrow energy bands are stacked in detector space and spatially fit with a model representing the stray light and uniform pattern expected from the cosmic X-ray background and the instrumental background, respectively. We establish baseline flux values from Earth-occulted data and validate the fitting method on stray light observations of the Crab, which further serve to calibrate the resulting spectra. We present independent spectra of the cosmic X-ray background with the FPMA and FPMB detector arrays, which are in excellent agreement with the canonical characterization by HEAO 1 and are $10\%$ lower than most subsequent measurements; $F_{\rm{3-20~keV}}^{FPMA} = 2.63 \times 10^{-11}~\rm{erg~s^{-1}~cm^{-2}~deg^{-2}}$ and $F_{\rm{3-20~keV}}^{FPMB} = 2.58 \times 10^{-11}~\rm{erg~s^{-1}~cm^{-2}~deg^{-2}}$. We discuss these results in light of previous measurements of the cosmic X-ray background and consider the impact of systematic uncertainties on our spectra. △ Less

Submitted 19 April, 2023; v1 submitted 16 April, 2023; originally announced April 2023.

Comments: 23 pages, 15 figures, 4 tables

LBT-MODS spectroscopy of high-redshift candidates in the Chandra J1030 field. A newly discovered z$\sim$2.8 large scale structure

Authors: Stefano Marchesi, Marco Mignoli, Roberto Gilli, Giovanni Mazzolari, Matilde Signorini, Marisa Brienza, Susanna Bisogni, Micol Bolzonella, Olga Cucciati, Quirino D'Amato, Alessandro Peca, Isabella Prandoni, Paolo Tozzi, Cristian Vignali, Fabio Vito, Andrea Comastri

Abstract: We present the results of a spectroscopic campaign with the Multi-Object Double Spectrograph (MODS) instrument mounted on the Large Binocular Telescope (LBT), aimed at obtaining a spectroscopic redshift for seven Chandra J1030 sources with a photometric redshift >=2.7 and optical magnitude r_AB=[24.5-26.5]. We obtained a spectroscopic redshift for five out of seven targets: all of them have z_spec… ▽ More We present the results of a spectroscopic campaign with the Multi-Object Double Spectrograph (MODS) instrument mounted on the Large Binocular Telescope (LBT), aimed at obtaining a spectroscopic redshift for seven Chandra J1030 sources with a photometric redshift >=2.7 and optical magnitude r_AB=[24.5-26.5]. We obtained a spectroscopic redshift for five out of seven targets: all of them have z_spec>=2.5, thus probing the reliability of the Chandra J1030 photometric redshifts. The spectroscopic campaign led to the serendipitous discovery of a z~2.78 large scale structure (LSS) in the J1030 field: the structure contains four X-ray sources (three of which were targeted in the LBT-MODS campaign) and two non-X-ray detected galaxies for which a VLT-MUSE spectrum was already available. The X-ray members of the LSS are hosted in galaxies that are significantly more massive (log(M_*/M_sun)=[10.0-11.1]) than those hosting the two MUSE-detected sources (log(M_*/M_sun)<10). Both observations and simulations show that massive galaxies, and particularly objects having log(M_*/M_sun)>10, are among the best tracers of large scale structures and filaments in the cosmic web. Consequently, our result can explain why X-ray-detected AGN have also been shown to be efficient tracers of large scale structures. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Comments: 16 pages, 9 Figures. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 673, A97 (2023)

Origin of the diffuse 4-8 keV emission in M82

Authors: K. Iwasawa, C. Norman, R. Gilli, P. Gandhi, M. A. Perez-Torres

Abstract: We present the first spatially resolved, X-ray spectroscopic study of the 4-8 keV diffuse emission found in the central part of the nearby starburst galaxy M82 on a few arcsecond scales. The new details that we see allow a number of important conclusions to be drawn on the nature of the hot gas and its origin as well as feedback on the ISM. We use archival data from Chandra with an exposure time o… ▽ More We present the first spatially resolved, X-ray spectroscopic study of the 4-8 keV diffuse emission found in the central part of the nearby starburst galaxy M82 on a few arcsecond scales. The new details that we see allow a number of important conclusions to be drawn on the nature of the hot gas and its origin as well as feedback on the ISM. We use archival data from Chandra with an exposure time of 570 ks. The Fexxv emission at 6.7 keV, expected from metal-enriched hot gas, is enhanced only in a limited area close to the starburst disc and is weak or almost absent over the rest of the diffuse emission, resulting in spatial variations in EW from <0.1 keV to 1.9 keV. This shows the presence of non-thermal emission due to inverse Compton scattering of the FIR photons by radio emitting cosmic ray electrons. The morphological resemblance between the diffuse X-ray, radio, and FIR emission maps supports this concept. Our decomposition of the diffuse emission spectrum indicates that ~70% of the 4-8 keV luminosity originates from the inverse Compton emission. The metal-rich hot gas with kT~5 keV makes a minor contribution to the 4-8 keV continuum, but it accounts for the majority of the observed Fexxv line. This hot gas appears to emerge from the circumnuclear starburst ring and fill the galactic chimneys identified as mid-infrared and radio emission voids. The energetics argument suggests that much of the supernova energy in the starburst site has gone into creating the chimneys and is transported to the halo. We argue that a hot, rarefied environment produced by strong supernova feedback results in displacing the brightest X-ray and radio supernova remnants which are instead found to reside in GMC. We find a faint X-ray source with a radio counterpart, close to the kinematic centre of the galaxy and we carefully examine the possibility that this source is a low-luminosity AGN in ADAF phase. △ Less

Submitted 16 March, 2023; originally announced March 2023.

Comments: 19 pages, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 674, A77 (2023)

X-ray Unveiling Events in a z~1.6 Active Galactic Nucleus in the 7 Ms Chandra Deep Field-South

Authors: Li-Ming Yu, Bin Luo, William N. Brandt, Franz E. Bauer, Demetra De Cicco, Andy C. Fabian, Roberto Gilli, Anton M. Koekemoer, Maurizio Paolillo, Donald P. Schneider, Ohad Shemmer, Paolo Tozzi, Jonathan R. Trump, Cristian Vignali, Fabio Vito, Jun-Xian Wang, Yongquan Xue

Abstract: We investigate the extreme X-ray variability of a z = 1.608 active galactic nucleus in the 7 Ms Chandra Deep Field-South (XID 403), which showed two significant X-ray brightening events. In the first event, XID 403 brightened by a factor of $>2.5$ in $\lesssim6.1$ rest-frame days in the observed-frame 0.5-5 keV band. The event lasted for $\approx5.0\textrm{-}7.3$ days, and then XID 403 dimmed by a… ▽ More We investigate the extreme X-ray variability of a z = 1.608 active galactic nucleus in the 7 Ms Chandra Deep Field-South (XID 403), which showed two significant X-ray brightening events. In the first event, XID 403 brightened by a factor of $>2.5$ in $\lesssim6.1$ rest-frame days in the observed-frame 0.5-5 keV band. The event lasted for $\approx5.0\textrm{-}7.3$ days, and then XID 403 dimmed by a factor of $>6.0$ in $\lesssim6.1$ days. After $\approx1.1\textrm{-}2.5$ years in the rest frame (including long observational gaps), it brightened again with the 0.5-5 keV flux increasing by a factor of $>12.6$. The second event lasted over 251 days and the source remained bright until the end of the 7 Ms exposure. The spectrum is a steep power law (photon index $Γ=2.8\pm0.3$) without obscuration during the second outburst, and the rest-frame 2-10 keV luminosity reaches $1.5^{+0.8}_{-0.5}\times10^{43}$ erg s$^{-1}$; there is no significant spectral evolution within this epoch. The infrared-to-UV spectral energy distribution of XID 403 is dominated by the host galaxy. There is no significant optical/UV variability and $R$-band (rest-frame $\approx2500$ $\unicode{xC5}$) brightening contemporaneous with the X-ray brightening. The extreme X-ray variability is likely due to two X-ray unveiling events, where the line of sight to the corona is no longer shielded by high-density gas clumps in a small-scale dust-free absorber. XID 403 is probably a high-redshift analog of local narrow-line Seyfert 1 galaxies, and the X-ray absorber is a powerful accretion-disk wind. On the other hand, we cannot exclude the possibility that XID 403 is an unusual candidate for tidal disruption events. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Comments: 23 pages, 9 figures, accepted for publication in ApJ

AGN feedback in an infant galaxy cluster: the LOFAR-Chandra view of the giant FRII radio galaxy J103025+052430 at z=1.7

Authors: M. Brienza, R. Gilli, I. Prandoni, Q. D'Amato, K. Rajpurohit, F. Calura, M. Chiaberge, A. Comastri, K. Iwasawa, G. Lanzuisi, E. Liuzzo, S. Marchesi, M. Mignoli, G. Miley, C. Norman, A. Peca, M. Raciti, T. Shimwell, P. Tozzi, C. Vignali, F. Vitello, F. Vito

Abstract: In the nearby universe jets from AGN are observed to have a dramatic impact on their surrounding extragalactic environment. Their effect at the `cosmic noon' (z>1.5), the epoch when star formation and AGN activity peak, is instead much less constrained. Here we present a study of the giant (750 kpc) radio galaxy 103025+052430 located at the centre of a protocluster at redshift z=1.7, with a focus… ▽ More In the nearby universe jets from AGN are observed to have a dramatic impact on their surrounding extragalactic environment. Their effect at the `cosmic noon' (z>1.5), the epoch when star formation and AGN activity peak, is instead much less constrained. Here we present a study of the giant (750 kpc) radio galaxy 103025+052430 located at the centre of a protocluster at redshift z=1.7, with a focus on its interaction with the external medium. We present new LOFAR observations at 144 MHz, which we combine with VLA 1.4 GHz and 0.5-7 keV Chandra archival data. The new map at 144 MHz confirms that the source has a complex morphology, possibly consistent with the `hybrid morphology' classification. The large size of the source gave us the possibility to perform a resolved radio spectral index analysis, a very unique opportunity for a source at such high redshift. This reveals a tentative flattening of the radio spectral index at the edge of the backflow in the Western lobe, which might be indicating plasma compression. The spatial coincidence between this region and the thermal X-ray bubble C suggests a causal connection between the two. Contrary to previous estimates for the bright X-ray component A, we find that inverse Compton scattering between the radio-emitting plasma of the Eastern lobe and the CMB photons can account for a large fraction (~45%-80%) of its total 0.5-7 keV measured flux. Finally, the X-ray bubble C, which is consistent with a thermal origin, is found to be significantly overpressurised with respect to the ambient medium. This suggests that it will tend to expand and release its energy in the surroundings, contributing to the overall intracluster medium heating. Overall, 103025+052430 gives us the chance to investigate the interaction between AGN jets and the surrounding gas in a system that is likely the predecessor of the rich galaxy clusters we all well know at z=0. △ Less

Submitted 24 February, 2023; v1 submitted 23 February, 2023; originally announced February 2023.

Comments: 12 pages, 5 figures, 2 tables, accepted for publication in A&A

The universal shape of the X-ray variability power spectrum of AGN up to $z\sim 3$

Authors: M. Paolillo, I. E. Papadakis, W. N. Brandt, F. E. Bauer, G. Lanzuisi, V. Allevato, O. Shemmer, X. C. Zheng, D. De Cicco, R. Gilli, B. Luo, M. Thomas, P. Tozzi, F. Vito, Y. Q. Xue

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

Submitted 16 February, 2023; originally announced February 2023.

Comments: In press on Astronomy and Astrophysics

Journal ref: A&A 673, A68 (2023)

Supermassive Black Hole Winds in X-rays -- SUBWAYS. I. Ultra-fast outflows in QSOs beyond the local Universe

Authors: G. A. Matzeu, M. Brusa, G. Lanzuisi, M. Dadina, S. Bianchi, G. Kriss, M. Mehdipour, E. Nardini, G. Chartas, R. Middei, E. Piconcelli, V. Gianolli, A. Comastri, A. L. Longinotti, Y. Krongold, F. Ricci, P. O. Petrucci, F. Tombesi, A. Luminari, L. Zappacosta, G. Miniutti, M. Gaspari, E. Behar, M. Bischetti, S. Mathur , et al. (26 additional authors not shown)

Abstract: We present a new X-ray spectroscopic study of $22$ luminous ($2\times10^{45}\lesssim L_{\rm bol}\rm /erg\,s^{-1} \lesssim 2\times10^{46}$) active galactic nuclei (AGNs) at intermediate-redshift ($0.1 \lesssim z \lesssim 0.4$), as part of the SUpermassive Black hole Winds in the x-rAYS (SUBWAYS) sample, mostly composed of quasars (QSOs) and type\,1 AGN. Here, 17 targets were observed with \textit{X… ▽ More We present a new X-ray spectroscopic study of $22$ luminous ($2\times10^{45}\lesssim L_{\rm bol}\rm /erg\,s^{-1} \lesssim 2\times10^{46}$) active galactic nuclei (AGNs) at intermediate-redshift ($0.1 \lesssim z \lesssim 0.4$), as part of the SUpermassive Black hole Winds in the x-rAYS (SUBWAYS) sample, mostly composed of quasars (QSOs) and type\,1 AGN. Here, 17 targets were observed with \textit{XMM-Newton} between 2019--2020 and the remaining 5 are from previous observations. The aim of this large campaign ($1.45\,\rm Ms$ duration) is to characterise the various manifestations of winds in the X-rays driven from supermassive black holes in AGN. In this paper we focus on the search and characterization of ultra-fast outflows (UFOs), which are typically detected through blueshifted absorption troughs in the Fe\,K band ($E>7\,\rm keV$). By following Monte Carlo procedures, we confirm the detection of absorption lines corresponding to highly ionised iron (e.g., Fe\,\textsc{xxv}\,H$α$, Fe\,\textsc{xxvi}\,Ly$α$) in 7/22 sources at the $\gtrsim95\%$ confidence level (for each individual line). The global combined probability of such absorption features in the sample is $>99.9\%$. The SUBWAYS campaign extends at higher luminosity and redshifts than previous local studies on Seyferts, obtained using \xmm and \suzaku observations. We find a UFO detection fraction of $\sim30\%$ on the total sample that is in agreement with the previous findings. This work independently provides further support for the existence of highly-ionised matter propagating at mildly relativistic speed ($\gtrsim0.1c$) in a considerable fraction of AGN over a broad range of luminosities, which is expected to play a key role in the self-regulated AGN feeding-feedback cycle, as also supported by hydrodynamical multiphase simulations. △ Less

Submitted 6 December, 2022; originally announced December 2022.

Comments: 30 pages, 10 figures, accepted for publication in A&A

The cavity of 3CR 196.1: H$α$ emission spatially associated with an X-ray cavity

Authors: A. Jimenez-Gallardo, E. Sani, F. Ricci, C. Mazzucchelli, B. Balmaverde, F. Massaro, A. Capetti, W. R. Forman, R. P. Kraft, G. Venturi, M. Gendron-Marsolais, M. A. Prieto, A. Marconi, H. A. Peña-Herazo, S. A. Baum, C. P. O'Dea, L. Lovisari, R. Gilli, E. Torresi, A. Paggi, V. Missaglia, G. R. Tremblay, B. J. Wilkes

Abstract: We present a multifrequency analysis of the radio galaxy 3CR 196.1 ($z = 0.198$), associated with the brightest galaxy of the cool core cluster CIZAJ0815.4-0303. This nearby radio galaxy shows a hybrid radio morphology and an X-ray cavity, all signatures of a turbulent past activity, potentially due to merger events and AGN outbursts. We present results of the comparison between $Chandra$ and VLT/… ▽ More We present a multifrequency analysis of the radio galaxy 3CR 196.1 ($z = 0.198$), associated with the brightest galaxy of the cool core cluster CIZAJ0815.4-0303. This nearby radio galaxy shows a hybrid radio morphology and an X-ray cavity, all signatures of a turbulent past activity, potentially due to merger events and AGN outbursts. We present results of the comparison between $Chandra$ and VLT/MUSE data for the inner region of the galaxy cluster, on a scale of tens of kpc. We discovered H$α$ + [N II]$\lambda6584$ emission spatially associated with the X-ray cavity (at $\sim$10 kpc from the galaxy nucleus) instead of with its rim. This result differs from previous discoveries of ionized gas surrounding X-ray cavities in other radio galaxies harbored in galaxy clusters and could represent the first reported case of ionized gas filling an X-ray cavity, either due to different AGN outbursts or to the cooling of warm ($10^4<T\leq10^7$ K) AGN outflows. We also found that the H$α$, [N II]$λ\lambda6548,6584$ and [S II]$λ\lambda6718,6733$ emission lines show an additional redward component, at $\sim$1000 km$\,$s$^{-1}$ from rest frame, with no detection in H$β$ or [O III]$λ\lambda4960,5008$. We believe the most likely explanation for this redward component is the presence of a background gas cloud since there appears to be a discrete difference in velocities between this component and the rest frame. △ Less

Submitted 8 November, 2022; originally announced November 2022.

Comments: 15 pages, 8 figures, ApJ accepted, pre-proof version

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

Authors: Q. D'Amato, I. Prandoni, R. Gilli, C. Vignali, M. Massardi, E. Liuzzo, P. Jagannathan, M. Brienza, R. Paladino, M. Mignoli, S. Marchesi, A. Peca, M. Chiaberge, G. Mazzolari, C. Norman

Abstract: We present deep L-Band observations of the equatorial field centered on the z=6.3 SDSS QSO, reaching a 1 sigma sensitivity of ~2.5 uJy at the center of the field. We extracted a catalog of 1489 radio sources down to a flux density of ~12.5 uJy (5 sigma) over a field of view of ~ 30' diameter. We derived the source counts accounting for catalog reliability and completeness, and compared them with o… ▽ More We present deep L-Band observations of the equatorial field centered on the z=6.3 SDSS QSO, reaching a 1 sigma sensitivity of ~2.5 uJy at the center of the field. We extracted a catalog of 1489 radio sources down to a flux density of ~12.5 uJy (5 sigma) over a field of view of ~ 30' diameter. We derived the source counts accounting for catalog reliability and completeness, and compared them with others available in the literature. Our source counts are among the deepest available so far, and, overall, are consistent with recent counts' determinations and models. We detected for the first time in the radio band the SDSS J1030+0524 QSO (26 +/- 5 uJy). We derived its optical radio loudness R_O = 0.62 +/- 0.12, which makes it the most radio quiet AGN at z >~ 6 discovered so far and detected at radio wavelengths. We unveiled extended diffuse radio emission associated with the lobes of a bright FRII radio galaxy located close to the center of the J1030 field, which is likely to become the future BCG of a protocluster at z=1.7. The lobes' complex morphology, coupled with the presence of X-ray diffuse emission detected around the FRII galaxy lobes, may point toward an interaction between the radio jets and the external medium. We also investigated the relation between radio and X-ray luminosity for a sample of 243 X-ray-selected objects obtained from 500 ks Chandra observations of the same field, and spanning a wide redshift range (0 ~< z ~< 3). Focused on sources with a spectroscopic redshift and classification, we found that sources hosted by ETG and AGN follow Log(L_R)/Log(L_X) linear correlations with slopes of ~0.6 and ~0.8, respectively. This is interpreted as a likely signature of different efficiency in the accretion process. Finally, we found that most of these sources (>~87%) show a radio-to-X-ray radio loudness R_X < -3.5, classifying these objects as radio quiet. △ Less

Submitted 27 October, 2022; originally announced October 2022.

Comments: 18 pages, 15 figures, 5 table. Accepted for publication in A&A

Journal ref: A&A 668, A133 (2022)

The 700 ks Chandra Spiderweb Field II: Evidence for inverse-Compton and thermal diffuse emission in the Spiderweb galaxy

Authors: P. Tozzi, R. Gilli, A. Liu, S. Borgani, M. Lepore, L. Di Mascolo, A. Saro, L. Pentericci, C. Carilli, G. Miley, T. Mroczkowski, M. Pannella, E. Rasia, P. Rosati, C. S. Anderson, A. Calabro', E. Churazov, H. Dannerbauer, C. Feruglio, F. Fiore, R. Gobat, S. Jin, M. Nonino, C. Norman, H. J. A. Rottgering

Abstract: We present the X-ray imaging and spectral analysis of the diffuse emission around the Spiderweb galaxy at z=2.16 and of its nuclear emission, based on a deep (700 ks) Chandra observation. We characterize the nuclear emission and computed the contamination in the surrounding regions due to the wings of the instrument PSF. Then, we quantified the extended emission within 12". We find that the Spider… ▽ More We present the X-ray imaging and spectral analysis of the diffuse emission around the Spiderweb galaxy at z=2.16 and of its nuclear emission, based on a deep (700 ks) Chandra observation. We characterize the nuclear emission and computed the contamination in the surrounding regions due to the wings of the instrument PSF. Then, we quantified the extended emission within 12". We find that the Spiderweb galaxy hosts a mildly absorbed quasar, with modest yet significant variability on a timescale of ~1 year. We find that the emission in the jet regions is well described by a power law with Gamma~2-2.5, and it is consistent with IC upscattering of the CMB photons by the relativistic electrons. We also find a roughly symmetric, diffuse emission within a radius of ~100 kpc. This emission is consistent with thermal bremsstrahlung from a hot ICM with a temperature of kT=2.0_{-0.4}^{+0.7} keV, and a metallicity of Z<1.6Z_sun. The average electron density within 100 kpc is n_e=(1.51+-0.24+-0.14)E-2 cm^{-3}, corresponding to an upper limit for the total ICM mass of <=(1.76+-0.30+-0.17)E+12 M_sun (where error bars are 1 sigma statistical and systematic, respectively). If we apply hydrostatic equilibrium to the ICM, we measure a total gravitational mass M(<100 kpc)=(1.5^{+0.5}_{-0.3})E+13 M_sun and, extrapolating at larger radii, we estimate a total mass M_{500}=(3.2^{+1.1}_{-0.6})E+13 M_sun within a radius of r_{500}=(220+-30) kpc. We conclude that the Spiderweb protocluster shows significant diffuse emission within a radius of 12 arcsec, whose major contribution is provided by IC scattering associated with the radio jets. Outside the jet regions, we also identified thermal emission within a radius of ~100 kpc, revealing the presence of hot, diffuse baryons that may represent the embryonic virialized halo of the forming cluster. △ Less

Submitted 30 September, 2022; originally announced September 2022.

Comments: 19 pages, 16 figures, accepted in Astronomy & Astrophysics

Compton-Thick AGN in the NuSTAR era VIII: A joint NuSTAR-XMM-Newton monitoring of the changing-look Compton-thick AGN NGC 1358

Authors: Stefano Marchesi, Xiurui Zhao, Núria Torres-Albà, Marco Ajello, Massimo Gaspari, Andrealuna Pizzetti, Johannes Buchner, Elena Bertola, Andrea Comastri, Anna Feltre, Roberto Gilli, Giorgio Lanzuisi, Gabriele Matzeu, Francesca Pozzi, Francesco Salvestrini, Dhrubojyoti Sengupta, Ross Silver, Francesco Tombesi, Alberto Traina, Cristian Vignali, Luca Zappacosta

Abstract: We present the multi-epoch monitoring with NuSTAR and XMM-Newton of NGC 1358, a nearby Seyfert 2 galaxy whose properties made it a promising candidate X-ray changing look AGN, i.e., a source whose column density could transition from its 2017 Compton-thick (CT-, having line-of-sight Hydrogen column density NH,los>10^24 cm^-2) state to a Compton-thin (NH,los<10^24 cm^-2) one. The multi-epoch X-ray… ▽ More We present the multi-epoch monitoring with NuSTAR and XMM-Newton of NGC 1358, a nearby Seyfert 2 galaxy whose properties made it a promising candidate X-ray changing look AGN, i.e., a source whose column density could transition from its 2017 Compton-thick (CT-, having line-of-sight Hydrogen column density NH,los>10^24 cm^-2) state to a Compton-thin (NH,los<10^24 cm^-2) one. The multi-epoch X-ray monitoring confirmed the presence of significant NH,los variability over time-scales as short as weeks, and allowed us to confirm the "changing look" nature of NGC 1358, which has most recently been observed in a Compton-thin status. Multi-epoch monitoring with NuSTAR and XMM-Newton is demonstrated to be highly effective in simultaneously constraining three otherwise highly degenerate parameters: the torus average column density and covering factor, and the inclination angle between the torus axis and the observer. We find a tentative anti-correlation between column density and luminosity, which can be understood in the framework of Chaotic Cold Accretion clouds driving recursive AGN feedback. The monitoring campaign of NGC 1358 has proven the efficiency of our newly developed method to select candidate NH,los-variable, heavily obscured AGN, which we plan to soon extend to a larger sample to better characterize the properties of the obscuring material surrounding accreting supermassive black holes, as well as constrain AGN feeding models. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 24 pages, 13 figures: main results in Figures 4, 7, 8 and 9. Accepted for publication in the Astrophysical Journal

Unresolved z~8 point sources and their impact on the bright end of the galaxy luminosity function

Authors: Yuzo Ishikawa, Takahiro Morishita, Massimo Stiavelli, Nicha Leethochawalit, Harry Ferguson, Roberto Gilli, Charlotte Mason, Michele Trenti, Tommaso Treu, Colin Norman

Abstract: The distribution and properties of the first galaxies and quasars are critical pieces of the puzzle in understanding galaxy evolution and cosmic reionization. Previous studies have often excluded unresolved sources as potential low redshift interlopers. We combine broadband color and photometric redshift analysis with morphological selections to identify a robust sample of candidates consistent wi… ▽ More The distribution and properties of the first galaxies and quasars are critical pieces of the puzzle in understanding galaxy evolution and cosmic reionization. Previous studies have often excluded unresolved sources as potential low redshift interlopers. We combine broadband color and photometric redshift analysis with morphological selections to identify a robust sample of candidates consistent with unresolved point sources at redshift $z\sim8$ using deep Hubble Space Telescope images. We also examine G141 grism spectroscopic data to identify and eliminate dwarf star contaminants. From these analyses, we identify three, bright ($M_{UV}\lesssim-22$ ABmag) dropout point sources at $7.5<z<8.1$. Spectral energy distribution analyses suggest that these sources are either quasars or compact star-forming galaxies. The flux captured by the IRAC 4.5 $μ$m channel suggests that they have moderate $Hβ$+$[OIII]$ equivalent widths. We calculate the number density of point sources at $z\sim7$-8, and find that a double powerlaw model well describes the point source distribution. We then extend our analysis to estimate the combined point source + galaxy luminosity function and find that the point sources have a non-negligible contribution to the bright-end excess. The fact that the point sources dominate only at $M_{UV}\lesssim-22$ suggests that their contribution to cosmic reionization is likely limited. While spectroscopic follow-up is needed to confirm the nature of these point sources, this work demonstrates that the inclusion of Lyman dropout point sources is necessary for a complete census of the early galaxies at the epoch of cosmic reionization. △ Less

Submitted 13 July, 2022; originally announced July 2022.

Comments: 21 pages. 10 figures

XXL-HSC: The link between AGN activity and star formation in the Early Universe ($z\geqslant3.5$)

Authors: E. Pouliasis, G. Mountrichas, I. Georgantopoulos, A. Ruiz, R. Gilli, E. Koulouridis, M. Akiyama, Y. Ueda, C. Garrel, T. Nagao, S. Paltani, M. Pierre, Y. Toba, C. Vignali

Abstract: In this work, we aimed at investigating the star formation rate of active galactic nuclei host galaxies in the early Universe. To this end, we constructed a sample of 149 luminous ($\rm L_{2-10keV} > 10^{44}\,erg\,s^{-1}$) X-ray AGNs at $\rm z \geq3.5$ selected in three fields with different depths and observed areas (Chandra COSMOS Legacy survey, XMM-XXL North and eFEDS). We built their spectral… ▽ More In this work, we aimed at investigating the star formation rate of active galactic nuclei host galaxies in the early Universe. To this end, we constructed a sample of 149 luminous ($\rm L_{2-10keV} > 10^{44}\,erg\,s^{-1}$) X-ray AGNs at $\rm z \geq3.5$ selected in three fields with different depths and observed areas (Chandra COSMOS Legacy survey, XMM-XXL North and eFEDS). We built their spectral energy distributions (SED) using available multi-wavelength photometry from X-rays up to far-IR. Then, we estimated the stellar mass, M$_{*}$, and the SFR of the AGNs using the X-CIGALE SED fitting algorithm. After applying several quality criteria, we ended up with 89 high-z sources. More than half (55\%) of the X-ray sample have spectroscopic redshifts. Based on our analysis, our high-z X-ray AGNs live in galaxies with median $\rm M_{*}=5.6 \times10^{10}~M_\odot$ and $\rm SFR_{*}\approx240\,M_\odot yr^{-1}$. The majority of the high-z sources ($\sim89$\%) were found inside or above the main sequence (MS) of star-forming galaxies. Estimation of the normalised SFR, $\rm SFR_{NORM}$, defined as the ratio of the SFR of AGNs to the SFR of MS galaxies, showed that the SFR of AGNs is enhanched by a factor of $\sim 1.8$ compared to non-AGN star-forming systems. Combining our results with previous studies at lower redshifts, we confirmed that $\rm SFR_{NORM}$ does not evolve with redshift. Using the specific BHAR (i.e., $\rm L_X$ divided by $\rm M_{*}$), $\rm λ_{BHAR}$, that can be used as a tracer of the Eddington ratio, we found that the bulk of AGNs that lie inside or above the MS have higher specific accretion rates compared to sources below the MS. Finally, we found indications that the SFR of the most massive AGN host galaxies ($\rm log\,(M_{*}/ M_\odot) >10^{11.5-12}$) remains roughly constant as a function of M$_*$, in agreement with the SFR of MS star-forming galaxies. △ Less

Submitted 28 June, 2022; originally announced June 2022.

Comments: Accepted for publication in A&A. 19 pages, 12 figures (+6 in Appendix), 4 tables

Journal ref: A&A 667, A56 (2022)

Sub-parsec resolution cosmological simulations of star-forming clumps at high redshift with feedback of individual stars

Authors: F. Calura, A. Lupi, J. Rosdahl, E. Vanzella, M. Meneghetti, P. Rosati, E. Vesperini, E. Lacchin, R. Pascale, R. Gilli

Abstract: We introduce a new set of zoom-in cosmological simulations with sub-pc resolution, intended to model extremely faint, highly magnified star-forming stellar clumps, detected at z=6.14 thanks to gravitational lensing. The simulations include feedback from individual massive stars (in both the pre-supernova and supernova phases), generated via stochastic, direct sampling of the stellar initial mass f… ▽ More We introduce a new set of zoom-in cosmological simulations with sub-pc resolution, intended to model extremely faint, highly magnified star-forming stellar clumps, detected at z=6.14 thanks to gravitational lensing. The simulations include feedback from individual massive stars (in both the pre-supernova and supernova phases), generated via stochastic, direct sampling of the stellar initial mass function. We adopt a modified 'delayed cooling' feedback scheme, specifically created to prevent artificial radiative loss of the energy injected by individual stars in very dense gas (n~10^3-10^5 cm^{-3}). The sites where star formation ignites are characterised by maximum densities of the order of 10^5 cm^{-3} and gravitational pressures P/k>10^7 K/cm^3, corresponding to the values of the local, turbulent regions where the densest stellar aggregates form. The total stellar mass at z=6.14 is 3.4x10^7 M_sun, in satisfactory agreement with the observed stellar mass of the observed systems. The most massive clumps have masses of ~10^6 M_sun and half-mass sizes of ~100 pc. These sizes are larger than the observed ones, including also other samples of lensed high-redshift clumps, and imply an average density one order of magnitude lower than the observed one. In the size-mass plane, our clumps populate a sequence that is intermediate between the ones of observed high-redshift clumps and local dSph galaxies. △ Less

Submitted 16 August, 2022; v1 submitted 27 June, 2022; originally announced June 2022.

Comments: MNRAS, accepted for publication. 24 pages, 17 figures

Quasars as high-redshift standard candles

Authors: A. Sacchi, G. Risaliti, M. Signorini, E. Lusso, E. Nardini, G. Bargiacchi, S. Bisogni, F. Civano, M. Elvis, G. Fabbiano, R. Gilli, B. Trefoloni, C. Vignali

Abstract: The non-linear relation between the X-ray and ultraviolet (UV) luminosity in quasars has been used to derive quasar distances and to build a Hubble diagram at redshifts up to $z\sim$ 7. This cosmological application is based on the assumption of independence of the relation on redshift and luminosity. We want to test the reliability of this hypothesis by studying the spectroscopic properties of hi… ▽ More The non-linear relation between the X-ray and ultraviolet (UV) luminosity in quasars has been used to derive quasar distances and to build a Hubble diagram at redshifts up to $z\sim$ 7. This cosmological application is based on the assumption of independence of the relation on redshift and luminosity. We want to test the reliability of this hypothesis by studying the spectroscopic properties of high-redshift quasars in the X-ray and UV bands. We performed a one-by-one analysis of a sample of 130 quasars at $z>$ 2.5 with high-quality X-ray and UV spectroscopic observations. We found that not only the X-ray to UV correlation still holds at these redshifts, but its intrinsic dispersion is as low as 0.12 dex (previous works reached 0.20$-$0.22 dex). For a sample of quasars at $z\sim$ 3 with particularly high-quality observations the dispersion further drops to 0.09 dex, a value entirely accountable for by intrinsic variability and source geometry effects. The composite spectra of these quasars, in both the X-rays and the UV, do not show any difference with respect to the average spectra of quasars at lower redshifts. The absence of any spectral difference between high- and low-$z$ quasars and the tightness of the X-ray to UV relation suggests that no evolutionary effects are present in the relation. Therefore, it can be safely employed to derive quasar distances. Under this assumption, we obtain a measurement of the luminosity distance at $z\sim$ 3 with 15 % uncertainty, and in a 4$σ$ tension with the concordance model. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: Accepted for publication in Astronomy & Astrophysics Letters, 8 pages, 6 figures

Journal ref: A&A 663, L7 (2022)

An X-ray fading, UV brightening QSO at $z\approx6$

Authors: Fabio Vito, Marco Mignoli, Roberto Gilli, William Nielsen Brandt, Ohad Shemmer, Franz Erik Bauer, Susanna Bisogni, Bin Luo, Stefano Marchesi, Riccardo Nanni, Gianni Zamorani, Andrea Comastri, Felice Cusano, Simona Gallerani, Cristian Vignali, Giorgio Lanzuisi

Abstract: Explaining the existence of $\gtrsim10^8\,\mathrm{M_\odot}$ SMBHs at $z>6$ is a persistent challenge to modern astrophysics. Multi-wavelength observations of $z\gtrsim6$ QSOs reveal that, on average, their accretion physics is similar to that of their counterparts at lower redshift. However, QSOs showing properties that deviate from the general behavior can provide useful insights into the physica… ▽ More Explaining the existence of $\gtrsim10^8\,\mathrm{M_\odot}$ SMBHs at $z>6$ is a persistent challenge to modern astrophysics. Multi-wavelength observations of $z\gtrsim6$ QSOs reveal that, on average, their accretion physics is similar to that of their counterparts at lower redshift. However, QSOs showing properties that deviate from the general behavior can provide useful insights into the physical processes responsible for the rapid growth of SMBHs in the early universe. We present X-ray (XMM-Newton, 100 ks) follow-up observations of a $z\approx6$ QSO, J1641+3755, which was found to be remarkably X-ray bright in a 2018 Chandra dataset. J1641+3755 is not detected in the 2021 XMM-Newton observation, implying that its X-ray flux decreased by a factor $\gtrsim7$ on a notably short timescale (i.e., $\approx115$ rest-frame days), making it the $z>4$ QSO with the largest variability amplitude. We also obtained rest-frame UV spectroscopic and photometric data with textit{LBT}, and compared them with archival datasets. Surprisingly, we found that J1641+3755 became brighter in the rest-frame UV band from 2003 to 2016, while no strong variation occurred from 2016 to 2021. Multiple narrow absorption features are detected in its rest-frame UV spectrum, and several of them can be associated with an intervening system at $z=5.67$. The variability properties of J1641+3755 can be due to intrinsic variations of the accretion rate, a small-scale obscuration event, gravitational lensing due to an intervening object, or an unrelated X-ray transient in a foreground galaxy in 2018. Accounting for all of the $z>6$ QSOs with multiple X-ray observations separated by $>10$ rest-frame days, we found an enhancement of strongly (i.e., by a factor $>3$) X-ray variable objects compared to QSOs at later cosmic times. This finding may be related to the physics of fast accretion in high-redshift QSOs. △ Less

Submitted 10 June, 2022; originally announced June 2022.

Comments: 14 pages. Accepted for publication on A&A

Journal ref: A&A 663, A159 (2022)

Supermassive Black Holes at High Redshift are Expected to be Obscured by their Massive Host Galaxies' Inter Stellar Medium

Authors: R. Gilli, C. Norman, F. Calura, F. Vito, R. Decarli, S. Marchesi, K. Iwasawa, A. Comastri, G. Lanzuisi, F. Pozzi, Q. D'Amato, C. Vignali, M. Brusa, M. Mignoli, P. Cox

Abstract: We combine results from deep ALMA observations of massive ($M_*>10^{10}\;M_{\odot}$) galaxies at different redshifts to show that the column density of their inter stellar medium (ISM) rapidly increases towards early cosmic epochs. Our analysis includes objects from the ASPECS and ALPINE large programs, as well as individual observations of $z\sim 6$ QSO hosts. When accounting for non-detections a… ▽ More We combine results from deep ALMA observations of massive ($M_*>10^{10}\;M_{\odot}$) galaxies at different redshifts to show that the column density of their inter stellar medium (ISM) rapidly increases towards early cosmic epochs. Our analysis includes objects from the ASPECS and ALPINE large programs, as well as individual observations of $z\sim 6$ QSO hosts. When accounting for non-detections and correcting for selection effects, we find that the median surface density of the ISM of the massive galaxy population evolves as $\sim(1+z)^{3.3}$. This means that the ISM column density towards the nucleus of a $z>3$ galaxy is typically $>100$ times larger than locally, and it may reach values as high as Compton-thick at $z\gtrsim6$. Remarkably, the median ISM column density is of the same order of what is measured from X-ray observations of large AGN samples already at $z\gtrsim2$. We develop a simple analytic model for the spatial distribution of ISM clouds within galaxies, and estimate the total covering factor towards active nuclei when obscuration by ISM clouds on the host scale is added to that of pc-scale circumnuclear material (the so-called 'torus'). The model includes clouds with a distribution of sizes, masses, and surface densities, and also allows for an evolution of the characteristic cloud surface density with redshift, $Σ_{c,*}\propto(1+z)^γ$. We show that, for $γ=2$, such a model successfully reproduces the increase of the obscured AGN fraction with redshift that is commonly observed in deep X-ray surveys, both when different absorption thresholds and AGN luminosities are considered. Our results suggest that 80-90\% of supermassive black holes in the early Universe ($z>6-8$) are hidden to our view, primarily by the ISM in their hosts. [abridged] △ Less

Submitted 22 July, 2022; v1 submitted 7 June, 2022; originally announced June 2022.

Comments: 22 pages, 15 figures, accepted for publication in A&A

Journal ref: A&A 666, A17 (2022)

Feedback effect on the observable properties of $z>6$ AGN

Authors: Fabio Vito, Fabio Di Mascia, Simona Gallerani, Tommaso Zana, Andrea Ferrara, Stefano Carniani, Roberto Gilli

Abstract: Active galactic nuclei (AGN) feedback has a major impact onto the supermassive black-hole (SMBH) growth, the properties of the host galaxies, and their cosmic evolution. We investigate the effects of different kinetic feedback prescriptions on the observable properties of AGN and their host galaxies at $z>6$ in a suite of zoom-in cosmological simulations. We find that kinetic feedback decreases th… ▽ More Active galactic nuclei (AGN) feedback has a major impact onto the supermassive black-hole (SMBH) growth, the properties of the host galaxies, and their cosmic evolution. We investigate the effects of different kinetic feedback prescriptions on the observable properties of AGN and their host galaxies at $z>6$ in a suite of zoom-in cosmological simulations. We find that kinetic feedback decreases the column density of the interstellar medium (ISM) in the host galaxy by up to a factor of $\approx10$, especially when the SMBHs reach high accretion rates ($\approx10-30\,\mathrm{M_\odot\,yr^{-1}}$). In particular, kinetic feedback is required to extend the ISM size to $>1$ kpc and match the observed sizes of the gas reservoirs in $z>6$ AGN host galaxies. Moreover, it produces unobscured lines of sight along which the AGN can be detected in the rest-frame UV band with magnitudes consistent with observed values of $z>6$ AGN. The assumed geometry of the outflow plays an important role in shaping the observed properties of high-redshift AGN. We find that a biconical geometry is favored over a spherical one to reproduce the observed properties, but it overestimates the number of multiple AGN systems detectable in X-ray observations. This result suggests that simplistic BH seeding recipes widely employed in cosmological simulations produce too many X-ray detectable multiple AGN at $z=6-7$, thus soliciting the adoption of more physically motivated seeding prescriptions. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: Accepted for publication on MNRAS, 18 pages

The properties of the X-ray corona in the distant ($z=3.91$) quasar APM 08279+5255

Authors: E. Bertola, C. Vignali, G. Lanzuisi, M. Dadina, M. Cappi, R. Gilli, G. A. Matzeu, G. Chartas, E. Piconcelli, A. Comastri

Abstract: We present new joint XMM-Newton and NuSTAR observations of APM 08279+5255, a gravitationally-lensed, broad-absorption line quasar ($z=3.91$). After showing a fairly stable flux ($f_{\rm2-10}\simeq4-5.5\times10^{-13}\rm~erg~s^{-1}$) from 2000 to 2008, APM 08279+5255 was found in a fainter state in the latest X-ray exposures ($f_{\rm2-10}\simeq2.7\times10^{-13}\rm~erg~s^{-1}$), which can likely be a… ▽ More We present new joint XMM-Newton and NuSTAR observations of APM 08279+5255, a gravitationally-lensed, broad-absorption line quasar ($z=3.91$). After showing a fairly stable flux ($f_{\rm2-10}\simeq4-5.5\times10^{-13}\rm~erg~s^{-1}$) from 2000 to 2008, APM 08279+5255 was found in a fainter state in the latest X-ray exposures ($f_{\rm2-10}\simeq2.7\times10^{-13}\rm~erg~s^{-1}$), which can likely be ascribed to a lower X-ray activity. Moreover, the 2019 data present a prominent Fe K$α$ emission line and do not show any significant absorption line. This fainter state, coupled to the first hard X-ray sampling of APM 08279+5255, allowed us to measure X-ray reflection and the high-energy cutoff in this source for the first time. From the analysis of previous XMM-Newton and Chandra observations, X-ray reflection is demonstrated to be a long-lasting feature of this source, but less prominent prior to 2008, possibly due to a stronger primary emission. The estimated high-energy cutoff ($E_{\rm cut}=99_{-35}^{+91}$ keV) sets a new redshift record for the farthest ever measured and places APM 08279+5255 in the allowed region of the compactness-temperature diagram of X-ray coronae, in agreement with previous results on high-$z$ quasars. △ Less

Submitted 2 May, 2022; originally announced May 2022.

Comments: Accepted for publication in A&A. 14 pages, 4 figures, 6 tables

Journal ref: A&A 662, A98 (2022)

The MURALES survey. VI. Properties and origin of the extended line emission structures in radio galaxies

Authors: Barbara Balmaverde, Alessandro Capetti, R. D. Baldi, S. Baum, M. Chiaberge, R. Gilli, Ana Jimenez-Gallardo, Alessandro Marconi, Francesco Massaro, E. Meyer, C. O'Dea, G. Speranza, E. Torresi, Giacomo Venturi

Abstract: This is the sixth paper presenting the results of the MUse RAdio Loud Emission line Snapshot survey (MURALES). We observed 37 radio sources from the 3C sample with z<0.3 and declination <20 degrees with the MUSE optical integral field spectrograph at the VLT. We here focus on the properties of the extended emission line regions (EELRs) that can be studied with unprecedented detail thanks to the de… ▽ More This is the sixth paper presenting the results of the MUse RAdio Loud Emission line Snapshot survey (MURALES). We observed 37 radio sources from the 3C sample with z<0.3 and declination <20 degrees with the MUSE optical integral field spectrograph at the VLT. We here focus on the properties of the extended emission line regions (EELRs) that can be studied with unprecedented detail thanks to the depth of these observations. Line emission in the 10 FRIs is, in most cases, confined to within <4 kpc, while large-scale (>4 kpc) ionized gas is seen in all but two of the 26 FRIIs. It usually takes the form of elongated or filamentary structures, typically extending between 10 and 30 kpc, but also reaching distances of ~80 kpc. We find that 1) the large-scale ionized gas structures show a tendency to be oriented at large angles from the radio axis, and 2) the gas on a scale of a few kpc from the nucleus often shows ordered rotation with a kinematical axis forming a median angle of 65 degrees with the radio axis. We also discuss the velocity field and ionization properties of the EELRs. The observed emission line structures appear to be associated with gaseous "superdisks" formed after a gas rich merger. The different properties of the EELR can be explained with a combination of the source evolutionary state and the orientation of the "superdisk" with respect to the radio axis. The general alignment between the superdisks and radio axis might be produced by stable and coherent accretion maintained over long timescales. △ Less

Submitted 1 April, 2022; originally announced April 2022.

Comments: Pre-proofs version - Accepted for publication in A&A

Journal ref: A&A 662, A23 (2022)

The 700 ks Chandra Spiderweb Field I: evidence for widespread nuclear activity in the Protocluster

Authors: P. Tozzi, L. Pentericci, R. Gilli, M. Pannella, F. Fiore, G. Miley, M. Nonino, H. J. A. Rottgering, V. Strazzullo, C. S. Anderson, S. Borgani, A. Calabro', C. Carilli, H. Dannerbauer, L. Di Mascolo, C. Feruglio, R. Gobat, S. Jin, A. Liu, T. Mroczkowski, C. Norman, E. Rasia, P. Rosati, A. Saro

Abstract: (Abridged) We present an analysis of the 700 ks Chandra ACIS-S observation of the field around the Spiderweb Galaxy at z=2.156, focusing on the nuclear activity in the associated large-scale environment. We identify unresolved X-ray sources down to flux limits of 1.3X10^{-16} and 3.9X10^{-16} erg/s/cm^2 in the soft and hard band, respectively. We search for counterparts in the optical, NIR and sub… ▽ More (Abridged) We present an analysis of the 700 ks Chandra ACIS-S observation of the field around the Spiderweb Galaxy at z=2.156, focusing on the nuclear activity in the associated large-scale environment. We identify unresolved X-ray sources down to flux limits of 1.3X10^{-16} and 3.9X10^{-16} erg/s/cm^2 in the soft and hard band, respectively. We search for counterparts in the optical, NIR and submm bands to identify X-ray sources belonging to the protocluster. We detect 107 X-ray unresolved sources within 5 arcmin (corresponding to 2.5 Mpc) of J1140-2629, among which 13 have optical counterparts with spectroscopic redshift 2.11<z<2.20, and 1 source with photometric redshift consistent with this range. Our X-ray spectral analysis shows that their intrinsic spectral slope is consistent with an average <Γ>~1.84+-0.04. The best-fit intrinsic absorption for 5 protocluster X-ray members is N_H>10^{23} cm^{-2}, while other 6 have upper limits of the order of fewX10^{22} cm^{-2}. Two sources can only be fitted with very flat Γ<=1, and are therefore considered Compton-thick candidates. Their 0.5-10 keV rest frame luminosities are larger than 2X10^{43} erg/s, significantly greater than X-ray luminosities expected from star formation activity. The X-ray luminosity function of AGN in the volume associated to the Spiderweb protocluster in the range 10^{43}<L_X<10^{44.5} erg/s, is at least 10 times higher than that in the field at the same redshift and significantly flatter. The X-ray AGN fraction is measured to be (25.5+-4.5)% in the stellar mass range log(M*/M_sun)>10.5, corresponding to an enhancement of 6.0^{+9.0}_{-3.0} with respect to the COSMOS field at comparable redshifts and stellar mass range. We conclude that the galaxy population in the Spiderweb Protocluster is characterized by enhanced X-ray nuclear activity triggered by environmental effects on Mpc scales. △ Less

Submitted 4 April, 2022; v1 submitted 4 March, 2022; originally announced March 2022.

Comments: Astronomy & Astrophysics. Minor changes in this version. Version accepted in A&A

Journal ref: A&A 662, A54 (2022)

Multi-Wavelength Study of a Proto-BCG at z = 1.7

Authors: Q. D'Amato, I. Prandoni, M. Brienza, R. Gilli, C. Vignali, R. Paladino, F. Loi, M. Massardi, M. Mignoli, S. Marchesi, A. Peca, P. Jagannathan

Abstract: In this work we performed a spectral energy distribution (SED) analysis in the optical/infrared band of the host galaxy of a proto-brightest cluster galaxy (BCG, NVSS J103023+052426) in a proto-cluster at z = 1.7. We found that it features a vigorous star formation rate (SFR) of ${\sim}$570 $\mathrm{M_{\odot}}$/yr and a stellar mass of $M_{\ast} \sim 3.7 \times 10^{11}$ $\mathrm{M_{\odot}}$; the h… ▽ More In this work we performed a spectral energy distribution (SED) analysis in the optical/infrared band of the host galaxy of a proto-brightest cluster galaxy (BCG, NVSS J103023+052426) in a proto-cluster at z = 1.7. We found that it features a vigorous star formation rate (SFR) of ${\sim}$570 $\mathrm{M_{\odot}}$/yr and a stellar mass of $M_{\ast} \sim 3.7 \times 10^{11}$ $\mathrm{M_{\odot}}$; the high corresponding specific SFR = $1.5 \pm 0.5$ $\mathrm{Gyr^{-1}}$ classifies this object as a starburst galaxy that will deplete its molecular gas reservoir in $\sim$ $3.5 \times 10^8$ yr. Thus, this system represents a rare example of a proto-BCG caught during the short phase of its major stellar mass assembly. Moreover, we investigated the nature of the host galaxy emission at 3.3 mm. We found that it originates from the cold dust in the interstellar medium, even though a minor non-thermal AGN contribution cannot be completely ruled out. Finally, we studied the polarized emission of the lobes at 1.4 GHz. We unveiled a patchy structure where the polarization fraction increases in the regions in which the total intensity shows a bending morphology; in addition, the magnetic field orientation follows the direction of the bendings. We interpret these features as possible indications of an interaction with the intracluster medium. This strengthens the hypothesis of positive AGN feedback, as inferred in previous studies of this object on the basis of X-ray/mm/radio analysis. In this scenario, the proto-BCG heats the surrounding medium and possibly enhances the SFR in nearby galaxies. △ Less

Submitted 13 December, 2021; v1 submitted 7 December, 2021; originally announced December 2021.

Comments: 13 pages, 3 figures, 1 table. Accepted for publication in Galaxies

XXL-HSC: An updated catalogue of high-redshift ($z\geqslant3.5$) X-ray AGN in the XMM-XXL northern field. Constraints on the bright end of the soft logN-logS

Authors: E. Pouliasis, I. Georgantopoulos, A. Ruiz, R. Gilli, E. Koulouridis, M. Akiyama, Y. Ueda, L. Chiappetti, C. Garrel, C. Horellou, T. Nagao, S. Paltani, M. Pierre, Y. Toba, C. Vignali

Abstract: X-rays offer a reliable method to identify Active Galactic Nuclei (AGN). However, in the high-redshift Universe X-ray AGN are poorly sampled due to their relatively low space density and the small areas covered by X-ray surveys. In addition to wide-area X-ray surveys, it is important to have deep optical data in order to locate the optical counterparts and determine their redshifts. In this work,… ▽ More X-rays offer a reliable method to identify Active Galactic Nuclei (AGN). However, in the high-redshift Universe X-ray AGN are poorly sampled due to their relatively low space density and the small areas covered by X-ray surveys. In addition to wide-area X-ray surveys, it is important to have deep optical data in order to locate the optical counterparts and determine their redshifts. In this work, we build a high-redshift ($z\geqslant3.5$) X-ray selected AGN sample in the XMM-XXL northern field using the most updated [0.5-2 keV] catalogue along with a plethora of new spectroscopic and multi-wavelength catalogues, including the deep optical Subaru Hyper Suprime-Cam (HSC) data reaching magnitude limits i$\sim$26 mag. We select all the spectroscopically confirmed AGN and complement this sample with high-redshift candidates that are HSC g- and r-band dropouts. To confirm the dropouts, we derive their photometric redshifts using spectral energy distribution techniques. We end up with a sample of 54 high-z sources (28 with spec-z), the largest in this field so far (almost three times larger than in previous studies), and we estimate the possible contamination and completeness. We calculate the number counts (logN-logS) in different redshift bins and compare our results with previous studies and models. We provide the strongest high-redshift AGN constraints yet at bright fluxes ($\rm f_{0.5-2~keV} > 10^{-15} erg~s^{-1}~cm^{-2}$). The samples of $z\geqslant3.5$, $z\geqslant4$ and $z\geqslant5$ are in agreement with an exponential decline model similar to that witnessed at optical wavelengths. Our work emphasizes the importance of using wide area X-ray surveys with deep optical data to uncover high redshift AGN. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Comments: Accepted for publication in A&A. 16 pages, 12 figures, 3 tables

Journal ref: A&A 658, A175 (2022)

The MURALES survey. V. Jet-induced star formation in 3C 277.3 (Coma A)

Authors: A. Capetti, B. Balmaverde, C. Tadhunter, A. Marconi, G. Venturi, M. Chiaberge, R. D. Baldi, S. Baum, R. Gilli, P. Grandi, Eileen T. Meyer, G. Miley, C. O'Dea, W. Sparks, E. Torresi, G. Tremblay

Abstract: We present observations obtained with the VLT/MUSE optical integral field spectrograph of the radio source 3C277.3, located at a redshift of 0.085 and associated with the galaxy Coma A. An emission line region fully enshrouds the double-lobed radio source, which is ~60 kpc x 90 kpc in size. Based on the emission line ratios, we identified five compact knots in which the gas ionization is powered b… ▽ More We present observations obtained with the VLT/MUSE optical integral field spectrograph of the radio source 3C277.3, located at a redshift of 0.085 and associated with the galaxy Coma A. An emission line region fully enshrouds the double-lobed radio source, which is ~60 kpc x 90 kpc in size. Based on the emission line ratios, we identified five compact knots in which the gas ionization is powered by young stars located as far as ~60 kpc from the host. The emission line filaments surrounding the radio emission are compatible with ionization from fast shocks (with a velocity of 350-500 km/s), but a contribution from star formation occurring at the edges of the radio source is likely. Coma A might be a unique example in the local Universe in which the expanding outflow triggers star formation throughout the whole radio source. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: Pre-proofs version - Accepted for publication in A&A

Redshift identification of X-ray selected active galactic nuclei in the J1030 field: searching for large-scale structures and high-redshift sources

Authors: Stefano Marchesi, Marco Mignoli, Roberto Gilli, Alessandro Peca, Micol Bolzonella, Riccardo Nanni, Marianna Annunziatella, Barbara Balmaverde, Marcella Brusa, Francesco Calura, Letizia P. Cassarà, Marco Chiaberge, Andrea Comastri, Felice Cusano, Quirino D'Amato, Kazushi Iwasawa, Giorgio Lanzuisi, Danilo Marchesini, Takahiro Morishita, Isabella Prandoni, Andrea Rossi, Paolo Tozzi, Cristian Vignali, Fabio Vito, Giovanni Zamorani , et al. (1 additional authors not shown)

Abstract: We publicly release the spectroscopic and photometric redshift catalog of the sources detected with Chandra in the field of the $z$=6.3 quasar SDSS J1030+0525. This is currently the fifth deepest X-ray field, and reaches a 0.5-2 keV flux limit $f_{\rm 0.5-2}$=6$\times$10$^{-17}$ erg s$^{-1}$ cm$^{-2}$. By using two independent methods, we measure a photometric redshift for 243 objects, while 123 (… ▽ More We publicly release the spectroscopic and photometric redshift catalog of the sources detected with Chandra in the field of the $z$=6.3 quasar SDSS J1030+0525. This is currently the fifth deepest X-ray field, and reaches a 0.5-2 keV flux limit $f_{\rm 0.5-2}$=6$\times$10$^{-17}$ erg s$^{-1}$ cm$^{-2}$. By using two independent methods, we measure a photometric redshift for 243 objects, while 123 (51%) sources also have a spectroscopic redshift, 110 of which coming from an INAF-Large Binocular Telescope (LBT) Strategic Program. We use the spectroscopic redshifts to determine the quality of the photometric ones, and find it in agreement with that of other X-ray surveys which used a similar number of photometric data-points. In particular, we measure a sample normalized median absolute deviation $σ_{NMAD}$=1.48||$z_{phot}$-$z_{spec}$||/(1+$z_{spec}$)=0.065. We use these new spectroscopic and photometric redshifts to study the properties of the Chandra J1030 field. We observe several peaks in our spectroscopic redshift distribution between $z$=0.15 and $z$=1.5, and find that the sources in each peak are often distributed across the whole Chandra field of view. This evidence confirms that X-ray selected AGN can efficiently track large-scale structures over physical scales of several Mpc. Finally, we computed the Chandra J1030 $z>$3 number counts: while the spectroscopic completeness at high-redshift of our sample is limited, our results point towards a potential source excess at $z\geq$4, which we plan to either confirm or reject in the near future with dedicated spectroscopic campaigns. △ Less

Submitted 16 September, 2021; originally announced September 2021.

Comments: 21 pages, 15 figures. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 656, A117 (2021)

Compton-Thick AGN in the NuSTAR era VI: The observed Compton-thick fraction in the Local Universe

Authors: N. Torres-Albà, S. Marchesi, X. Zhao, M. Ajello, R. Silver, T. T. Ananna, M. Baloković, P. B. Boorman, A. Comastri, R. Gilli, G. Lanzuisi, K. Murphy, C. M. Urry, C. Vignali

Abstract: We present the analysis of simultaneous NuSTAR and XMM-Newton data of 8 Compton-thick (CT-) active galactic nuclei (AGN) candidates selected in the Swift-Burst Alert Telescope (BAT) 100 month survey. This work is part of an ongoing effort to find and characterize all CT-AGN in the local ($z\leq$0.05) Universe. We used two physically motivated models, MYTorus and borus02, to characterize the source… ▽ More We present the analysis of simultaneous NuSTAR and XMM-Newton data of 8 Compton-thick (CT-) active galactic nuclei (AGN) candidates selected in the Swift-Burst Alert Telescope (BAT) 100 month survey. This work is part of an ongoing effort to find and characterize all CT-AGN in the local ($z\leq$0.05) Universe. We used two physically motivated models, MYTorus and borus02, to characterize the sources in the sample, finding 5 of them to be confirmed CT-AGN. These results represent an increase of $\sim19$% over the previous NuSTAR-confirmed, BAT-selected CT-AGN at $z\leq0.05$, bringing the total number to 32. This corresponds to an observed fraction of $\sim 8$\% of all AGN within this volume-limited sample, although it increases to $20\pm5$% when limiting the sample to $z\leq0.01$. Out of a sample of 48 CT-AGN candidates, selected using BAT and soft (0.3$-$10 keV) X-ray data, only 24 are confirmed as CT-AGN with the addition of the NuSTAR data. This highlights the importance of NuSTAR when classifying local obscured AGN. We also note that most of the sources in our full sample of 48 Seyfert 2 galaxies with NuSTAR data have significantly different line-of-sight and average torus column densities, favouring a patchy torus scenario. △ Less

Submitted 1 September, 2021; originally announced September 2021.

Comments: 31 pages, 13 figures, accepted for publication in ApJ

Compton-Thick AGN in the NuSTAR era VII. A joint NuSTAR, Chandra and XMM-Newton analysis of two nearby, heavily obscured sources

Authors: Alberto Traina, Stefano Marchesi, Cristian Vignali, Núria Torres-Albà, Marco Ajello, Andrealuna Pizzetti, Ross Silver, Xiurui Zhao, Tonima Tasnim Ananna, Mislav Baloković, Peter Boorman, Poshak Gandhi, Roberto Gilli, Giorgio Lanzuisi

Abstract: We present the joint Chandra, XMM-Newton and NuSTAR analysis of two nearby Seyfert galaxies, NGC 3081 and ESO 565-G019. These are the only two having Chandra data in a larger sample of ten low redshift ($z \le 0.05$), candidates Compton-thick Active Galactic Nuclei (AGN) selected in the 15-150 keV band with Swift-BAT that were still lacking NuSTAR data. Our spectral analysis, performed using physi… ▽ More We present the joint Chandra, XMM-Newton and NuSTAR analysis of two nearby Seyfert galaxies, NGC 3081 and ESO 565-G019. These are the only two having Chandra data in a larger sample of ten low redshift ($z \le 0.05$), candidates Compton-thick Active Galactic Nuclei (AGN) selected in the 15-150 keV band with Swift-BAT that were still lacking NuSTAR data. Our spectral analysis, performed using physically-motivated models, provides an estimate of both the line-of-sight (l.o.s.) and average (N$_{H,S}$) column densities of the two torii. NGC 3081 has a Compton-thin l.o.s. column density N$_{H,z}$=[0.58-0.62] $\times 10^{24}$cm$^{-2}$, but the N$_{H,S}$, beyond the Compton-thick threshold (N$_{H,S}$=[1.41-1.78] $\times 10^{24}$cm$^{-2}$), suggests a "patchy" scenario for the distribution of the circumnuclear matter. ESO 565-G019 has both Compton-thick l.o.s. and N$_{H,S}$ column densities (N$_{H,z}>$2.31 $\times 10^{24}$cm$^{-2}$ and N$_{H,S} >$2.57 $\times 10^{24}$cm$^{-2}$, respectively). The use of physically-motivated models, coupled with the broad energy range covered by the data (0.6-70 keV and 0.6-40 keV, for NGC 3081 and ESO 565-G019, respectively) allows us to constrain the covering factor of the obscuring material, which is C$_{TOR}$=[0.63-0.82] for NGC 3081, and C$_{TOR}$=[0.39-0.65] for ESO 565-G019. △ Less

Submitted 1 September, 2021; originally announced September 2021.

Comments: 15 pages, 22 figures, Accepted for publication on ApJ

Lower-Luminosity Obscured AGN Host Galaxies are Not Predominantly in Major-Merging Systems at Cosmic Noon

Authors: Erini Lambrides, Marco Chiaberge, Timothy Heckman, Allison Kirkpatrick, Eileen T. Meyer, Andreea Petric, Kirsten Hall, Arianna Long, Duncan J. Watts, Roberto Gilli, Raymond Simons, Kirill Tchernyshyov, Vicente Rodriguez-Gomez, Fabio Vito, Alexander De La Vega, Jeffrey R. Davis, Dale D Kocevski, Colin Norman

Abstract: For over 60 years, the scientific community has studied actively growing central super-massive black holes (active galactic nuclei -- AGN) but fundamental questions on their genesis remain unanswered. Numerical simulations and theoretical arguments show that black hole growth occurs during short-lived periods ($\sim$ 10$^{7}$ -10$^{8}$ yr) of powerful accretion. Major mergers are commonly invoked… ▽ More For over 60 years, the scientific community has studied actively growing central super-massive black holes (active galactic nuclei -- AGN) but fundamental questions on their genesis remain unanswered. Numerical simulations and theoretical arguments show that black hole growth occurs during short-lived periods ($\sim$ 10$^{7}$ -10$^{8}$ yr) of powerful accretion. Major mergers are commonly invoked as the most likely dissipative process to trigger the rapid fueling of AGN. If the AGN-merger paradigm is true, we expect galaxy mergers to coincide with black hole accretion during a heavily obscured AGN phase (N$_H$ $ > 10^{23}$ cm$^{-2}$). Starting from one of the largest samples of obscured AGN at 0.5 $<$ $z$ $<$ 3.1, we select 40 non-starbursting lower-luminosity obscured AGN. We then construct a one-to-one matched redshift- and near-IR magnitude-matched non-starbursting inactive galaxy control sample. Combining deep color \textit{Hubble Space Telescope} imaging and a novel method of human classification, we test the merger-AGN paradigm prediction that heavily obscured AGN are strongly associated with galaxies undergoing a major merger. On the total sample of 80 galaxies, we estimate each individual classifier's accuracy at identifying merging galaxies/post-merging systems and isolated galaxies. We calculate the probability of each galaxy being in either a major merger or isolated system, given the accuracy of the human classifiers and the individual classifications of each galaxy. We do not find statistically significant evidence that obscured AGN at cosmic noon are predominately found in systems with evidence of significant merging/post-merging features. △ Less

Submitted 15 July, 2021; originally announced July 2021.

Comments: 19 pages, 11 figures, accepted in ApJ

The most luminous blue quasars at 3.0<z<3.3. II. CIV/X-ray emission and accretion disc physics

Authors: Elisabeta Lusso, Emanuele Nardini, Susanna Bisogni, Guido Risaliti, Roberto Gilli, Gordon T. Richards, Francesco Salvestrini, Cristian Vignali, Giada Bargiacchi, Francesca Civano, Martin Elvis, Giuseppina Fabbiano, Alessandro Marconi, Andrea Sacchi, Matilde Signorini

Abstract: We analyse the properties of the CIV broad emission line in connection with the X-ray emission of 30 bright SDSS quasars at z~3.0-3.3 with pointed XMM-Newton observations, which were selected to test the suitability of AGN as cosmological tools. In our previous work, we found that a large fraction (~25%) of the quasars in this sample are X-ray underluminous by factors of >3-10. As absorbing column… ▽ More We analyse the properties of the CIV broad emission line in connection with the X-ray emission of 30 bright SDSS quasars at z~3.0-3.3 with pointed XMM-Newton observations, which were selected to test the suitability of AGN as cosmological tools. In our previous work, we found that a large fraction (~25%) of the quasars in this sample are X-ray underluminous by factors of >3-10. As absorbing columns of >10$^{23}$ cm$^{-2}$ can be safely ruled out, their weakness is most likely intrinsic. Here we explore possible correlations between the UV and X-ray features of these sources to investigate the origin of X-ray weakness. We fit their UV SDSS spectra and analyse their CIV properties (e.g., equivalent width, EW; line peak velocity, $\upsilon_{\rm peak}$) as a function of the X-ray photon index and 2-10 keV flux. We confirm the trends of CIV $\upsilon_{\rm peak}$ and EW with UV luminosity at 2500 angstrom for both X-ray weak and X-ray normal quasars, as well as the correlation between X-ray weakness and CIV EW. In contrast to some recent work, we do not observe any clear relation between the 2-10 keV luminosity and $\upsilon_{\rm peak}$. We find a correlation between the hard X-ray flux and the integrated CIV flux for X-ray normal quasars, whilst X-ray weak quasars deviate from the main trend by more than 0.5 dex. We argue that X-ray weakness might be interpreted in a starved X-ray corona picture associated with an ongoing disc-wind phase. If the wind is ejected in the vicinity of the black hole, the extreme-UV radiation that reaches the corona will be depleted, depriving the corona of seeds photons and generating an X-ray weak quasar. Yet, at the largest UV luminosities (>10$^{47}$ erg s$^{-1}$), there will still be an ample reservoir of ionising photons that can explain the excess CIV emission observed in the X-ray weak quasars with respect to normal sources of similar X-ray luminosities. △ Less

Submitted 6 July, 2021; originally announced July 2021.

Comments: 22 pages, 15 figures (with 3 more figures in the Appendix), abstract abridged. Accepted for publication in A&A

The XMM-SERVS survey: XMM-Newton point-source catalogs for the W-CDF-S and ELAIS-S1 fields

Authors: Q. Ni, W. N. Brandt, C. -T. Chen, B. Luo, K. Nyland, G. Yang, F. Zou, J. Aird, D. M. Alexander, F. E. Bauer, M. Lacy, B. D. Lehmer, L. Mallick, M. Salvato, D. P. Schneider, P. Tozzi, I. Traulsen, M. Vaccari, C. Vignali, F. Vito, Y. Xue, M. Banerji, K. Chow, A. Comastri, A. Del Moro , et al. (8 additional authors not shown)

Abstract: We present the X-ray point-source catalogs in two of the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS) fields, W-CDF-S (4.6 deg$^2$) and ELAIS-S1 (3.2 deg$^2$), aiming to fill the gap between deep pencil-beam X-ray surveys and shallow X-ray surveys over large areas. The W-CDF-S and ELAIS-S1 regions were targeted with 2.3 Ms and 1.0 Ms of XMM-Newton observations, respectively;… ▽ More We present the X-ray point-source catalogs in two of the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS) fields, W-CDF-S (4.6 deg$^2$) and ELAIS-S1 (3.2 deg$^2$), aiming to fill the gap between deep pencil-beam X-ray surveys and shallow X-ray surveys over large areas. The W-CDF-S and ELAIS-S1 regions were targeted with 2.3 Ms and 1.0 Ms of XMM-Newton observations, respectively; 1.8 Ms and 0.9 Ms exposures remain after flare filtering. The survey in W-CDF-S has a flux limit of 1.0 $\times$ 10$^{-14}$ erg cm$^{-2}$ s$^{-1}$ over 90% of its area in the 0.5-10 keV band; 4053 sources are detected in total. The survey in ELAIS-S1 has a flux limit of 1.3 $\times$ 10$^{-14}$ erg cm$^{-2}$ s$^{-1}$ over 90% of its area in the 0.5-10 keV band; 2630 sources are detected in total. Reliable optical-to-IR multiwavelength counterpart candidates are identified for $\approx$ 89% of the sources in W-CDF-S and $\approx$ 87% of the sources in ELAIS-S1. 3186 sources in W-CDF-S and 1985 sources in ELAIS-S1 are classified as AGNs. We also provide photometric redshifts for X-ray sources; $\approx$ 84% of the 3319/2001 sources in W-CDF-S/ELAIS-S1 with optical-to-NIR forced photometry available have either spectroscopic redshifts or high-quality photometric redshifts. The completion of the XMM-Newton observations in the W-CDF-S and ELAIS-S1 fields marks the end of the XMM-SERVS survey data gathering. The $\approx$ 12,000 point-like X-ray sources detected in the whole $\approx$ 13 deg$^2$ XMM-SERVS survey will benefit future large-sample AGN studies. △ Less

Submitted 19 June, 2021; originally announced June 2021.

Comments: 40 pages, 35 figures. Accepted for publication in ApJS. Data products available at: https://personal.psu.edu/wnb3/xmmservs/xmmservs.html