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.
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.
In high density environments, the gas content of galaxies is stripped, leading to a rapid quenching of their star formation activity. This dramatic environmental effect, which is not related to ...typical passive evolution, is generally not taken into account in the star formation histories (SFHs) usually assumed to perform spectral energy distribution (SED) fitting of these galaxies, yielding a poor fit of their stellar emission and, consequently, biased estimate of the star formation rate (SFR). In this work, we aim at reproducing this rapid quenching using a truncated delayed SFH that we implemented in the SED fitting code CIGALE. We show that the ratio between the instantaneous SFR and the SFR just before the quenching (rSFR) is well constrained as long as rest-frame UV data are available. This SED modeling is applied to the Herschel Reference Survey (HRS) containing isolated galaxies and sources falling in the dense environment of the Virgo cluster. The latter are Hi-deficient because of ram pressure stripping. We show that the truncated delayed SFH successfully reproduces their SED, while typical SFH assumptions fail. A good correlation is found between rSFR and Hi−def, the parameter that quantifies the gas deficiency of cluster galaxies, meaning that SED fitting results can be used to provide a tentative estimate of the gas deficiency of galaxies for which Hi observations are not available. The HRS galaxies are placed on the SFR-M∗ diagram showing that the Hi-deficient sources lie in the quiescent region, thus confirming previous studies. Using the rSFR parameter, we derive the SFR of these sources before quenching and show that they were previously on the main sequence relation. We show that the rSFR parameter is also recovered well for deeply obscured high redshift sources, as well as in the absence of IR data. SED fitting is thus a powerful tool for identifying galaxies that underwent a rapid star formation quenching.
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.
Detailed studies of the spectral energy distribution (SED) of normal galaxies have increasingly been used to understand the physical mechanism dominating their integrated emission, mainly owing to ...the availability of high quality multi-wavelength data from the UV to the far-infrared (FIR). However, systems hosting dust-enshrouded nuclear starbursts and/or an accreting supermassive black hole (an active galactic nucleus or AGN) are especially challenging to study. This is due to the complex interplay between the heating by massive stars and the AGN, the absorption and emission of radiation from dust, as well as the presence of the underlying old stellar population. We used the latest release of CIGALE, a fast state-of-the-art galaxy SED-fitting model relying on energy balance, to study the influence of an AGN in a self consistent manner in estimating both the star formation rate (SFR) and stellar mass in galaxies, as well as to calculate the contribution of the AGN to the power output of the host. Using the semi-analytical galaxy formation model galform, we created a suite of mock galaxy SEDs using realistic star formation histories (SFH). We also added an AGN of Type-1, Type-2, or intermediate-type whose contribution to the bolometric luminosity can be variable. We performed an SED-fitting of these catalogues with CIGALE, assuming three different SFHs: a single-exponentially-decreasing (1τ-dec), a double-exponentially-decreasing (2τ-dec), and a delayed SFH. Constraining the overall contribution of an AGN to the total infrared luminosity (fracAGN) is very challenging for fracAGN< 20%, with uncertainties of ~5–30% for higher fractions depending on the AGN type, while FIR and sub-mm are essential. The AGN power has an impact on the estimation of M∗ in Type-1 and intermediate-type AGNs but has no effect on galaxies hosting Type-2 AGNs. We find that in the absence of AGN emission, the best estimates of M∗ are obtained using the 2τ-dec model but at the expense of realistic ages of the stellar population. The delayed SFH model provides good estimates of M∗ and SFR, with a maximum offset of 10% as well as better estimates of the age. Our analysis shows that the under-estimation of the SFR increases with fracAGN for Type-1 systems, as well as for low contributions of an intermediate AGN type, but it is quite insensitive to the emission of Type-2 AGNs up to fracAGN ~ 45%. A lack of sampling the FIR, or sub-mm domain systematically over-estimates the SFR (<20%), independent of the contribution of the AGN. Similarly, the UV emission is critical in accurately retrieving both the M∗ for Type-1 and intermediate- type AGN and the SFR of all three AGN types. We show that the presence of AGN emission introduces a scatter to the SFR-M∗ main sequence relation derived from SED-fitting, which is driven by the uncertainties on M∗. Finally, we used our mock catalogues to test the popular IR SED-fitting code DecompIR and show that fracAGN is under-estimated but that the SFR is recovered well for Type-1 and intermediate-types of AGN. The fracAGN, SFR, and LIR estimates of Type-2 AGNs are more problematic owing to a FIR emission disagreement between predicted and observed models.
We present IRAM-30 m Telescope 12CO and 13CO observations of a sample of 55 luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) in the local universe. This sample is a subset of the Great ...Observatory All-Sky LIRG Survey (GOALS), for which we use ancillary multi-wavelength data to better understand their interstellar medium and star formation properties. Fifty-three (96%) of the galaxies are detected in 12CO, and 29 (52%) are also detected in 13CO above a 3σ level. The median full width at zero intensity (FWZI) velocity of the CO line emission is 661 km s−1, and ∼54% of the galaxies show a multi-peak CO profile. Herschel photometric data is used to construct the far-IR spectral energy distribution of each galaxy, which are fit with a modified blackbody model that allows us to derive dust temperatures and masses, and infrared luminosities. We make the assumption that the gas-to-dust mass ratio of (U)LIRGs is comparable to local spiral galaxies with a similar stellar mass (i.e., gas/dust of mergers is comparable to their progenitors) to derive a CO-to-H2 conversion factor of ⟨α⟩ = 1.8−0.8+1.3 M⊙ ⟨ α ⟩ = 1 . 8 − 0.8 + 1.3 M ⊙ $ \langle\alpha\rangle=1.8^{+1.3}_{-0.8}\,M_\odot $ (K km s−1 pc2)−1; such a value is comparable to that derived for (U)LIRGs based on dynamical mass arguments. We derive gas depletion times of 400 − 600 Myr for the (U)LIRGs, compared to the 1.3 Gyr for local spiral galaxies. Finally, we re-examine the relationship between the 12CO/13CO ratio and dust temperature, confirming a transition to elevated ratios in warmer systems.
Abstract
We present the results of a Hubble Space Telescope WFC3 near-UV and Advanced Camera for Surveys Wide Field Channel optical study into the star cluster populations of a sample of 10 luminous ...infrared galaxies (LIRGs) in the Great Observatories All-Sky LIRG Survey. Through integrated broadband photometry we have derived ages, masses, and extinctions for a total of 1027 star clusters in galaxies with
d
L
< 110 Mpc in order to avoid issues related to cluster bending. The measured cluster age distribution slope of
dN
/
d
τ
∝
τ
−
0.5
+
/
−
0.12
is steeper than what has been observed in lower-luminosity star-forming galaxies. Further, differences in the slope of the observed cluster age distribution between inner- (
dN
/
d
τ
∝
τ
−
1.07
+
/
−
0.12
) and outer-disk (
dN
/
d
τ
∝
τ
−
0.37
+
/
−
0.09
) star clusters provide evidence of mass-dependent cluster destruction in the central regions of LIRGs driven primarily by the combined effect of strong tidal shocks and encounters with massive giant molecular clouds. Excluding the nuclear ring surrounding the Seyfert 1 nucleus in NGC 7469, the derived cluster mass function (CMF;
dN
/
dM
∝
M
α
) offers marginal evidence for a truncation in the power law at
M
t
∼ 2×10
6
M
⊙
for our three most
cluster-rich
sources, which are all classified as early stage mergers. Finally, we find evidence of a flattening of the CMF slope of
dN
/
dM
∝
M
−
1.42
±
0.1
for clusters in late-stage mergers relative to early stage (
α
= −1.65 ± 0.02), which we attribute to an increase in the formation of massive clusters over the course of the interaction.
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.
We present the results from our analysis of
NuSTAR
data of the luminous infrared galaxy Mrk 266, which contains two nuclei, south-western (SW) and north-eastern (NE), which were resolved in previous
...Chandra
imaging. Combining this with the
Chandra
data, we intepret the hard X-ray spectrum obtained from a
NuSTAR
observation to result from a steeply rising flux from a Compton-thick active galactic nuclei (AGN) in the SW nucleus which is very faint in the
Chandra
band, confirming the previous claim. This hard X-ray component is dominated by reflection, and its intrinsic 2–10 keV luminosity is likely to be ∼1 × 10
43
erg s
−1
. Although it is bright in soft X-ray, only a moderately absorbed NE nucleus has a 2–10 keV luminosity of 4 × 10
41
erg s
−1
, placing it in the low-luminosity AGN class. These results have implications for understanding the detectability and duty cycles of emission from dual AGN in heavily obscured mergers.