We present an analysis of , O iii88, N ii122, and far-infrared (FIR) fine-structure line observations obtained with Herschel/PACS, for ∼240 local luminous infrared galaxies (LIRGs) in the Great ...Observatories All-sky LIRG Survey. We find pronounced declines ("deficits") of line-to-FIR continuum emission for N ii122, , and as a function of FIR color and infrared luminosity surface density, . The median electron density of the ionized gas in LIRGs, based on the N ii122/N ii205 ratio, is = 41 cm−3. We find that the dispersion in the deficit of LIRGs is attributed to a varying fractional contribution of photodissociation regions (PDRs) to the observed emission, f( ) = / , which increases from ∼60% to ∼95% in the warmest LIRGs. The / ratio is tightly correlated with the PDR gas kinetic temperature in sources where is not optically thick or self-absorbed. For each galaxy, we derive the average PDR hydrogen density, , and intensity of the interstellar radiation field, G, in units of and find G/ ratios of ∼0.1-50 cm3, with ULIRGs populating the upper end of the distribution. There is a relation between G/ and , showing a critical break at 5 × 1010 L kpc−2. Below , G/ remains constant, 0.32 cm3, and variations in are driven by the number density of star-forming regions within a galaxy, with no change in their PDR properties. Above , G/ increases rapidly with , signaling a departure from the typical PDR conditions found in normal star-forming galaxies toward more intense/harder radiation fields and compact geometries typical of starbursting sources.
Abstract
Mergers of galaxies are thought to cause significant gas inflows to the inner parsecs, which can activate rapid accretion on to supermassive black holes (SMBHs), giving rise to active ...galactic nuclei (AGN). During a significant fraction of this process, SMBHs are predicted to be enshrouded by gas and dust. Studying 52 galactic nuclei in infrared-selected local luminous and ultraluminous infrared galaxies in different merger stages in the hard X-ray band, where radiation is less affected by absorption, we find that the amount of material around SMBHs increases during the last phases of the merger. We find that the fraction of Compton-thick (CT, N H ≥ 1024 cm− 2) AGN in late-merger galaxies is higher ($f_{\rm \,CT}=65^{+12}_{-13}{\rm per\, cent}$) than in local hard X-ray selected AGN (f CT = 27 ± 4 per cent), and that obscuration reaches its maximum when the nuclei of the two merging galaxies are at a projected distance of D12 ≃ 0.4–10.8 kpc ($f_{\rm \,CT}=77_{-17}^{+13}{\rm per\, cent}$). We also find that all AGN of our sample in late-merger galaxies have N H > 1023 cm− 2, which implies that the obscuring material covers $95^{+4}_{-8}{\rm per\, cent}$ of the X-ray source. These observations show that the material is most effectively funnelled from the galactic scale to the inner tens of parsecs during the late stages of galaxy mergers, and that the close environment of SMBHs in advanced mergers is richer in gas and dust with respect to that of SMBHs in isolated galaxies, and cannot be explained by the classical AGN unification model in which the torus is responsible for the obscuration.
Enhanced emission from the dense gas tracer HCN (relative to HCO+) has been proposed as a signature of active galactic nuclei (AGN). In a previous single-dish millimeter line survey we identified ...galaxies with HCN/HCO+ (1-0) intensity ratios consistent with those of many AGN but whose mid-infrared spectral diagnostics are consistent with little to no ( 15%) contribution of an AGN to the bolometric luminosity. To search for putative heavily obscured AGN, we present and analyze NuSTAR hard X-ray (3-79 keV) observations of four such galaxies from the Great Observatories All-sky LIRG Survey. We find no X-ray evidence for AGN in three of the systems and place strong upper limits on the energetic contribution of any heavily obscured ( ) AGN to their bolometric luminosity. The upper limits on the X-ray flux are presently an order of magnitude below what XDR-driven chemistry models predict are necessary to drive HCN enhancements. In a fourth system we find a hard X-ray excess consistent with the presence of an AGN, but contributing only ∼3% of the bolometric luminosity. It is also unclear if the AGN is spatially associated with the HCN enhancement. We further explore the relationship between HCN/HCO+ (for several Jupper levels) and / for a larger sample of systems in the literature. We find no evidence for correlations between the line ratios and the AGN fraction derived from X-rays, indicating that HCN/HCO+ intensity ratios are not driven by the energetic dominance of AGN, nor are they reliable indicators of ongoing supermassive black hole accretion.
ABSTRACT
The merger of two or more galaxies can enhance the inflow of material from galactic scales into the close environments of active galactic nuclei (AGNs), obscuring and feeding the ...supermassive black hole (SMBH). Both recent simulations and observations of AGN in mergers have confirmed that mergers are related to strong nuclear obscuration. However, it is still unclear how AGN obscuration evolves in the last phases of the merger process. We study a sample of 60 luminous and ultra-luminous IR galaxies (U/LIRGs) from the GOALS sample observed by NuSTAR. We find that the fraction of AGNs that are Compton thick (CT; $N_{\rm H}\ge 10^{24}\rm \, cm^{-2}$) peaks at $74_{-19}^{+14}{{\ \rm per\ cent}}$ at a late merger stage, prior to coalescence, when the nuclei have projected separations (dsep) of 0.4–6 kpc. A similar peak is also observed in the median NH $(1.6\pm 0.5)\times 10^{24}\rm \, cm^{-2}$. The vast majority ($85^{+7}_{-9}{{\ \rm per\ cent}}$) of the AGNs in the final merger stages (dsep ≲ 10 kpc) are heavily obscured ($N_{\rm H}\ge 10^{23}\rm \, cm^{-2}$), and the median NH of the accreting SMBHs in our sample is systematically higher than that of local hard X-ray-selected AGN, regardless of the merger stage. This implies that these objects have very obscured nuclear environments, with the $N_{\rm H}\ge 10^{23}\rm \, cm^{-2}$ gas almost completely covering the AGN in late mergers. CT AGNs tend to have systematically higher absorption-corrected X-ray luminosities than less obscured sources. This could either be due to an evolutionary effect, with more obscured sources accreting more rapidly because they have more gas available in their surroundings, or to a selection bias. The latter scenario would imply that we are still missing a large fraction of heavily obscured, lower luminosity ($L_{2-10}\lesssim 10^{43}\rm \, erg\, s^{-1}$) AGNs in U/LIRGs.
ABSTRACT We report the detection of a heavily obscured active galactic nucleus (AGN) in the luminous infrared galaxy (LIRG) NGC 6286 identified in a 17.5 ks Nuclear Spectroscopic Telescope Array ...observation. The source is in an early merging stage and was targeted as part of our ongoing NuSTAR campaign observing local luminous and ultra-luminous infrared galaxies in different merger stages. NGC 6286 is clearly detected above 10 keV and by including the quasi-simultaneous Swift/XRT and archival XMM-Newton and Chandra data, we find that the source is heavily obscured (NH (0.95−1.32) × 1024 cm−2) with a column density consistent with being Compton-thick (CT, ). The AGN in NGC 6286 has a low absorption-corrected luminosity (L2−10 keV ∼ 3−20 × 1041 erg s−1) and contributes 1% to the energetics of the system. Because of its low luminosity, previous observations carried out in the soft X-ray band (<10 keV) and in the infrared did not notice the presence of a buried AGN. NGC 6286 has multiwavelength characteristics typical of objects with the same infrared luminosity and in the same merger stage, which might imply that there is a significant population of obscured low-luminosity AGNs in LIRGs that can only be detected by sensitive hard X-ray observations.
Abstract
We use a suite of cosmological zoom galaxy formation simulations and dust radiative transfer calculations to explore the use of the monochromatic 850
μ
m luminosity (
L
ν
,850
) as a ...molecular gas mass (
M
mol
) estimator in galaxies between 0 <
z
< 9.5 for a broad range of masses. For our fiducial simulations, where we assume that the dust mass is linearly related to the metal mass, we find that empirical
L
ν
,850
–
M
mol
calibrations accurately recover the molecular gas mass of our model galaxies and that the
L
ν
,850
-dependent calibration is preferred. We argue that the major driver of scatter in the
L
ν
,850
–
M
mol
relation arises from variations in the molecular gas-to-dust mass ratio, rather than variations in the dust temperature, in agreement with the previous study of Liang et al. Emulating a realistic measurement strategy with ALMA observing bands that are dependent on the source redshift, we find that estimating
S
ν
,850
from continuum emission at a different frequency contributes 10%–20% scatter to the
L
ν
,850
–
M
mol
relation. This additional scatter arises from a combination of mismatches in assumed
T
dust
and
β
values, as well as the fact that the SEDs are not single-temperature blackbodies. However, this observationally induced scatter is a subdominant source of uncertainty. Finally, we explore the impact of a dust prescription in which the dust-to-metals ratio varies with metallicity. Though the resulting mean dust temperatures are ∼50% higher, the dust mass is significantly decreased for low-metallicity halos. As a result, the observationally calibrated
L
ν
,850
–
M
mol
relation holds for massive galaxies, independent of the dust model, but below
L
ν
,850
≲ 10
28
erg s
−1
(metallicities
) we expect that galaxies may deviate from literature observational calibrations by ≳0.5 dex.
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.
Modeling of the spectral line energy distribution (SLED) of the CO molecule can reveal the physical conditions (temperature and density) of molecular gas in Galactic clouds and other galaxies. ...Recently, the Herschel Space Observatory and ALMA have offered, for the first time, a comprehensive view of the rotational J = 4−3 through J = 13−12 lines, which arise from a complex, diverse range of physical conditions that must be simplified to one, two, or three components when modeled. Here we investigate the recoverability of physical conditions from SLEDs produced by galaxy evolution simulations containing a large dynamical range in physical properties. These simulated SLEDs were generally fit well by one component of gas whose properties largely resemble or slightly underestimate the luminosity-weighted properties of the simulations when clumping due to nonthermal velocity dispersion is taken into account. If only modeling the first three rotational lines, the median values of the marginalized parameter distributions better represent the luminosity-weighted properties of the simulations, but the uncertainties in the fitted parameters are nearly an order of magnitude, compared to approximately 0.2 dex in the "best-case" scenario of a fully sampled SLED through J = 10−9. This study demonstrates that while common CO SLED modeling techniques cannot reveal the underlying complexities of the molecular gas, they can distinguish bulk luminosity-weighted properties that vary with star formation surface densities and galaxy evolution, if a sufficient number of lines are detected and modeled.
ABSTRACT
Accreting supermassive black holes (SMBHs), also known as active galactic nuclei (AGN), are generally surrounded by large amounts of gas and dust. This surrounding material reprocesses the ...primary X-ray emission produced close to the SMBH and gives rise to several components in the broadband X-ray spectra of AGN, including a power-law possibly associated with Thomson-scattered radiation. In this work, we study the properties of this scattered component for a sample of 386 hard-X-ray-selected, nearby ($z\sim0.03$) obscured AGN from the 70-month Swift/BAT catalogue. We investigate how the fraction of Thomson-scattered radiation correlates with different physical properties of AGN, such as line-of-sight column density, X-ray luminosity, black hole mass, and Eddington ratio. We find a significant negative correlation between the scattering fraction and the column density. Based on a large number of spectral simulations, we exclude the possibility that this anticorrelation is due to degeneracies between the parameters. The negative correlation also persists when considering different ranges of luminosity, black hole mass, and Eddington ratio. We discuss how this correlation might be either due to the angle dependence of the Thomson cross-section or to more obscured sources having a higher covering factor of the torus. We also find a positive correlation between the scattering fraction and the ratio of O iii λ5007 to X-ray luminosity. This result is consistent with previous studies and suggests that the Thomson-scattered component is associated with the narrow-line region.