ABSTRACT
Code Investigating GALaxy Emission (cigale) is a powerful multiwavelength spectral energy distribution (SED) fitting code for extragalactic studies. However, the current version of cigale is ...not able to fit X-ray data, which often provide unique insights into active galactic nucleus (AGN) intrinsic power. We develop a new X-ray module for cigale, allowing it to fit SEDs from the X-ray to infrared (IR). We also improve the AGN fitting of cigale from UV-to-IR wavelengths. We implement a modern clumpy two-phase torus model, skirtor. To account for moderately extincted type 1 AGNs, we implement polar-dust extinction. We publicly release the source code (named ‘x-cigale’). We test x-cigale with X-ray detected AGNs in SDSS, COSMOS, and AKARI-NEP. The fitting quality (as indicated by reduced χ2) is good in general, indicating that x-cigale is capable of modelling the observed SED from X-ray to IR. We discuss constrainability and degeneracy of model parameters in the fitting of AKARI-NEP, for which excellent mid-IR photometric coverage is available. We also test fitting a sample of AKARI-NEP galaxies for which only X-ray upper limits are available from Chandra observations, and find that the upper limit can effectively constrain the AGN SED contribution for some systems. Finally, using x-cigale, we assess the ability of Athena to constrain the AGN activity in future extragalactic studies.
Aims.
Dust is observed in the polar regions of nearby active galactic nuclei (AGN) and it is known to contribute substantially to their mid-IR emission and to the obscuration of their UV to optical ...emission. We aim to carry out a statistical test to check whether this component is a common feature based on an analysis of the integrated spectral energy distributions of these composite sources.
Methods.
We selected a sample of 1275 broad-line AGN in the XMM-XXL field, with optical to infrared photometric data. These AGN are seen along their polar direction and we expect a maximal impact of dust located around the poles when it is present. We used X-CIGALE, which introduces a dust component to account for obscuration along the polar directions, modeled as a foreground screen, and an extinction curve that is chosen as it steepens significantly at short wavelengths or is much grayer. By comparing the results of different fits, we are able to define subsamples of sources with positive statistical evidence in favor of or against polar obscuration (if present) and described using the gray or steep extinction curve.
Results.
We find a similar fraction of sources with positive evidence for and against polar dust. Applying statistical corrections, we estimate that half of our sample could contain polar dust and among them, 60% exhibit a steep extinction curve and 40% a flat extinction curve; although these latter percentages are found to depend on the adopted extinction curves. The obscuration in the
V
-band is not found to correlate with the X-ray column density, while
A
V
/
N
H
ratios span a large range of values and higher dust temperatures are found with the flat, rather than with the steep extinction curve. Ignoring this polar dust component in the fit of the spectral energy distribution of these composite systems leads to an overestimation of the stellar contribution. A single fit with a polar dust component described with an SMC extinction curve efficiently overcomes this issue but it fails at identifying all the AGN with polar dust obscuration.
Fundamental differences in the radio properties of red quasars (QSOs), as compared to blue QSOs, have been recently discovered, positioning them as a potential key population in the evolution of ...galaxies and black holes across cosmic time. To elucidate the nature of these objects, we exploited a rich compilation of broad-band photometry and spectroscopic data to model their spectral energy distributions (SEDs) from the ultraviolet to the far-infrared and characterise their emission-line properties. Following a systematic comparison approach, we characterise the properties of the QSO accretion, obscuration, and host galaxies in a sample of ∼1800 QSOs at 0.2 <
z
< 2.5, classified into red and control QSOs and matched in redshift and luminosity. We find no strong differences in the average multiwavelength SEDs of red and control QSOs, other than the reddening of the accretion disk expected by the colour selection. Additionally, no clear link can be recognised between the reddening of QSOs and the interstellar medium as well as star formation properties of their host galaxies. Our modelling of the infrared emission using dusty torus models suggests that the dust distributions and covering factors in red QSOs are strikingly similar to those of the control sample, inferring that the reddening is not related to the torus and orientation effects. Interestingly, we detect a significant excess of infrared emission at rest-frame 2−5 μm, which shows a direct correlation with optical reddening. To explain its origin, we investigated the presence of outflow signatures in the QSO spectra, discovering a higher incidence of broad O
III
wings and high C
IV
velocity shifts (> 1000 km s
−1
) in red QSOs as compared to the control sample. We find that red QSOs that exhibit evidence for high-velocity wind components present a stronger signature of the infrared excess, suggesting a causal connection between QSO reddening and the presence of hot dust distributions in QSO winds. We propose that dusty winds at nuclear scales are potentially the physical ingredient responsible for the optical colours in red QSOs, as well as a key parameter for the regulation of accretion material in the nucleus.
ABSTRACT
There is a severe tension between the observed star formation rate (SFR)–stellar mass (M⋆) relations reported by different authors at z = 1–4. In addition, the observations have not been ...successfully reproduced by state-of-the-art cosmological simulations that tend to predict a factor of 2–4 smaller SFRs at a fixed M⋆. We examine the evolution of the SFR–M⋆ relation of z = 1–4 galaxies using the skirt simulated spectral energy distributions of galaxies sampled from the Evolution and Assembly of GaLaxies and their Environments simulations. We derive SFRs and stellar masses by mimicking different observational techniques. We find that the tension between observed and simulated SFR–M⋆ relations is largely alleviated if similar methods are used to infer the galaxy properties. We find that relations relying on infrared wavelengths (e.g. 24 ${\rm \, \mu m}$, MIPS – 24, 70, and 160 ${\rm \, \mu m}$ or SPIRE – 250, 350, and 500 ${\rm \, \mu m}$) have SFRs that exceed the intrinsic relation by 0.5 dex. Relations that rely on the spectral energy distribution fitting technique underpredict the SFRs at a fixed stellar mass by −0.5 dex at z ∼ 4 but overpredict the measurements by 0.3 dex at z ∼ 1. Relations relying on dust-corrected rest-frame ultraviolet luminosities, are flatter since they overpredict/underpredict SFRs for low/high star-forming objects and yield deviations from the intrinsic relation from 0.10 to −0.13 dex at z ∼ 4. We suggest that the severe tension between different observational studies can be broadly explained by the fact that different groups employ different techniques to infer their SFRs.
Abstract
We study the relation between obscuration and supermassive black hole (SMBH) accretion using a large sample of hard X-ray selected active galactic nuclei (AGNs). We find a strong decrease in ...the fraction of obscured sources above the Eddington limit for dusty gas (
log
λ
Edd
≳
−
2
) confirming earlier results, and consistent with the radiation-regulated unification model. This also explains the difference in the Eddington ratio distribution functions (ERDFs) of type 1 and type 2 AGNs obtained by a recent study. The break in the ERDF of nearby AGNs is at
log
λ
Edd
*
=
−
1.34
±
0.07
. This corresponds to the
λ
Edd
where AGNs transition from having most of their sky covered by obscuring material to being mostly devoid of absorbing material. A similar trend is observed for the luminosity function, which implies that most of the SMBH growth in the local universe happens when the AGN is covered by a large reservoir of gas and dust. These results could be explained with a radiation-regulated growth model, in which AGNs move in the
N
H
–
λ
Edd
plane during their life cycle. The growth episode starts with the AGN mostly unobscured and accreting at low
λ
Edd
. As the SMBH is further fueled,
λ
Edd
,
N
H
and the covering factor increase, leading the AGN to be preferentially observed as obscured. Once
λ
Edd
reaches the Eddington limit for dusty gas, the covering factor and
N
H
rapidly decrease, leading the AGN to be typically observed as unobscured. As the remaining fuel is depleted, the SMBH goes back into a quiescent phase.
Polarization is an important tool to further the understanding of interstellar dust and the sources behind it. In this paper we describe our implementation of polarization that is due to scattering ...of light by spherical grains and electrons in the dust Monte Carlo radiative transfer code SKIRT. In contrast to the implementations of other Monte Carlo radiative transfer codes, ours uses co-moving reference frames that rely solely on the scattering processes. It fully supports the peel-off mechanism that is crucial for the efficient calculation of images in 3D Monte Carlo codes. We develop reproducible test cases that push the limits of our code. The results of our program are validated by comparison with analytically calculated solutions. Additionally, we compare results of our code to previously published results. We apply our method to models of dusty spiral galaxies at near-infrared and optical wavelengths. We calculate polarization degree maps and show them to contain signatures that trace characteristics of the dust arms independent of the inclination or rotation of the galaxy.
ABSTRACT
We present a detailed study of ionized outflows in a large sample of ∼650 hard X-ray-detected active galactic neuclei (AGNs). Using optical spectroscopy from the BAT AGN Spectroscopic Survey ...(BASS), we are able to reveal the faint wings of the O iii emission lines associated with outflows covering, for the first time, an unexplored range of low AGN bolometric luminosity at low redshift (z ∼0.05). We test if and how the incidence and velocity of ionized outflow is related to AGN physical parameters: black hole mass ($\rm \mathit{ M}_{BH}$), gas column density ($\rm \mathit{ N}_{H}$), Eddington ratio ($\rm \lambda _{Edd}$), O iii, X-ray, and bolometric luminosities. We find a higher occurrence of ionized outflows in type 1.9 (55 per cent) and type 1 AGNs (46 per cent) with respect to type 2 AGNs (24 per cent). While outflows in type 2 AGNs are evenly balanced between blue and red velocity offsets with respect to the O iii narrow component, they are almost exclusively blueshifted in type 1 and type 1.9 AGNs. We observe a significant dependence between the outflow occurrence and accretion rate, which becomes relevant at high Eddington ratios log($\rm \lambda _{Edd}$) ≳ −1.7. We interpret such behaviour in the framework of covering factor-Eddington ratio dependence. We do not find strong trends of the outflow maximum velocity with AGN physical parameters, as an increase with bolometric luminosity can be only identified when including samples of AGNs at high luminosity and high redshift taken from literature.
ABSTRACT
We present the ionized gas outflow morphology in the Circinus galaxy using the narrow-field mode (NFM) of the Multi Unit Spectroscopic Explorer (MUSE) instrument onboard the Very Large ...Telescope (VLT). The NFM observations provide a spatial resolution of ∼0.1 arcsec, corresponding to a physical scale of ∼2 pc, one of the highest spatial resolution achievable using ground-based adaptive optics-assisted observations in the optical wavelengths. The MUSE observations reveal a collimated clumpy outflow profile originating near the active galactic nucleus (AGN) location and extending up to 1.5 arcsec (∼30 pc) in the north-west direction. The collimated structure then fragments into two filaments, giving the entire outflowing gas a ‘tuning-fork’ morphology. These structures remain undetected in the lower spatial resolution MUSE wide-field mode data. We explain the origin of this tuning-fork structure to the interaction of the outflow with a dense clump in the interstellar medium (ISM) as the outflow propagates outwards. The origin of the collimated structure itself could be from jet–ISM interactions on small scales. These observations also provide evidence to the origin of the ionized gas filaments previously observed in the Circinus galaxy out to kiloparsec scales. We find instantaneous and time-averaged mass outflow rates of 10−2 and 10−4 M⊙ yr−1, respectively. Based on the star formation rate in the Circinus galaxy reported in the literature, the observed ionized outflows are not expected to regulate star formation within the ∼100 pc scales probed by the NFM data.
A crucial aspect of 3D Monte Carlo radiative transfer is the choice of the spatial grid used to partition the dusty medium. We critically investigate the use of octree grids in Monte Carlo dust ...radiative transfer, with two different octree construction algorithms (regular and barycentric subdivision) and three different octree traversal algorithms (top-down, neighbour list, and the bookkeeping method). In general, regular octree grids need higher levels of subdivision compared to the barycentric grids for a fixed maximum cell mass threshold criterion. The total number of grid cells, however, depends on the geometry of the model. Surprisingly, regular octree grid simulations turn out to be 10 to 20% more efficient in run time than the barycentric grid simulations, even for those cases where the latter contain fewer grid cells than the former. Furthermore, we find that storing neighbour lists for each cell in an octree, ordered according to decreasing overlap area, is worth the additional memory and implementation overhead: using neighbour lists can cut down the grid traversal by 20% compared to the traditional top-down method. In conclusion, the combination of a regular node subdivision and the neighbour list method results in the most efficient octree structure for Monte Carlo radiative transfer simulations.
Abstract
Accreting supermassive black holes (SMBHs) located at the centers of galaxies are typically surrounded by large quantities of gas and dust. The structure and evolution of this circumnuclear ...material can be studied at different wavelengths, from the submillimeter to the X-ray. Recent X-ray studies have shown that the covering factor of the obscuring material tends to decrease with increasing Eddington ratio, likely due to radiative feedback on dusty gas. Here we study a sample of 549 nearby (
z
≲ 0.1) hard X-ray (14–195 keV) selected nonblazar active galactic nuclei (AGN) and use the ratio between the AGN infrared and bolometric luminosity as a proxy of the covering factor. We find that, in agreement with what has been found by X-ray studies of the same sample, the covering factor decreases with increasing Eddington ratio. We also confirm previous findings that showed that obscured AGN typically have larger covering factors than unobscured sources. Finally, we find that the median covering factors of AGN located in different regions of the column density–Eddington ratio diagram are in good agreement with what would be expected from a radiation-regulated growth of SMBHs.