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.
ALMA observations of the long wavelength dust continuum are used to estimate the interstellar medium (ISM) masses in a sample of 708 galaxies at z = 0.3 to 4.5 in the COSMOS field. The galaxy sample ...has known far-infrared luminosities and, hence, star formation rates (SFRs) and stellar masses ( ) from the optical-infrared spectrum fitting. The galaxies sample SFRs from the main sequence (MS) to 50 times above the MS. The derived ISM masses are used to determine the dependence of gas mass on redshift, , and specific SFR (sSFR) relative to the MS. The ISM masses increase approximately with the 0.63 power of the rate of increase in SFRs with redshift and the 0.32 power of the sSFR/sSFRMS. The SF efficiencies also increase as the 0.36 power of the SFR redshift evolution and the 0.7 power of the elevation above the MS; thus the increased activities at early epochs are driven by both increased ISM masses and SF efficiency. Using the derived ISM mass function, we estimate the accretion rates of gas required to maintain continuity of the MS evolution ( yr−1 at z > 2.5). Simple power-law dependencies are similarly derived for the gas accretion rates. We argue that the overall evolution of galaxies is driven by the rates of gas accretion. The cosmic evolution of total ISM mass is estimated and linked to the evolution of SF and active galactic nucleus activity at early epochs.
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.
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.
The C ii 158 m fine-structure line is the brightest emission line observed in local star-forming galaxies. As a major coolant of the gas-phase interstellar medium, C ii balances the heating, ...including that due to far-ultraviolet photons, which heat the gas via the photoelectric effect. However, the origin of C ii emission remains unclear because C+ can be found in multiple phases of the interstellar medium. Here we measure the fractions of C ii emission originating in the ionized and neutral gas phases of a sample of nearby galaxies. We use the N ii 205 m fine-structure line to trace the ionized medium, thereby eliminating the strong density dependence that exists in the ratio of C ii/N ii 122 m. Using the FIR C ii and N ii emission detected by the KINGFISH (Key Insights on Nearby Galaxies: a Far- Infrared Survey with Herschel) and Beyond the Peak Herschel programs, we show that 60%-80% of C ii emission originates from neutral gas. We find that the fraction of C ii originating in the neutral medium has a weak dependence on dust temperature and the surface density of star formation, and has a stronger dependence on the gas-phase metallicity. In metal-rich environments, the relatively cooler ionized gas makes substantially larger contributions to total C ii emission than at low abundance, contrary to prior expectations. Approximate calibrations of this metallicity trend are provided.
Interstellar dust and starlight are modeled for the galaxies of the project "Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel." The galaxies were observed by the Infrared Array ...Camera and the Multiband Imaging Photometer for Spitzer on Spitzer Space Telescope, and the Photodetector Array Camera and Spectrometer and the Spectral and Photometric Imaging Receiver on Herschel Space Observatory. With data from 3.6 to 500 m, dust models are strongly constrained. Using a physical dust model, for each pixel in each galaxy we estimate (1) dust surface density, (2) dust mass fraction in polycyclic aromatic hydrocarbons (PAHs), (3) distribution of starlight intensities heating the dust, (4) total infrared (IR) luminosity emitted by the dust, and (5) IR luminosity originating in subregions with high starlight intensity. The dust models successfully reproduce the observed global and resolved spectral energy distributions. With the angular resolution of Herschel, we obtain well-resolved maps (available online) for the dust properties. As in previous studies, we find the PAH fraction to be an increasing function of metallicity, with a threshold oxygen abundance Z/Z 0.1, but we find the data to be fitted best with increasing linearly with above a threshold value of 0.15(O/H) . We obtain total dust masses for each galaxy by summing the dust mass over the individual map pixels; these "resolved" dust masses are consistent with the masses inferred from a model fit to the global photometry. The global dust-to-gas ratios obtained from this study are found to correlate with galaxy metallicities. Systems with Z/Z 0.5 have most of their refractory elements locked up in dust, whereas in systems with Z/Z 0.3 most of these elements tend to remain in the gas phase. Within galaxies, we find that is suppressed in regions with unusually warm dust with . With knowledge of one long-wavelength flux density ratio (e.g., f160/f500), the minimum starlight intensity heating the dust ( ) can be estimated to within ∼50%, despite a variation in of more than two orders of magnitude. For the adopted dust model, dust masses can be estimated to within ∼0.2 dex accuracy using the f160/f500 flux ratio and the integrated dust luminosity, and to ∼0.07 dex accuracy using the 500 m luminosity alone. There are additional systematic errors arising from the choice of dust model, but these are hard to estimate. These calibrated prescriptions for estimating starlight heating intensity and dust mass may be useful for studies of high-redshift galaxies.
We present a new study of stellar mass in a sample of ~70 submillimeter-selected galaxies (SMGs) with accurate spectroscopic redshifts. We fit combinations of stellar population synthesis models and ...power laws to the galaxies' observed-frame optical through mid-IR spectral energy distributions (SEDs) to separate stellar emission from non-stellar near-IR continuum. The availability of spectroscopic redshifts significantly enhances our ability to determine unambiguously not only the mass and luminosity of SMGs, but also the presence and contribution of non-stellar emission to their SEDs. By separating the stellar emission from the non-stellar near-IR continuum, we find that ~50% of our sample have non-stellar contributions of less than 10% in rest-frame H band and ~10% of our sample have non-stellar contributions greater than 50%. We find that the K-band luminosity of the non-stellar continuum emission is correlated with hard X-ray luminosity, indicating an active galactic nucleus (AGN) origin of the emission. Upon subtracting this AGN-contributed continuum component from all of the galaxies in our sample, we determine a lower median stellar mass for SMGs than previous studies, ~7 X 1010 M . We use constraints of the starburst timescale from molecular gas studies to estimate the amount of fading our sample would undergo if they passively evolve after the starburst terminates. The results suggest that typical SMGs, while among the most massive galaxies at z ~ 2, are likely to produce descendants of similar mass and luminosity to L* galaxies in the local universe.
We present HST narrowband near-infrared imaging of Pa and Paβ emission of 48 local luminous infrared galaxies (LIRGs) from the Great Observatories All-Sky LIRG Survey. These data allow us to measure ...the properties of 810 spatially resolved star-forming regions (59 nuclei and 751 extranuclear clumps) and directly compare their properties to those found in both local and high-redshift star-forming galaxies. We find that in LIRGs the star-forming clumps have radii ranging from ∼90 to 900 pc and star formation rates (SFRs) of ∼1 × 10−3 to 10 M yr−1, with median values for extranuclear clumps of 170 pc and 0.03 M yr−1. The detected star-forming clumps are young, with a median stellar age of 8.7 Myr, and have a median stellar mass of 5 × 105 M . The SFRs span the range of those found in normal local star-forming galaxies to those found in high-redshift star-forming galaxies at z = 1-3. The luminosity function of the LIRG clumps has a flatter slope than found in lower-luminosity, star-forming galaxies, indicating a relative excess of luminous star-forming clumps. In order to predict the possible range of star-forming histories and gas fractions, we compare the star-forming clumps to those measured in the MassiveFIRE high-resolution cosmological simulation. The star-forming clumps in MassiveFIRE cover the same range of SFRs and sizes found in the local LIRGs and have total gas fractions that extend from 10% to 90%. If local LIRGs are similar to these simulated galaxies, we expect that future observations with ALMA will find a large range of gas fractions, and corresponding star formation efficiencies, among the star-forming clumps in LIRGs.
We present models of deeply buried ultraluminous infrared galaxy (ULIRG) spectral energy distributions (SEDs) and use them to construct a three-dimensional diagram for diagnosing the nature of ...observed ULIRGs. Our goal is to construct a suite of SEDs for a very simple model ULIRG structure, and to explore how well this simple model can (by itself) explain the full range of observed ULIRG properties. We use our diagnostic to analyze archival Spitzer Space Telescope Infrared Spectrograph data of ULIRGs and find that: (1) in general, our model does provide a comprehensive explanation of the distribution of mid-IR ULIRG properties; (2) >75% (in some cases 100%) of the bolometric luminosities of the most deeply buried ULIRGs must be powered by a dust-enshrouded active galactic nucleus; (3) an unobscured "keyhole" view through 10% of the obscuring medium surrounding a deeply buried ULIRG is sufficient to make it appear nearly unobscured in the mid-IR; (4) the observed absence of deeply buried ULIRGs with large polycyclic aromatic hydrocarbon (PAH) equivalent widths is naturally explained by our models, showing that deep absorption features are "filled-in" by small quantities of foreground unobscured PAH emission (e.g., from the host galaxy disk) at the level of ∼1% the bolometric nuclear luminosity. The modeling and analysis we present will also serve as a powerful tool for interpreting the high angular resolution spectra of high-redshift sources to be obtained with the James Webb Space Telescope.